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SECTION 1. 


GENERAL 


1.1 Purpose of the Program Maintenance Manual . 

The objective for writing this Program Maintenance Manual for project 
Nickel-Cadmium Battery Expert System (NICBES), Contract Number : 
NAS8-35922, is to provide the maintenance programmer personnel with the 
information necessary to effectively maintain or enhance the system. 


1.2 Project References . 

NICBES is an Expert System for 
Nickel-Cadmium Batteries found in 
Electrical Power System (EPS) Tes 
Center (MSFC). NICBES resides on 
two modes. The first mode is 

The second mode is 
a logical programming 
manuals serve as reference materials 
NICBES User’s Manual - September 
IBM-PC AT Manuals 

ARITY PROLOG Manuals - Version 4. 
MICROSOFT C Manuals 


MICROSOFT C. 
ARITY PROLOG, 


fault diagnosis and 
the Hubble Space 
tbed located at Marsh 
a dedicated IBM-PC AT 
the Data-Handler whi 
the Expert System whi 
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1.3 Terms and Abbreviations. 


AHI 

- 

ampere hours in 

AHO 

- 

ampere hours out 

bprc 

- 

battery protection and reconditioning 

CCC 

- 

charge current controllers 

DOD 

- 

depth of discharge 

EOC 

- 

end of charge 

EOD 

- 

end of discharge 

EOF 

- 

end of file 

EPS 

- 

electrical power system 

HS T 

- 

Hubble Space Telescope 

MMDA 

- 

Martin Marietta Denver Aerospace 

MSFC 

- 

Marshall Space Flight Center 

NICBES 

- 

Nickel-Cadmium Battery Expert System 

SOW 

- 

Statement of Work 

SPA 

- 

solar panel array 
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SECTION 2. 


SYSTEM 


DESCRIPTION 


2.1 System Application . 

NICBES was developed as an assistant for engineers working on the 
Testbed to aid in decision making with regard to the Nickel- 
Batteries. NICBES analysis depends on the particular 

configuration at MSFC, see Figure 1, and the particular 

manufacturer. 


HST EPS 
Cadmium 
tes tbed 
Battery 


2 . 2 Securi tv . 

There are no security requirements. 


2.3 General Description . 

NICBES, as programmed for the IBM-PC AT, a single tasking computer, 
requires two independent processes. The first is the Data-Handler which 
processes incoming telemetry every one minute. Input to the Data-Handler 
comes from the DEC LSI-11 over a RS232 to the IBM-PC AT. Each telemetry 
burst contains 370 integer and floating point values preceeded by the 
character ’A’. There are 96 minutes in one orbit. An orbit is composed 
of a discharge and charge phase. An orbit starts at the beginning of the 
discharge phase and ends at the completion of the charge phase. Once 
this process is completed (for 12 orbits total), the Expert System of 
NICBES can be run. Input to the Expert System consists of processed data 
files output by the Data-Handler. Output from the Expert System includes 
fault diagnosis, battery status and advice, plus decision support. Any 
Expert System screen displays can be routed to the STAR-SD-15 Printer for 
hardcopy. 


2.4 Program Description . 

Program Descriptions for NICBES will be given in two sets. First the 
Data-Handler will be described and then the Expert System. 


2.4.1 Data-Handler . 

The Data-Handler, written in MICROSOFT C, is installed on the IBM-PC AT 
according to the procedures listed in the MICROSOFT C Manual. All the 
following programs can be found under C:\USR . This is also where they 
should be executed. The data output files are also written to this 

directory. See Figure 1 for Data-Handler Flow Diagram. 

2.4. 1.1 datacl.bat 

datacl.bat creates the executable for the Data-Handler. It compiles 

data hdl.c, the main control routine, and then links all the other object 

files needed. The result is data hdl.exe which invokes the 

Data-Handler. The user simply types ’data_hdl’ from the DOS prompt to 
execute the Data-Handler. If changes are made to any of the ’C’ 
programs, they must be recompiled (’msc f i lename . c ; * ) and then relinked 
(’datacl*) to create a new executable. ’printf’ statements used in 
development have been left in the programs but commented out. These 
statements are can be reinstated for debugging purposes. 
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DATA-HANDLER 
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2.4. 1.2 hst.h 

hst.h contains all the header items for the Data-Handl er . It starts with 
define statements for errors, phases and data files. It also lists the 
needed system include files. Last are the structure and matrix 

definitions for global arrays. 

Define Statements: 

SUCCESS 1 

FAIL 0 

EOC 2 

EOD 3 

EVMIN 1 

DEFLT -9999.0 

SHOWF ( N) D AT N where N = 1 to 13 

CURF ( N ) DAT N where N = 14 to 16 

FAULTDAT 17 

Include files: 

stdio.h, stdlib.h, process. h and errno.h 

Global Structure for storing telemetry run, hid[0] = last run, 
hid[l] = current run. 


hid [ i ] . year 

Int 

i = 0 to 1 


hid [ i ] . day 

Int 

col 0 = last telemetry run 

hid [ i ] . hour 

Int 

col 1 = current 

telemetry run 

hid [ i ] .min 

Int 



hid [ i ] . sec 

Int 



hid[i] .orbit 

Int 



hid[ i] . phase 

Int 

0 if discharge, 

1 if charge 

hid [ i ] . day_min 

Int 

minute in charge 

phase 

hid[i] . night _m in 

Int 

minute in discharge phase 

hid[i] . batd [ j ] .battno 

Int 

j = 0 to 5, for 

6 batteries 

hid f i ] . batd [ j ] .cellv[k] 

Real 

k = 0 to 22, for 

23 cells 

hid [ i ] . batd [ j ] . cel lp [ k j 

Real 

per battery 


hid [ i ] . batd f j ] .batv 

Real 

p => pressure, v 

= > volts, 

hid f i ] . batd f j j .bate 

Real 

c => current 


hid [ i ] . batd [ j J .bprcc 

Real 

batt reconditioning current 

hidfi] . bat d [ j ] . bat emp [ k ] 

Real 

6 temp sensors per battery 

hid f i ] .batd[j] .batrecond 

Int 

int flag for reconditioning 

hid [ i ] . spac [ k j 

Real 

k = 0 to 12, for 

13 SPAs 

hid [ i j . bd [k j .'dusv 
hidfi] . bd f k J . busc 

Real 

k = 0 to 2, for 

3 busses 


The global matrices and arrays which support the data files are 
documented in Appendix A of the code listing for hst.h 
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2-4.1. 3 dat 

dat a hdl . c 

descript ion 
r ead_in i t ( ! , 


a hdl . c 

is the main driver 
of all the ’C’ routines 
finish(), sig_catch(), 


for the Data-Handler. Following is a 
in this program ( m a i n ( ) , in i t . ' , 
process!) and incomplete!) ). 


Identification: main 


Function: main is the Data-Handler driver. It controls the program 

flow and determines the times for events to happen. 


Input: No Input, although argc and argv make it possible to easily 

add inputs. 


Processing: Keeps track of Orbit and Phase. Controls program flow by 

its routine calling sequence. 

Output: Error Message written to screen if more than CILIMIT consecutive 

incomplete telemetry runs. 

Error Message if fault is detected. 

Message written to screen at the completion of each orbit. 


Local Variables: 

orbitno - Counter for number of orbits 
err - recieves return from read_data() 
phase_signal - EOC, EOD or EVMIN 
phasel - phase of previous telemetry run 
phase2 - phase of current telemetry run 

Global Variables: 

EOC, EOD, EVMIN 
FAULTDAT 
fault [ ] 

FAIL 
hid[ ] 

Interfaces: Calls init(), process!), finish!), sig_catch(), incomplete!) 

read_dat.c - read_data() 
writ_fil.c - wf ( ) . 
process. c - check_f aul t ( ) 

Error Handling: If read_data encounterd an EOF while reading the 

telemetry stream, FAIL is returned. If there are 
CILIMIT consecutive incomplete runs, Data-Handler 
is shutdown after the fault flag is set to 1. 

If a fault was detected the fault flag is set to 1 
and the Data-Handler is exited. 
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Identification:, init 


Function: Sets the communication port, initializes interrupt handler, 

and calls buffer initialization and read telemetry 
initialization routines. 

Input: No input. 

Processing: No processing. 

Outputs: Error message stating that an EOF was found while reading 
the telemetry stream and that the system is shutting down. 

Local Variables: 

err - receives return value from read_init(). 
port - communication port 

Global Variables: None 

Interfaces: Called by main. 

Calls system signal() and set_port(), read_init(), 
interrupt handler - serini() and process. c - df_init(). 

Error Handling: If read_init() returns FAIL, system will shutdown. 


Identification: read_init 

Function: Initializes telemetry reading to first full orbit. 

An orbit starts at the beginning of* the discharge phase. 

Input: No input. 

Processing: Reads the telemetry stream until night_min = 1, start of 

discharge phase. This signifies the start of an orbit. 
process() is then called to process this first set of data. 

Output: Prints message to the screen "Starting first full orbit”. 

Local Variables: err - receives return value from read_data(). 

Global Variables: 

FAIL, DEFLT , EVMIN 
hid [ ] 

Interfaces: Called by init(). 

Calls process() , incomplete () , 
read_dat.c - read_data(). 

Error Handling: If read_data() returns FAIL, control is passed back to 

init() also with a FAIL message. 
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Identification: finish 

Function: Exits Data Handler. 

Input : No input . 

Processing: Ends interrupt handler, writes data output files and calls 

exit ( ) . 

Output: No output. 

Local Variables: None. 

Global Variables: EOC, EOD 

Interfaces: Called by main(), init(), sig_catch(), incomplete ( ) or 

read_dat.c - get2(). 

Calls writ fil.c - write_f ile ( ) , system exit() 
and interrupt handler - serrst(). 

Error Handling: None. 

Identification: sig_catch 

Function: Catches ’~C’ signal input by operator to halt Data-Handler . 

Input: User signal input ,/ 'C’. 

Processing: No processing. 

Output: Writes message to screen, "Interrupt caught; exiting!" 

Local Variables: None. 

Global Variables: None 

Interfaces: Calls finish(). 


Error Handling: None. 



Identification: process 


Function: Calls the routines necessary to process the telemetry. 

Input: phase_signal - EOC, EOD or EVMIN 

orbitno - 0 to N, N is the number of orbits. 

Processing: No processing. 

Output: No output. 

Local Variables: None. 

Global Variables: None. 

Interfaces: Calls process. c - process_data( ) , df_initl(), 

writ_fil.c - write_f ile( ) . 

Error Handling: None. 


Identification: incomplete() 

Function: Checks for CILIMIT number of consecutive incomplete telemetry 

bursts . 

Input : None . 

Processing: If there are CILIMIT consecutive incomplete telemetry runs 

fault flags are set, fault.dat is written and finish() is 
called. 

Output: Error message. 

Local Variables: None. 

Global Variables: 

CILIMIT - Limit for consecutive incomplete telemetry runs 
conseq_incmp It - counter for incomplete telemetry bursts 
FAULTDAT 
f ault [ ] 

Called by main() and read_init(). 

Calls finishO and writ_fil.c - wf ( ) . 


Interfaces : 


2.4. 1.4 read dat.c 

read dat.c is the routine which reads the telemetry from the interrupt 

handler buffer. Following is a description of all the ’C’ routines in 

this file (read data(), proc sync(), proc head(), proc_bat ( ; , proc 

solarO, read_buf f er ( ) , getl() and get2() ). 


Identification: read data 


Function: Read telemetry from DEC LSI-11 over RS232 every 1 minute. 


Input: No input, however telemetry from HST EPS Testbed, 370 values 

preceded by ’A’, is utilized. 

Processing: Calls routines to read telemetry. 


Output: Returns FAIL or SUCCESS. 

Telemetry is placed in structured arrary hid[l] for later 
processing. A description of the structure of hid follows: 



Local 


Variables 


err - return value from proc_head(), proc_bat(), 
proc_so lar ( ) . 


Global Variables: 

T1 LIMIT , T2LIMIT 

buf[] - intermediate character storage array. 

count - keeps count of data values per telemetry burst. 

FAIL, SUCCESS 


Interfaces: Called by main(). 

Calls proc_sync(), proc_head(), proc_bat(), proc_solar ( ) . 

Error Handling: If proc_head(), proc_bat() or proc_solar() return a 

FAIL, reading is stopped and control is returned to the 
read_data() call in main(). 
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Ident if icat ion : 


proc_sync 


Function: Synchronize reading data to start of telemetry burst. 

Input: No input, but uses telemetry burst, first character is ’A’. 

Processing: Checks input characters until character ’A" is found. 

Also checks for a shutdown signal from the DEC LSI-11. 


Output: Can write message to screen "Sync Received!" 

Can write message to screen "Received shutdown signal from 
DEC LSI-11!" 


Local Variables: 

cc - character found in the interrupt handler’s buffer. 

Global Variables: 

f ault [ ] 

FAULTDAT 



Called by read_data(). 

Calls data_hdl.c - finish(), writ_fil.c 
get 1 ( ) and get2 ( ) . 


wf ( ) , 


Error Handling: Error Message if shutdown signal from DEC LSI-11. 


Identification: proc_head 

Function: Read header data from input buffer. 

Input: No input, but uses telemetry stream. 

Processing: Read header data from input buffer into buf[j. Use sscanf 

to put characters into integer format. 

Output: Integer header data put into global structured array hidfl]. 

Returns FAIL or SUCCESS. 

Local Variables err - return value of read_buf fer ( ) . 

Global Variables: 

FAIL, SUCCESS 
hid [ ] 
buf [] 



Called by read_data(). 
Calls read_buf f er ( ) . 


Error Handling: If read_buf f er ( ) returns a FAIL, proc_head() stops and 

returns FAIL to read_data(). 
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Identification: proc_bat 

Function: Head battery data from input buffer. 

Input: No input, but uses telemetry, 57 values for each of the 6 

batteries . 

Processing: Read battery data from input buffer into buf[j. Use sscanf 

to put characters into integer and floating point format. 


Output: Integer and real battery data put into global structured array 

hid[l]. Returns FAIL or SUCCESS. 


Local Variables err - return value from r ead_buf f er ( ) . 


Global Variables: 

FAIL, SUCCESS 
hid [ ] 
buf [ ] 



In 


terf aces : 


Called by read_data(). 
Calls read_buf f er ( ) . 


Error Handling: If read_buf f er ( ) returns a FAIL, proc_bat() stops and 

returns FAIL to read_data(). 


Identification: proc_solar 

Function: Read SPA and bus data from input buffer. 

Input: No input, but uses telemetry, 13 SPA and 6 bus values. 

Processing: Read SPA and bus data from input buffer into buf[]. Use 

sscanf to put characters into floating point format. 

Output: SPA and bus data put into global structured array hid[lj. 

Returns FAIL or SUCCESS. 

Local Variables err - return value from read_buf f er ( ) . 


Global 


Variables: 
FAIL, 
hid[] 
buf [ ] 


SUCCESS 


Interfaces 


Called by read_data() 
Calls read_buffer( ) . 


Error Handling: If read_buffer( ) returns a FAIL, proc_solar() stops and 

returns FAIL to read_data(). 
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Identification: read_buffer 

Function: Puts characters from interrupt handler buffer into buffj. 

Input: k is the actual number of telemetry values to be read. 

Processing: Get next character in interrupt handler buffer, this 

includes newlines and <CR>s. Characters are put in buf[j. 

Output: Returns found characters in buf[] to calling routine. 

Returns FAIL or SUCCESS. 

Local Variables nl_count - counts telemetry values read by counting 

newlines. 


Global Variables: 
count 
buf [ ] 

EOF 

FAIL, SUCCESS 


Interfaces : 


Called by proc_head(), proc_bat(), proc_solar(). 
Calls get 1 ( ) . 


Error Handling: 



If EOF is encountered whi le reading the telemetry 
stream, read_buf f er ( ) is stopped and FAIL is returned 
to the calling routine. 
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Identification getl() 

Function: Reads characters from interrupt handler buffer during 

telemetry burst. 

Input: None. 

Processing: Gets characters from interrupt handler buffer. Checks 

time and count to insure complete telemetry runs are read. 

Output: Returns EOF or character. 

Local Variables: 

cc - character read from interrrupt handler buffer. 
tO - start time for timer, 
tn - current time 

timer - difference between tn and tO. 

Global Variables: 
count 
EOF 

T1LIMIT - maximum time to wait for next character 

Interfaces: Called by r ead_buf f er ( ) and proc_sync(). 

Calls interrupt handler b_char(). 

Error Handling: If a character is not read within a given time limit 

it is assumed to be an incomplete telemetry burst. 
Passes EOF back to r ead_buf f er ( ) . 
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Identification get2() 

Function: Reads characters from interrupt handler buffer at the start 

of a telemetry run. 


Input : . None . 


Processing: Gets characters from interrupt handler buffer. Checks 

time to insure the beginning of telemetry is read within 
T2LIMIT time limit. 

Output: Returns character. Writes Error Message if time limit exceeded. 


Local Variables: 

cc - character read from interrrupt handler buffer. 
tO - start time for timer, 
tn - current time 

timer - difference between tn and tO. 


Global Variables: 

EOF 

T2LIMIT - maximum time to wait for next character 

fault [1 

FAULTDAT 


Interfaces 


Called by proc_sync(). 

Calls interrupt handler b_char() 


Error Handling: 


If a character is not read within a given time limit 
it is assumed that a telemetry run is missed. 

Sets fault flag to 1, writes fault.dat, writes Error 
Message "No communication in 3 minutes; exiting!", and 
exits . 
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2.4. 1.5 process. c 

process. c processes the current telemetry which has been stored in hid 

by read dat.c. Following is a description of all the ’€’ routines 

this file (df init(), df initlf), hidtohid(), move_buf f ers •; ) , move( 

process_data() , check_f aul t ( ) ). 

Identification: df_init 

Function: Start up initialization counters and globally defined 

structures, matrices and arrays. 

Input: No input, but uses globally defined structures and arrays. 

Processing: Sets counters, matrices and arrays to 0 or DEFLT. 

Output: No output. 

Local Variables: None. 

Global Variables: 


DEFLT - default value 
no_druns - number of 
no_cruns - number of 
fault [ ] 

FAULTDAT 

DATA FILE 

= -9999.0. Signifies missing data, 
discharge runs in an orbit, 
charge runs in an orbit. 

CORRESPONDING BUFFER 

showf 1 . dat 



eod_vo 1 tage ( 6 , 12 ) 

showf 2 . dat 



hc_voltage ( 6 , 12 ) 




high_buf f er ( 6 ) 

showf 3 . dat 



rc rat io ( 6,12) 




ahi ( 6 ) 




ahoo ( 6 ) 

showf4 . dat 



cv_eod_hv(6, 12) 




cv_eod_lv ( 6 , 12 ) 




cv eod av ( 6 , 12 ) 

showf5 . dat 



cv_hc_hv ( 6 , 12) 




cv_hc_lv ( 6 , 12 ) 




cv_hc_av ( 6 , 12 ) 

showf 6 . dat 



cellv_eod(6, 12 ) 

showf 7 . dat 



cel lv_hc (6,23) 

showf 8 . dat 



avgt ( 6 , 48 ) 

showf 9 . dat 



avg_temp ( 6 , 12 ) 




avg_temp_Duf f er ( 6 ) 

showf 10 . dat 



cp_eod (6,23) 




cp_eoc(6,23) 

showf 11 . dat 



time_tc(6, 12) 




trickle(6) 

showf 12 . dat 



rc_orbit (6) 




be drc(6,48) 

showf 13 . dat 



aho(6, 12) 

curf2 . dat 



batt_avg( 6) 


15 



t 


Interfaces: Called by data_hdl.c - init(). 

Error Handling: None. 


Identification: df_initl 

Function: After each telemetry burst is read, certain counters, 

structures and arrays used for processing need to be 
re- in i t i al ized . There are 3 cases - EOC, EOD or EVMIN. 

Input: phase_signal - EOC, EOD or EVMIN. Also uses globally defined 

structures, arrays and matrices. 

Processing: Sets counters, variables and arrays to 0 or DEFLT (-9999. 

Output: No output. 

Local Variables: None. 

Global Variables: 

DEFLT - default value = -9999.0. Signifies missing data. 
no_druns - number of discharge runs in an orbit. 
no_cruns - number of charge runs in an orbit. 

EOC, EOD, EVMIN 

DATA FILE CORRESPONDING BUFFER 

high_buffer ( 6) 
ahoo ( 6 ) 
ah i ( 6 ) 

cel 1 v_eod ( 6 , 12 ) 
cel lv_hc (6,23) 
avgt (6,48) 
avg_temp_buf f er ( 6 ) 
cp_eod (6,23) 
cp_eoc (6,23) 
bat t_avg ( 6 ) 

Interfaces: Called by data_hdl.c - process. 

Error Handling: None. 


showf 2 . dat 
showf 3 . dat 

showf 6 . dat 
showf 7 . dat 
showf 8 . dat 
showf 9 . dat 
showf 10 . dat 

curf 2 . dat 
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Identification: hidtohid 

Function: After each telemetry burst is read and processed, the data 

in hid[ll is put in hid[0] to prepare for next data burst 
which will be stored in hid[lj. 

Input: No input. 

Processing: Puts hid[l] column into hid[0] column. 

Output: No output. 

Local Variables: None. 

Global Variables: 
hid [ ] 

DEFLT 

Interfaces: Called by df_initl(). 

Error Handling: None. 
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Identification: move_buf fers 

Function: Prepares the globally defined matrices for movet). 

The affected matrices are those associated with data files 
containing data for 12 orbits. When more than 12 orbits 
have been processed, the arrays need to be shifted so that 
they contain the only the last 12 orbit’s data. 

There are two cases - EOC and EOD. 

Input: phase_signal - tells whether charge or discharge phase. 

Processing: Prepares matrices and then calls move(). 

Output: No output. 

Local Variables None. 

Global Variables: 

EOC, EOD 
DATA FILE 
showf 1 . dat 
showf 2 . dat 
showf 3 . dat 
showf 4 . dat 


showf 5 . dat 


showf 9 . dat 
showf 1 1 . dat 
showf 13 . dat 

Interfaces: Called by process_data( ) . 

Calls move(). 

Error Handling: None. 


CORRESPONDING BUFFER 
eod_vol tage ( 6 , 12 ) 
hc_voltage ( 6 , 12) 
rc_rat io ( 6 , 12; 
cv_eod_hv (6, 12) 
cv_eod_lv ( 6 , 12 ) 
cv_eod_av( 6,12) 
cv_hc_hv (6,12) 
cv_hc_lv (6,12) 
cv_hc_av ( 6 , 12 ) 
avg_temp( 6,12) 
time_tc(6, 12) 
aho ( 6 , 12) 
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Iden t ification: 


move 


Function: When the Data-Handler continues after 12 completed orbits, 

the 6x12 matrices must loose their first column, the remaining’ 
data must be shifted one column to the left and the next 
orbit’s data will be put in the 12th column. You hence always 
have the latest 12 orbits. 

Input: buffer[] to be moved. 

Processing: Shift columns in buffer one column to the left, dropping the 

first column. Set the 12th column to DEFLT (-9999.). 

Output: No output. 

Local Variables: None. 

Global Variables: DEFLT - default value = -9999.0. 

Interfaces: Called by move_buf f ers ( ) . 


Error 


Handling: 


None . 


Identification: process_data 

Function: 96 minutes of telemetry making up each orbit, are summarized 

mathematically in preparation for writing the data to output 
files. There are 3 cases - EOC, EOD and EVMIN. ' 


Input: phase_signal - EOC, EOC or EVMIN. 

orbitno - 0 to N, where N is the number of orbits. 

Processing: Data is prepared for showf ( n ) . dat , n = 1 to 13 

and curf(n), n = 1 to 3. Following, in Section 3.3.2 on 
data bases, a description is given of each data file’s 
functional requirements. 

Output: No output. 


Local 



Variables : 

col - 0 to 11, matches orbit to column number of matrices, 
sum - variable used to sum 6 temperature sensors per battery, 
avg - sum / 6 to give average temperature of battery per min. 
current_min - sum of night_min and day_min. 
jj - flag to do processing on even minutes, 
xl - used to find maximums. 
x2 - used to find minumums. 
x3 - used to find averages, 
x - miscellaneous variable. 
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Global Variables: 


DEFLT - default value 
no_druns - number of 
no_cruns - number of 
DCHGLIMIT - necessary 
CHGLIMIT - necessary 

DATA FILE 

= -9999.0. Signifies missing dat 
discharge runs in an orbit, 
charge runs in an orbit. 

number of discharge runs per orbi 
number of charge runs per orbit. 

CORRESPONDING BUFFER 

showf 1 . dat 

eod_voltag'e ( 6 , 12 ) 

showf 2 . dat 

hc_vo 1 tage ( 6 , 12) 
high_buf f er ( 6 ) 

showf 3 . dat 

rc_rat io ( 6,12) 
ahi ( 6 ) 
ahoo ( 6 ) 

showf4 . dat 

cv_eod_hv ( 6 , 12 ) 
cv_eod_lv ( 6 , 12) 
cv_eod_av(6, 12) 

showf 5 . dat 

cv_hc_hv ( 6 , 12 ) 
cv_hc_lv (6,12) 
cv_hc_av (6,12) 

showf 6 . dat 

cellv eod( 6,12) 

showf 7 . dat 

cel lv_hc (6,23) 

showf 8 . dat 

avgt (6,48) 

showf 9 . dat 

avg_temp ( 6 , 12 ) 
avg_temp_buf f er ( 6 5 

showf 10 . dat 

cp_eod (6,23) 
cp_eoc (6,23) 

showf 1 1 . dat 

t ime_t c ( 6 , 12 ) 
trickle ( 6 ) 

showf 12 . dat 

rc_orb i t ( 6 ) 
be drc(6,48) 

showf 13 . dat 

aho(6, 12) 

curf 2 . dat 

bat t_avg ( 6 ) 


The indented buffers are working arrays which support the 
main arrays. The data files followed by blanks require no 
processing but directly use telemetry from h i d [ 1 ] . 

In addition process_data determines if it is necessary to 
call move_buffers and keeps count of number of charge and 
discharge runs per orbit. 

Interfaces: Called by main(). 

Error Handling: Checks that no division by zero occurs. 

Checks that enough charge and discharge runs have been 
recorded to validate the processed data. Else DEFLT is 
left in the global -matrices and process() stops. 
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Identification: check_f aul t ( ) 

Function: Detects faults in telemetry. 

Input : No Input . 

Processing: There are four fault categories that are checked: 

1. Power Supplies 

a. SPA current < 5 amps during first 5 minutes of charge phase. 

b. SPA current >= 8 amps for 1-SPAs (1,3,5,7,9,11). 

SPA current >= 16 amps for 2-SPAs (2,4,6,8,10,12,13). 

c. SPA current > 5 amps during discharge phase. 

2. Batteries 

a. Cell voltage <= 0 volts for any cell in any battery. 

b. Cell voltage > 1.55 volts for any cell in any battery. 

3. Load Banks 

a. Sum of 3 bus currents > 99 amps. 

b. Load < 5 amps on any single bus during discharge phase. 

4. Temperature 

a. Average of the 6 temperature sensors > 25 C or < -10 C. 
Output: Returns FAIL or SUCCESS 

Local Variables: 

x - miscellaneous floating point number, 
sum - sum of various arrays. 

Global Variables: 

FAIL, SUCCESS 
hid [ ] 

Interfaces: Called by data_hdl.c - main(). 

Error Handling: Returns FAIL if fault found, SUCCESS if not. 
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2.4. 1.6 writ fil.c 

writ fil.c contains the 

buffers to output files, 
routines in this file 
f i 1 e ( ) ) . 


C* routines needed to write the processed data 
Following is a description of all the ’C’ 
(wff), write fl(), write_f2 ( ) , wri te_f3( ) , write 


Identification: write_file 

Function: Determines which data files should be written. There are 3 

cases - EOC, EOD or EVMIN. 


Input: phase_signal - EOC, EOD or EVMIN. 


Processing: According to the time, wf is called to write the data 

output files. Following is a time chart: 


TIME 

EOC 


EOD 


EVMIN 


DATA FILES 
showf 2 . dat 
showf 3 . dat 
showf 5 . dat 
showf 6 . dat 
showf 7 . dat 
showf 8 . dat 
showf 9 . dat 
showf 1 1 . dat 
showf 12 . dat 
showf 1 . dat 
showf 4 . dat 
showf 6 . dat 
showf 1 0 . dat 
statf 1 . dat 
curf 1 . dat 
curf 2 . dat 
curf 3 . dat 


fault.dat is written initially with fault flag = 0, and then 
only after fault flag is set to 1. 

Output: Error message is written to the screen 

"Couldn’t open ’filename’!". 

Local Variables: err - return value from wf ( ) . 

Global Variables: 

EOC, EOD, EVMIN 

SHOWF (N) DAT for N = 1 to 13 

CURF (N) DAT for N = 1 to 3 
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Interfaces: Called by data_hdl.c - process() and finish!). 

Calls wf ( ) . 

Error Handling: If a data file can not be opened, a message is written 

to the screen. No other action is taken. 


Identification: wf 

Function: Write output files for Expert System from globally defined 

matrices and arrays containing summarized telemetry. There 
is a case statement for each data file. 

Input: Name of the data file to be written. 

Processing: Open output files, write processed data from matrices in 

list format, then close output file. 

Output: Data files to be used by the Expert System. 

Returns FAIL or SUCCESS. 


Local Variables: 

err - return value from write_fl(), write_f2(), write_f3() 
sfp - output file pointer. 


Global Variables: 

SHOWF ( N ) D AT for N = 1 
CURF ( N ) D AT for N = 1 
FAIL, SUCCESS 
hi d [ ] 

DATA FILE 
showf 1 . dat 
showf 2 . dat 
showf 3 . dat 
showf4 . dat 


showf 5 . dat 


showf 6 . dat 
showf 7 . dat 
showf 8 . dat 
showf 9 . dat 
showf 10 . dat 

showf 11 . dat 
showf 12 . dat 

showf 13 . dat 
curf2 . dat 


to 13 

to 3, and FAULTDAT 


CORRESPONDING BUFFER 
e o d_v oltage(6, 12) 
hc_voltage (6,12) 
rc_rat io ( 6,12) 
cv_eod_hv (6,12) 
cv_eod_lv (6,12) 
cv_eod_av(6, 12) 
cv_hc_hv ( 6 , 12 ) 
cv_hc_lv ( 6,12) 
cv_hc_av ( 6 , 12 ) 
cel lv_eod ( 6,12) 
cel 1 v_hc (6,23) 
avgt (6,48) 
avg_temp(6, 12) 
cp_eod ( 6,23) 
cp_eoc ( 6,23) 
t ime^t c ( 6,12) 
rc_orbit ( 6) 
bc_drc (6,48) 
aho(6, 12) 
batt_avg(6) 
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Interfaces: Called by wr i t e_f i le ( ) , 

data_hdl.c - main(), incompl et e ( ) , 
read_dat.c - get2().. 

Calls write_fl(), write_f2(), write_f3(j, 

system fopen(), fclose() and fprintf(). 

Error Handling: Returns FAIL if data file can not be opened. 

Note: List format means that data is contained in brackets and separated 
by commas. A main list holds all the sublists and ends with a 
peroid. Following is an example for showfl.dat: 

show (1, [ [al,bl,cl,dl,el,fl,gl,hl,il,jl,kl,ll] , 


[ a2 , b2 , c2 , d2 , e2 , f2 , g2 , ] , 

[a3, ] , 

l a 4 * ] i 

[ a5 ] , 

[a6, ]]). 


In this example there are 
listed in chronological order. 


columns, one per orbit. Orbital data is 


Identification: write_fl 

Function: Write output files for Expert System from globally defined 

matrices and arrays containing summarized telemetry. 

Input: filename, buffer and file number to be written to output file. 

Processing: Open output files, write processed data in list format, 

then close output file. 

Output: Data files to be used by the Expert System. 

Returns FAIL or SUCCESS. 

Local Variables file - output file pointer. 

Global Variables: FAIL, SUCCESS 

Interfaces: Called by wf() for showfl, showf2, showf3, showf9, statl. 

Calls system fopen(), fclose() and fprintf(). 

Error Handling: Returns FAIL if data file can not be opened. 
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Identification: write_f2 

Function: Write output files for Expert System from globally defined 

matrices and arrays containing summarized telemetry. 

Input: filename, buffer and file number to be written to output file. 

Processing: Open output files, write processed data in list format, 

then close output file. 

Output: Data files to be used by the Expert System. 

Returns FAIL or SUCCESS. 

Local Variables: file - output file pointer. 

Global Variables: FAIL, SUCCESS 

Interfaces: Called by wf() for showf4 and.showf5. 

Calls system fopen(), fclose() and fprintf(). 

Error Handling: Returns FAIL if data file can not be opened. 

Identification: write_f3 

Function: Write output files for Expert System from globally defined 

matrices and arrays containing summarized telemetry. 

Input: filename, buffer and file number to be written to output file. 

Processing: Open output files, write processed data in list format, 

then close output file. 

Output: Data files to be used by the Expert System. 

Returns FAIL or SUCCESS. 

Local Variables: file - output file pointer. 

Global Variables: FAIL, SUCCESS 

Interfaces: Called by wf() for showfS and showf7. 

Calls system fopen(), fclose() and fprintf(). 

Error Handling: Returns FAIL if data file can not be opened. 
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2.4. 1.7 Interrupt Handler 

The Interrupt Handler takes over control of the IBM poling technique for 
receiving' data over a communication net. Instead each incoming' character 
is retrieved in a buffer which can be accessed by tne Data- Handler 
programs. This is to insure that the telemetry is read accurately and not 
written over. 

These are the programs needed: 
serial. c and serial. obj 
com_cfns.c and com_cfns.obj 
com_fns . asm and com_fns.obj 
serset.asm and serset.obj 
fixup. asm and fixup. obj 

Functions directly called from the Data-Handler are: 
data_hdl.c - init() calls serini() 

- finish() calls serrst() 
read_dat.c - g'etl() calls b_char() 
get2 ( ) calls b_char() 

Include files required are: 
ser ial . h 
entry . h 
asment . h 
asmexi t . h 
sasment . h 
sasmexit . h 
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2.4.2 Expert System . 

The Expert System is written in ARITY PROLOG which is installed on the 
IBM-PC AT according to the procedures listed in the ARITY PROLOG Manual. 
All the following programs can be found under C:\PROLOG. This is also 
where they should be executed. The set of data files to be used for the 
Expert System analysis need to be copied to C:\PROLOG. Unlike the ’C’ 
programs, the PROLOG programs of made up of many predicates and control is 
implemented by predicate calls to other predicates. In a way this is 
similar to subroutines at a smaller level. The programs have been grouped 
to be modular. See Figure 2 for Expert System Flow Diagram. 

2.4.2. 1 start ■ prg 

Identification: start 

Function: Main control routine for the Expert System. It calls other 

segments of the Expert System via the utilization of user 
menus . 

Input: fault.dat, curfl.dat and user responses to menus. 

Processing: First fault.dat is checked to see if the fault flag has been 

set to 1. If so, faultd.prg is called to perform fault 
diagnosis. The user can then opt for more information. 

In this case or if no fault, the Main Menu is written to 
the screen from which the user can select from Plots and 
Graphs, Battery Status or Advice. Next the user is asked 
to select Battery. 

Control is passed to one of the above 3 choices with the 
selected Battery. Further menus are shown for Plots and 
Graphs and for Advice. The user can always opt for another 
Battery selection or to Quit to the Main Menu where they can 
opt to Quit NICBES. 

Output: Menus and control parameters which are passed to other 

portions of the Expert System telling what the user’s choices 
are in response to menus, and which Battery to view. 

Interfaces: Invoked by prolog.ini, the PROLOG initiation program. 

Start calls functions in faultd.prg, status. prg, advice. prg, 
showpak.prg and utility. prg . 
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EXPERT SYSTEM FLOW DIAGRAM 



FIGURE 2 












2.4. 2.2 faultd.prg 

Identification: faultd 

Function: Perform fault diagnosis for the HST EPS Testbed. 

Input: curf2.dat and curf3.dat 

Processing: Five conditions are checked to determine the source of the 

fault. See faultd.doc in Appendix B. 

. Output: Output is in the form of screen report detailing the fault 

cause(s) and advising on correctional procedures. 

Interfaces: Called by start. prg - fault_diag/0 . 

2.4. 2. 3 status. prg 

Identification: status 

Function: Status analysis is performed for Batteries 1 to 6. 

Input: Battery number (Bat), showf3.dat, showf4.dat, showf8.dat, 

showf 13 . dat 

Processing: Status checks reconditioning flag first. If battery is 

being reconditioned status stops because data would be 
misleading. If not, temperature, workload, charging 
scheme and divergence are checked using averages which are 
compared to thrfeshold values. See status.doc in Appendix B. 

Output: Output is in the form of a screen report detailing the 

condition of the battery with respect to the above checks. 

Interfaces: Called by start. prg - b a t t ery_s t at us ( Bat ) . 

Calls functions in utility. prg. 



2. 4. 2. 4 advice. prg 
Identification: advice 

Function: Advice uses trend analysis for voltage, recharge ratio, 

temperature and divergence to give further detail on three 
subjects: whether a battery needs recondit ionin , changes in 

charging scheme or changes in workload. 

Input: Battery number (Bat), Advice Menu Choice (1 to 3), showfl.dat, 

showf2.dat, showf3.dat, showf5.dat, showf9.dat and showfl3.dat. 

Processing: Depending on the Choice, data files are read, trends are 

derived using the difference of two weighting functions and 
deviation factors. These trends are then compared to 
conditions to tell whether a battery needs to be changed. 
Explanations are given to back up the resulting diagnosis. 
See advice.doc in Appendix B. 

Output: Output is in the form of a screen report detailing Battery 

Advice and explanations. 

Interfaces: Called by start. prg - advice ( Bat , Choice ) . 

Calls functions in utility. prg. 


2 . 4 . 2 . 5, showpak.prg 
Identification: showpak 

Function: Decision Support portion of the Expert System providing 

12 Plots to the user for each battery. 

Input: Plot (N), Battery Number (Bat), Orbit number (Orbit), 

showft.dat (# from 1 to 12). 

Processing: Data from the appropriate data file, for the appropriate 

Battery is read. The data structure show/ll, containing 
the parameters needed for plotting, is called. 

See showpak.doc in Appendix B. 

Output: Plotting parameters are passed to grafpak.prg 

Interfaces: Called by start. prg - show_view ( N , Bat , Orbi t ) . 

Calls functions in grafpak.prg and utility. prg. 
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2. 4. 2. 6 grafpak.prg 
Identification: grafpak 

Function: Draws to the screen the any of the 12 available plots, 

per battery. 

Inputs: List of points to be plotted and all plotting' parameters 

including captions. 

Processing: Uses graphics primitives to draw plots on the screen. 

Plots have X and Y axes, title, header and points displayed 
in color and symbol. Missing data and data out of range are 
also displayed. See grafpak.doc in Appendix B for more 
details and start.doc for a listing of the graphs. 

Outputs: Plots drawn to the screen. 

Interfaces: Called by showpak.prg - graphplus/6 and plot/8. 

Calls functions in utility. prg 


2. 4. 2. 7 utility. prg 
Identification: utility 

Function: Collection of miscellaneous Prolog functions used by one or 

more of the Prolog routines. 

Inputs: Parameters are passed for the particular function call. 

Processing: Depends on the function call. See utility.doc in Appendix B 

for detailed description of the functions as well as built-in 
Arity functions. Details on the handling of data files is 
also described there. 

Outputs: Sends requested values back to the calling function. 

Interfaces: Called by start. prg, faultd.prg, status. prg, advice. prg, 

showpak.prg, grafpak. prg. 



2. 4. 2. 8 prolog.ini 
Identification: prolog.ini 

Function: Consults the programs needed to run the Expert System 

Inputs: No inputs. 

Processing: start. prg, faultd.prg, status. prg, advice. prg, showpak.prg 

grafpak.prg and utility. prg are loaded at the initiation of 
Prolog, when the user enters ’api’ at the DOS prompt. 

All data files are copied to the NICBES directory. 

The currenet data files and fault.dat are also loaded. 

The Expert System is then called into operation. 

Outputs: No output. 

Interfaces: No interfaces. 
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SECTION 3. 


ENVIRONMENT 


3.1 Equipment Environment . 

The following computer equipment is needed for the execution of NICBES: 

DEC LSI-11 

RS232 cable connector 
IBM-PC AT 

STAR-SD-15 Printer 

3.2 Software Support 

The following computer software is needed for the execution of NICBES: 

DOS (IBM’s operating system) 

MICROSOFT C 

ARITY PROLOG - Version 4.1 

3.3 Data Base 

The following paragraphs will detail the data base utilized by NICBES. 

3.3.1 General Characteristics . 

As NICBES is actually two systems, a data base description will be given 
for each. 

First for the Data-Handler whose data base consists of telemetry, received 
and processed every one minute. This dynamic data base is not stored, but 
condensed and summarized by performing mathematical operations. The final 
historical data will be written to files for use by the Expert System. 
The only limitations for the telemetry are time constraints and reading 
and writing validity. 

The data base for the Expert System consists of the data output files 
written by the Data-Handler. These files are static and should not be 
modified. However, they can be stored in uniquely referrenced locations 

for later review. 

3.3.2 Organization and Detailed Description . 

Telemetry for the Data-Handler : 

Start of Telemetry Burst 
A 

Header Information (Integer) 

1. year 

2. day of year - 198X 

3. hour ^ 0 to 24 

4. minute - 0 to 60 

5. second - 0 to 60 

6. orbit - Positive Integer 

7. phase - 0 for discharge, 1 for charge 

8. day minute (minute in charge) - 0 to 70 

9. night minute (minute in discharge) - 0 to 37 
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Battery Information (for each of 6 batteries) 
.10 - 351, 57 values for each battery 


o < 

battery number 
cell voltage 
cell pressure 
battery voltage 
battery current 


bprc current 
temperature sensors 
battery reconditioning 


Integer 1 - 
23 Reals -2 
23 Integers 
Real 
Real 


Real 


6 

to 2 volts 
0 to 150 psi 
0 to 40 volts 
-30 to +25 amps 
negative for discharge phase 
positive for charge phase 
0 to 5 amps 
6 Reals -15 to 30 (degrees C) 
Integer 0 for no, 1 for yes 


Miscellaneous Information 


352 

- 364 

Solar Array current 

13 

Reals 

0 to 

20 

amps 

365 

- 367 

Bus 

Voltage 

3 

Reals 

0 to 

40 

volts 

368 

- 370 

Bus 

current 

3 

Reals 

0 to 

90 

amps 


Reals are five place floating point numbers. Each telemetry value is 
sent one per line with an associated new line and carriage return. 370 
values are sent every one minute, 96 minutes per orbit. 


Data Files for the Expert System : 

All data files are written in list format. All the show files have 6 


lists, one for each battery. 
PROLOG Expert System as facts, 
details. 

fault.dat - Contains a fault flag = 


All the data files are loaded into the 
See documentation in Appendix B for 


1 if there was a fault 
0 if no fault was detected. 


curfl.dat - Contains the current orbit number and 

a reconditioning flag for each battery = 1 for reconditioning 

= 0 no reconditioning. 


curf2.dat - Contains Phase (charge or discharge) 

D ay_min 

Current from 13 SPAs (Solar Panel Array) 

Current from 3 Busses 

Average Temperature for 6 Batteries 

curf3.dat - Contains 6 battery cell voltages (23 per battery) 

showfl.dat - File contains battery voltage at EOD for last 12 orbits, in 
chronological order. 

showf2.dat - File contains the battery voltage during high in-charge 
period, last 12 orbits. 
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showf 3 . dat 
orbits. 
showf4 . dat 
value and 
orbits 


File contains the recharge ratio = AHO/AHI per orbit for 12 

File contains cell voltages at EOD, with the high value, low 
average of all values, in this order for each of the last 12 


showf5.dat - File contains cell voltages at high-charge; high, low and 
average of all values, order H,L,A, for each of last 12 orbits, per 
battery . 


showf6.dat - File contains 23 cell voltages at EOD for each battery, from 
the latest orbit. 


showf7.dat - File contains 23 cell voltages at high-charge for each 
battery, from latest orbit. 

showf8.dat - File contains the averages of the six temperature sensors 
(degrees C), at two minute intervals over the latest orbit. The first 
value in this file is the minute into orbit, followed by the temperature 
readings for the batteries. 

showf 9 . dat - File contains the average battery temperatures per orbit for 
the last 12 orbits. 


showfl0.dat - File contains the 23 cell pressures taken at EOC and then 

EOD for each battery in the last full orbit. 

showfll.dat - File contains the time on trickle charge for each battery 

from last 12 orbits. 

showfl2.dat - File contains battery current during reconditioning, at 

2-minute intervals, for last reconditioning of each battery. It is 

recorded every 2 minutes, only when battery reconditioning is 1 and only 
for one orbit. The file contains zeroes until a battery is reconditioned. 

showfl3.dat - File contains AHO summed at EOD over last 12 orbits. 



SECTION 4. 


PROGRAM MAINTENANCE PROCEDURES 


4.1 Conventions . 

Each routine in the programs comprising the Data-Handler have headers as 
well as code documentation. The Expert System files are documented in 
separte files having the same name as the Prolog program but with ’doc’ as 
their extension. 

a. Conventional extensions to file names are designed as mnemonic 
identifiers (file and variable names) based upon descriptive 
abbreviations of function title. 

c ’C’ programs 

obj object files 

exe executable files 

asm assembler programs 

prg Prolog programs 

dat data files 

doc document files 

b. Refer to SAMSO EX 2.3.3 and MIL-STD-847 (Documentation). 

4.2 Verification Procedures . 

Any enhancements added to the Data-Handler should be verified by checking 
the data output files. It is always wise to test the changes on a small 
test case before running the whole procedure. One enhancement that could 
be made is to increase the error checking on incoming telemetry so that 
faults could be detected directly from the Data-Handler. Another 

enhancement would be to check the ranges on each tel eme t ry value as it is 
read in. 

Changes to the Expert System logic would have to be verified by a Nickel 
Cadmium Battery ’expert’ for validity. These could include changes to the 
deviation factors, threshold variations and adding conditions to be 
checked. Enhancements to the screen displays can obviously be checked by 

running the Expert System and viewing the screen. 

The Test Procedures listed in Appendix C will give you a baseline upon 
which to verify any changes. 

4.3 Error Conditions . 

There are no special provisons for operating system errors. Procedures to 
take at such instances would include rebooting the IBM-PC AT, checking to 
see that all files are intact and starting the NICBES system again. 

4.4 Special Maintenance Procedures . 

Data files written by the Data-Handler need to be archived for Later use. 
One way of doing this would be to create a data directory at the root. 
Then for each set of data files created, a sub-directory could be created 
into which the data file set could be copied. This sub-directory can then 
be referrenced by location and time. A command file - data.bat, has been 
written for just this purpose. It is located in C:\USR. To run data.bat 
simply enter ’data directory-name <CR>’. The directory-name can be a 
date as 01-16-87 for later referencing. 
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It is also wise to make periodic backups of the NICBES system as well as 
the accumulated data. There is no need to backup the MICROSOFT C 
directories (except C:\USR) or the ARITY PROLOG files as these can always 
be re-installed from their original disks. 

4.5 Special Maintenance Programs . 

There are no special maintenance programs. 

4.6 Listings . 

All NICBES program listings will accompany this Maintenance Manual. The 
’C’ programs are documented internally while the Prolog programs are 
documented in files with the same name as the Prolog routine but with 
’doc’ as their extension. APPENDIX A contains the Data-Handler listings. 
APPENDIX B contains the Expert System listings. 
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APPENDIX A 


DOCUMENTED CODE LISTINGS FOR THE DATA-HANDLER 



ORIGINAL PAGE IS 
OF POOR QUALITY 

HST. H 


/ * D e f i n e D a t a - H a n d . 1 . e r c on s t a n ts 
if d e f i n e S U C C E S S 1 

# define • FAIL 0 


if define 

EOC 

2 


/ * P h a s e s i g ri a 1 s 

ifde f j. n e 

EOD 




ifde f i. n e 

EVMIN 

i 



ifdef ine 

DEFLT 

9999 

.0 

/+• Missing data si 

#def ine 

SH0WF1DAT 

i 


/4 Numbers for dat 

ifde f i n e 

SH0WF2DAT 




ifde f i n e 

SH0WF3DAT 

*? 



ifdef i. ne 

SH0WF4DAT 

4 



ifde fine 

SH0WF5DAT 

5 



ifde fine 

SHGWF 6 DAT 




ifdef ine 

SH0WF7DAT 

“7 



ifdef ine 

SHOWFSDAT 

8 



ifdef ine 

SH0WF9DAT 

9 



ifdef ine 

SH0WF.10DAT 

10 



ifde fine 

SH0WF.1 1DAT 

11 



ifdef ine 

SH0WF12DAT 

12 



ifdef ine 

SH0WF13DAT 

13 



ifde f i n e 

CURF1DAT 

14 



ifdef ine 

CURF2DAT 

15 



ifde f .1 n e 

CURE 3D AT 

16 



ifde f i ri e 

FAULTDAT 

17 



if i n c lude 

<stdio . h > 



/+■ Include files f 

if inc lude 

< e r r ri o . h > 




if ine lude 

< st.dl ib . h > 




if inc lude 

< process . h > 




ifi ri elude 

. <time.h> 




if ine lude 

< signal . h > 




/:*: ***** 

Structures for 

HST 

incomming data file ***** 

struct 

bati ■{ 




i. ri t 

bat tno ; 

/■* 

Battery 

Number +/ 

f 1 oa t 

cellv [23] ; 

/* 

Ce 1 1 Vo 1 tage */ 

float 

cellp [23] ; 

/* 

Cell Pressure +■/ 

float 

batv ; 

/ * 

Battery 

Voltage */ 

float 

bate ; 

/* 

battery 

Current •+/' 

f 1 oa t 

bprcc ; 

/* 

BPRC Cur 

rent +■/ 

T" 1 0 ci t 

|— *, "■ +• .—v m P X- **I > 

u. - ci i.. c: m K L J t 

/* 

Cc :Ui +- T' tTi fc- V/ 
L.- 1 i.- c- v_. » ■/ 

T e m p e r a t u r e * / 

i ri t 

batrecond ; 

/■*■ Battery 

Recoridi tion i rig*/ 


}; 
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str uct 

bus 

/ * 

Bus data 

.•+.-/ 

t 

float 

bus V ; 

/ * 

Bus Voltage 

+/ 

f 1 oa t 

bUSC : 

/ ' * '■ 

Bus Current. 


str uc t 

'! S t .1 d 


/* HST Iricomming Data 


i n t 

year ; 

/ * 

Year 

4- / 

1 n t 

day ; 

/* 

Day 

4 / 

i. n t 

hour ; 

/:+: 

Hour 

; +: / 

.1 ri t 

in i n ; 

/* 

Minutes 

r/ 

i n t 

sec ; 

/*■ 

Seconds 

*■/ 

i ri t 

o r b i t ; 

/ ■ 

Orbit# 


i n t 

phase ; 

/* 

Phase=0 if night 

* / 



/* 

=1 if day 

+/ 

i n t 

day_miri ; 

/* 

M.in into charge peroic 

:+:/ 

i ri t 

night _m in ; 

/* 

Min into discharge 

*/ 

struct 

bati batd[6] ; 

/* 

Battery Data 


float 

spac [13] ; 

/.*. 

SPA Current. 

*/ 

str uct 

bus bd[3] ; 

/* 

Bus info for three bus 

*/ 


3 i'l i d [ 2 ] ; ’ 


/ : ’i ; 

1" .1. O ct 1 


S how F i 1 e # 1 ( .1 2 orbi t s f o r 
e o d _ v o 1 1 a g e [ 6 j [12] ; 


a battery and for 6 batteries) 
/* Voltage taken at EOD for 12 
/* orbits, 6 batteries. To be 
/* updated when phase changes 
/-*• from 0 to 1 


float 
f 1 oa t 


show File#2(12 orbits for a 
h c _ v o 1 1 a g e [ 6 ] [12] ; / *■ 

h i g h _ b u f f e r [ 6 ] ; / * 

/* 

/* 


battery and for 6 batteries) 
Max Batt Voltage taken for 1 
orbits, 6 batteries during 
phase 1. To be updated when 
phase changes from 1 to 0 


/* 

float 
float 
f 1 oa t 


how File#3(12 orbits 

for a 

rc_ratio[6] [12] ; 

/* 

ahoo [6] ; 

/ :+: 

ahi [6] ; 

/* 


/ * 
• ..L. 


i *'r* 
/* 


/ r 
/* 


b a 1 1 e r y a ri d f or 6 b a 1. 1 e r i. e s ) 
Recharge Ratio for 12 orbits 
6 batteries. Sum of battery 
current during phase 1 
divided by sum of battery 
current during phase 0. It 
i s a 1 w a y s > 1 . U p d a t e d w h e n 
phase changes from 0 to 1 


ft- 2 


• 

/* 
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Show File if 4 

*/ 

float 
f 1 oa t 
f .1 oa t 

c v _ e o d hi v [ 6 ] [ 1 2 ] ; 
cv eod iv[6l [12] ; 
c v _ e o d a v [ 6 J [ 1 2 j 

/*■ Cell Voltage for 6 batteries */ 
/* Calculate high, low and avg */ 
/ *■ out of 23 cells when phase > / 
/* changes to 1 . */ 

/ * 

f 1 oa t 
f 1 oa t 
r loat 

Show Fi. le#5 
c v _ h c _ h v [ 6 ] [ 1 2 1 ; 

cv .hc lv[6] [12] ; 

cv_hc_av[6] [12] ; 

r- / 

/■*■ Calculate high, low and avg +/' 
/* out of 23 cells during phase */ 
/ * 1 , a t h i gh-char g e . T o b e + : / 
/* update when phase changes */ 
/* from 1 to 0. i/ 

/* 

float 

Show File#6 

cellv_eod[6] [23] ; 

.+ / 

/* Cell Voltage taken at EOD. ■+/ 
/ + To be updated when phase ■*'-/ 
/*■ changes from 0 to .1 . */ 

/ * 

float 

• 

Show File#? 

cellv_hc [6] [23] ; 

1/ 

/* Cell Voltage taken during +/ 
/* phase 1 at high-charge. */ 
/+ To be updated when phase */ 
/* change from .1 to 0 ■+/' 

/ * 

1" .1. oa t 

Show F i 1 e#S 
avgt [6] [48] ; 

• |: / 

/* avg temperature taken from *7 
/*■ 6 sensors per battery, 6 */ 
/* batteries, at 2 min intervals*/ 

/* 

f .1. oa t 
float 

Show File#9 

avg_temp[6] [12] ; 
avg_temp_buf fer [6] ; 

+7 

/* Every Minute take an average +/ 
/* of 6 sensors’ temp for 96 +•/ 
/* minutes in the orbit, for .12 */' 
/* orbits, & for each of 6 batts +/ 
/* Avg over 96 minutes or entire*/ 
/* orbit. To be updated when */ 
/* when phase change from 1 to 0 +/ 

/.+. 

f 1 oa t 
f .1. oa t 

Show FilettlQ 
cp_eod[6] [23] ; 
cp_eoc [6] [23] ; 

• : F ; / 

/* Cell pressure at eod *7 

/* Cell pressure at eoc *■/ 


fir 3 
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Show File# 1.1 

ti.mcy._tc [6] [ 12 j 
t r i c k 1 e i" 6 1 ; 


When battery current is lower i / 
than 3 amps, add 1 to trickle*/ 
t i me c ou n t e r . On 1 y du r i n g * / 
c h a r g e ph a s e . * / 


Show Fi le#12 
r c o r b i t 1 6 1 ; 
bc_.drc [6] [48] 


0 r b i t a t w h i c h r e c o ri d i t i. o n i ri g * / 
t a k e s p 1 a c e . B a 1 1 e r y c u r r e n t * / 
during reconditioning at •+/ 

2 minute intervals *•/ 


Stat f ile#l 
aho[6] [12; 


/* Keep running sum of amp hours*/ 
/* out during phase 0. */ 


Cur File #2 
batt avqt 


fault . data 
fault [2] ; 


/* Average temp of six sensors */ 
/'■*■ for each of 6 batteries. ■*■/ 


/ * Fault flag, fault [0] =1 fault */ 
/* =0 no fault*/ 

/ * f a u 1 1 [ 1 ] = “1 when p r o b 1 e m : i- ; / 

/* with communication */ 
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DATA-HDL.C 


#de f i n e C I L I MI T 
#in elude "hst.h" 


i n t c o n s e q _ i nc m p 1 1 = 0 ; 


i i:4 ? 




**:+:*:* 





.********.*:*:* 

-PM 

+ ! 

■ i. !•' 

( a r gc 

, argv) /■+• maj 

n 

is it 

e Data- 

•Handle? 

driver. It 

C on t r o 

is 

i: 

urge - 


/* pre 

>g r 

am f 1 

ow ar 

id 

deter mi 

ne 

s the times 

for e 

0 

n t 

**a r 

gv ; 












i n t 

orbi 

tno = o.err; 











i n t 

phas 

el , phase2 ; 






/* 

phasel ~ .1 

ast ph 

as 

e 


phas 







/* 

phase2 = c 

u r r eri t. 

pha 

i n t 

e_signal ; 






/* 

EOD phase 

from 0 

17 

o 









/* 

EOC phase 

f r o/n 1 

f. 


.1. n i t 

O ; 







/ :+ 

initialize 

tiori c 

a I 

] 

wli i 1 

•i 

e ( 1 ) 












]_ 

con 

seq_incmpl t 

z 

0; 










wh i. 

le ((err = 

ea 

d_dat 

cl ( / ) 

= = 

FAIL) 

/* 

Read telemetry b 

,;r 

S t 



incomplete 

<> 

•> 




/:+ 

Check for 

c on sec 

.it 

.1 V 


if 







/•+ 

incmplt- da 

ta lour 

•V t 

s 


( (err = chec 

k_ 

fault 

( ) ) “ 


FAIL) 

/:* 

Check for 

fau i. ts 




[ 







/ * 

If fault s 

et fau 

i 

X l. 




fault [0] = 

1 

; 




/'i 

flag - 1 a 

n d f i n 

is 

h 



wf (FAULT'D A 

T) 

’> 










}■ 

finish ( ) ; 



















/ + 

deter mine 

phase 

::h 

an 


i f ( 

( phase 1 = hi 

d [ 

0] . ph 

ase) 

” - 

(phase 

9 

-- hid[l].ph 

ase ) ) 





phas 

e_ 

sign a 

1 = EVM 

IN; 







i f 

(iphasel && 

ph 

ase2) 

phas 

e_ 

s i gri a 1 

= 

EOD ; 





i f 

(phase! && i 

Ph 

ase2) 

phas 

e__ 

signal 

r. 

EOC ; 












/* 

EuD and EOC case; 

“V 



i f 

■{ 

( (phase_s.ign 

a 1 

== EOD) 1 

J 

(phase_ 

si. 

gna.1 == EOC 

) ) 





process (ph 

as 

©-Signal , 

or 

bit.no) ; 

/* 

process da 

ta EOC, 

'E 

3D 



if (phase_s 

ignal = 

= EOC 

) 









print 

f ( 

“ «V ..J 

/to Q . 

ORBIT 

= %d\n " 

, + 

+or bit.no, h i 

d [ 0 ] .o( 

.. 

i. t 


} 




















/ * 

EVMIN case 





p r a 

'-•ess (EVMIN, o 

rb 

i tn o ) 




/* 

process da 

ta EVM 

T i| 


}- / :fc: 

end 

while loop 

/ 
’ / 










f i. n .1 

sh ( ) 

f 






/* 

exiting routine 




H- / 

■■■r / 
■Y- / 
A / 


} /*■ end main */ 


.*/ 

*/ 

* / 
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4:4: i::+. 

+ 4 

4:4 

4: 

4: 4:4::+: 4 

:+::+: 

4: 4:4: 4::+:. +:*::+ 

4 : 4 : 4 : 

4: :+::+::+: 


4:4: 4: 4 

4: 4: 4: :+::+: 4:4::+ 

4:4:4 

4:4:4 +:4:4 

He 4 

4 : 4 : 4 ::+: / 

) 




/ + 

0 

ts the C 

ommunic 

at ion por 

t arid signal 

, i n i t i a .1. 

i z 

0 s x / 





/+: 

i n 

terrupt 

h an 

dler 

, global 

buff 

ers and 

tel 

erne try 

r 

ead4/ 

i n t 

er 

r , 

P 

:> r t ; 












i n t 

X 

Q 

r; 

atch < 

) ; 











s i gn 

al 

fc 

*r 

X 

HINT , 

si 

3 _ catch ) 




/ 

in i t ial i 

:: 0 

user i. 

n 0 

|J. f -t- / 











/ 4 

signal c 

a tc 

h e r 


.4. / 

ser i 

rii 

( ) 

; 







/:+: 

initiali 

ze 

i n t e r r 

up 

t 4 / 











/:+: 

: 

handler 




^ / 

por t 



0 

t_por 

t ( 

( i n t ) 0 ) 




/4: 

set comm 

port for 

R 

-> d j, .:p. 

df_i 

i"! i. 

t ( 

) 







/ 4: 

1 n i. t i. a .1. i 

Z 0 

data b 

u. f 

t v? y f 

r e a d 

X 

n i 

t 

( j .. 






/:+' 

ini t i a 1 i 

z e 

tel erne 

tr 

y ; i- / 











/:+: 

reading 

to 

orbit 

St 

ar t-t/ 


} /* end in it •+•/ 


/ * * * * * * * * * * 4 : * * 4 - ■¥ ¥ •+■ 4 - - 4 : :+: * :+: 4 = * 4 - 4 - 4 - ■¥ 4 - 4 : +' :+: * * 4 : 4 : :+: * 4 : 4 : :+: 4 - + 4 : :+: :*: * :+: * * * :+: :*::+: :+: :+: :+: :+: :+: :+: :+: / 

read_init ( ) /•+ Initializes telemetry reading to first full orbit. +/ 

-i 

i n t i , e r r ; 



i = DEFL.T ; 
whiled ! - 1 ) 

{ 

while ((err = read_data() ) == FAIL) 
incomplete ( ) ; 
i. -• h i d [ 1 ] . n i. g h t _ m .1 n ; 

}■ 

pr in tf ("Starting first full orbit\n") ; 
proc e s s ( E V M I N , 0 ) ; 

/* end read_in.it */ 


/'*■ Orbit starts at night ■*■/ 
/■*' minute 1, discharge */ 
/ + phase.. No action til •+••/ 
/ : + : s t a r t f o und . C h e c k s f o r / 
/ * c: o ri s e c u t i. v e i n c o mp 1 e t e +■ / 
/ : + : t e 1 e m e t r y b Li r s t s . 4; / 

/* Process first min data*/ 


finish () /:+: Exits Data Handler and writes files. •+/ 

l 


ser r st. ( ) ; 

/ * e ri d i n t e r r r u p t h a n d 1 e r + 

wr i t.e_f .i le ( EOC ) ; 

/4 ; write output files * 

wr i te_f i 1 e ( EOD ) ; 


exit (0) ; 

/* exit Data-Haridler 

end finish */ 



A- 2-fl 
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*«*:***« by operator to V 

: i 9 _c a t ch ( ) ' ! * h a j t n a t a - H a ndl e r - 

pri ntf (-interrupt caught: ex i ting ! \n " ) : 

finish ( ) ; 

****** 7 ^ 7 **'^^ :• 

proress<phase_sigr.al, orbitno) . ss the telemetry . 

int phase.signal, orbxtno; ' ' ' 


_ -i-it ( ohrr'P signal , orb itno ) ; 
p f o c e s s _ u a t. a v P •• a -> c _ ->• 

w r i t e f i i e ( Ph ase.s i. gn a 1 ) ; 

df iriitl ( phase_sigrtal ) ; 

} /* end process ■*■/ 


/ * 
7 * 

i 

/ 


process data 

write files 
r ein t i a 1 i z e bu f f et s 


} /* end process */ M 

/*****************^^ ' incomplete data bursts * 

inc o m p 1 e tev. ) • 


■ »H\n“ con sea inomplt, , ** ! 

pr in tf ( "consecutive - -o\n ,cun_e /;+: inc 


coriseq_incmplt++; __ p ILl M TT i 

if (conseq_mcmplt -- '-• LL - i 


inc consecutive in 
complete telemety i ui»s 
if CILIMIT runs, halt 

. , consecutive incomplete\n ,oILIMxT) , 

printf ("Recieve* ~ - ■ sh(Jtting down ! ! \n" > j 

printf ( telemetry uut_c_. / + Set faui 

fault [0] = 1; 

fault [1] = “1" 
wf (FAULTDAi ) ; 
fin ish ( ) ; 


f 1 ags 


> /* end incomplete */ 



READ-DAT. C 


#def ine T1LIMIT 2 
^define T2LIMIT 180 

#include "hst . h" 

char buf [2048] ; 
int count ; 
float cel 1 [ 23 ] ; 

/ t*******t*tttt ************* tt***tt*t^t**** **************** t***^*******t / 
read_data ( ) /* Read telemetry from DEC LSI-11 over RS232 every 1 min.* 

{ 

int err; 

/* Call proc_sync to read start character ’A’, call proc_head to read*/ 
/* header info, proc_bat() to read battery data, and proc_solar to */ 
/* to read SPA and bus data */ 

proc_sync( ) ; 
count. = 0; 

if ((err = proc_head()) == FAIL) return ( FAI L ) ; 
if ((err = proc_bat()) == FAIL) return( FAIL) ; 
if (.(err = proc__sol ar ( ) ) == FAIL) return ( FAIL) ; 
return (SUCCESS) ; 

} /* end read_data */ 

/***************************************************; ******************** / 
proc_sync( ) /* Synchronize reading data to start of telemetry burst.* 

{ 

char cc; 

while (cc = g'et2()) /* Get character from */ 

i /* interrupt hdlr buffer */ 

if (cc == ’A’) /* If ’ A ’ , t e 1 eme t ry start*/ 

{ 

/* printf("Sync Received ! \n" ) ; */ 

break ; 

} 

if(cc == ’B’) /* If ’B’, shutdown */ 

s 

printf ( "\n\nReceived shutdown signal from DEC LSI-ll!\n"); 
/* fault [0] = 1; 

fault [ 1 j = ? ; 
wf (FAULTDAT) ; */ 

f inish ( ) ; 

} 

} /* end while loop */ 


cc = get 1 ( ) ; 

/* 

Read newline and CR 

*/ 

cc = getl ( ) ; 

/* 

to position pointer 

at */ 

/* end proc_sync */ 

/* 

next character 

*/ 


ft- l -R 




/******^***********************JM****************************************/ 
proc_head() /* Read header data from input buffer. */ 

int err; 


if ((err = read_buffer(9)) == FAIL) 
return (FAIL) ; 


/* Read 9 header items #/ 


sscanf (buf , n %d%d%d%d%d%d%d%d%d" , /* Put header data in buf*/ 

&hid [ 1 ]. year , &hid [ 1 ]. day, /* into structured array */ 

&hid [ 1 ]. hour , &hid [ 1] . min, /* hid[l] in int format */ 

&hid [ 1 j .sec,&hid[l] .orbit, &hid[lj .phase, 

&hid [ 1 ] . day_min , &hid [ 1 ] .night_min) ; 

/* Print header data */ 

/* printf("time = %d\nday = &d\nhour = *d\nmin = S»d\nsec = &d\norbit = S»d\nphase = %d\n, 

day_min = %d\nnight_min = 5sd\n\n” , 

hid[l] .year,hid[l] .day,hid[l] .hour,hid[l] . min , hid [ 1 ] .sec, 

hid [ 1 ] . orbit, hid[l] . phase, hid [ 1 ] .day_min,hid[l] .night_min) : */ 


/* 




ret urn ( SUCCESS ) ; 
end proc_head */ 


A - 2 -t 



/***********************************************************************/ 
proc_bat() /* Read battery data from input buffer. */ 

{ 

int j , k , err ; 

f or < k = 0 ; k<6 ; k++ ) /* For 6 batteries * / 

r 

i 

if ((err = read_buf f er ( 1 ) ) == FAIL) /* Read battery no * / 

return (FAIL) ; 

s scanf (buf , "?«d " , &hi d [ 1 ] . b at d [ k ] . bat t no ) ; /* Put buf */ 

/* contents into h i d [ 1 j */ 

/* printf(”batt no = %d, count = Sd\n" , hid [ I ] . bat d [ k ] . bat tno , coun t ) ; * 

if ((err = read_buf f er ( 23 ) ) == FAIL)/* Get 23 cell voltages */ 
return ( FAIL ) ; 

sscanf (buf , , 

&cel 1 [ 0 ] , &cel 1 [ 1 ] , &ce 1 1 [ 2 J , &cell[3J , &cell[4] ,&cell[5] , 

&cell [6] , &cel 1 f 7 ] , &cell[8] ,&cell[9] , &cell[10] , &cell[il] , 

&cel 1 [ 12 J ,&cell[13],&cell[14] ,&cell[15] , &cell[16] , fccelin? ] , 
&cell [18] , &cel 1 [ 19] , &cel 1 [20 ] , ftcell [21 ] ,&cell [22] ) ; 

for ( J = 0 ; J < 23 ; j++) 

hid[l].batd[k].cellv[j] = c e 1 1 [ j ] ; 

if ((err = read_buf f er ( 23 ) ) == FAIL) /* Get 23 cell pressures */ 
return ( FAIL ) ; 

s s c anf ( buf , " S f % f % f % f % f %i % f% f% f% f % f f % f % f % f % f % f%f% f% f% f% f% f " , 

&ce 1 1 [ 0 ] , &cel 1 [ 1 ] , &ce 1 1 [ 2 ] ,&cell[3] , &cell[4] ,&cell r 5] , 

&cel 1 [ 6] , &cel 1 [ 7 ] ,&celi[8] , &cell[9] ,&cell[10] ,&celii 11] , 

&cel 1 [ 12 ] , &cell[ 13] ,&cell [14] ,fccell[15] ,tcell[16] ,&celi[17 1 , 
&cell [18] ,&cell[19] ,&cell[20] ,&cell [21] ,&cell[22] ) ; 

for ( j=0; j<23; j++) 

hid [ 1 ] . batd [k ] . cel lp [ j ] = cellfj] ; 

if ((err = read_buf f er ( 10 ) ) == FAIL)/* Get remaining 10 bat */ 
return ( FAIL ) ; 

sscanf (buf , , /* Put buf contents */ 

&hid[l].batd[k].batv, /* into hid[l] */ 

&hid [ 1 ] . batd [ k ] . bate, 8th i d [ 1 ] .batd[k] .bprcc, 

&hid [ 1 ] . batd [ k J .batemp[0] , &hid [ 1 ] . batd [ k ] . batemp [ 1 ] , 

&hid [ 1 ] . batd [k ] . bat emp [2] , &hid [ 1 ] . batd f k ] .batemp[3] , 

&hid[l] . batd [ k ] . bat emp [4 ] , &hid [ 1 ] . batd [ k ] . bat emp [5 j , 

&hid [ 1 ] .batd[k].batrecond) ; 

} /* end for loop, battery 1 to 6 */ 

return (SUCCESS) ; 

} /* end proc_bat */ 


(V- 3 



/ *************************************************************** ****** 
proc_solar( ) /* Read SPA and bus data from input buffer. * 

int j,err; 

if ((err = read_buffert. 13)} == FAIL) / * Get 13 SPA values * 

return(FAIL); /* values from buffer * 

/'* Put buf contents into * 
/* hid[ 1 ] * 

Sthidfl] .spac[0] , &hid [ 1 ] . spac [ 1 ] , &hid [ 1 ] .spac[2] , 

&h id [ 1 j ,spac[3] , &hi d [ 1 ] . spac [ 4 ] , & hi d [ 1 j .spac [5] , 

&hid f 1 ] .spac[6] , &hi d [ 1 ] . spac [ 7 ] , &hi d [ 1 ] .spac [8] , 

&hid [ 1 ] . spac [ 9 ] , &hid [ 1 ] . spac [ 10 j ,&hid[l] . spac [ 1 1 ] , 

&hid[l] .spac [12] ); 

if ((err = read_buf f er ( 6 ) ) == FAIL) /* 3 busses, volt and * 

return ( FAIL ) ; /* current for each * 

sscanf (buf , ”%fSsfSf!Sf£f£f " , /* Put buf contents * 

&hid[l] .bd[0] .busv, /* into hid[l] * 

&hid [ 1 ] .bd[0] .busc,&hid[l] .bd[l] .busv,&hid[l] .bd[l] .busc, 

&hid [ 1 ] . b d [ 2 ] . busv , &hid [ 1 ] . b d [ 2 ] .busc) ; 

/* printf ( "solar count = ?«d\n ", count ) ; */ 

/* for(j=0; j<3; j++ ) */ /* Print bus data * 

/ * { 

pr int f ( ”bd [ %d ] . bus v = Ssf\n" , j , hid [ 1 ]. bd [ j ]. busv ) ; 
printf ( "bd [%d] . busc = £f\n ", j , hid [ 1 J. bd [ j ]. busc ) ; 

} */ 

return ( SUCCESS ) ; 

} / * end proc_solar * / 


ft- 4 -i. 



/I***********************-************************************************/ 

read_buf f er ( k ) /* Retrieves telemetry values from interrupt*/ 

int k; /* handler buffer and puts them in buf[j */ 


int i=0, nl_count = 0; 


whi le ( nl_count ! = k) 

/* 

For k data points 

* / 

< 

if ( (buf [ i ] = get 1 ( ) ) == EOF) 

/* 

Get each character 

*/ 

return (FAIL) ; 

/* 

from input buffer and 

*/ 

if (buf [ i ] == ’ \r ’ ) 

/* 

put in buf, including 

* / 

{ 

/* 

newline and CR 

*/ 

count++ ; 

/* 

count values read per 

*/ 


/* 

telemetry run 

*/ 

nl count++; 

/* 

Returns FAIL if EOF 

*/ 

if ( (buf [ i ] = get 1 () ) == EOF) 

/* 

read before end of 

* / 

return (FAIL) ; 

•% 

} 

i++ ; 

} 

buf[i] = * \ 0 * ; 

/* 

telemetry burst 

*/ 

/* 

sets end of data in 

*/ 

return ( SUCCESS ) ; 

/* 

in buf [ i ] 

* / 


} /* end read_buffer */ 

/***********i| c***********************************************************' 
getl() /* get character from interrupt handler buffer * 

char cc; 

long tO , tn , t imer ; 

time(&t0); - /* initial time */ 

while ((cc = b_char()) == EOF) /* get char from interrupt*; 

{ /* handler buffer, EOFs */ 

/* start timer = tn - tO */ 
timer = time(&tn) - tO; /* tn is current time */ 

/* printf ( "get : count = Sid, timer = %ld\n" , count , timer ) ; */ 

if ((count < 370) && (timer > T1LIMIT)) /* check for incmplt*/ 
return(EOF); /* run, T1LIMIT exceeded */ 

} 

return (cc) ; 

} /* end get */ 
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/***********************************************************************/ 
get2() /* get character from interrupt handler buffer */ 


char 

cc ; 





long 

tO , tn , t imer ; 





t ime ( &t0 ) ; 


/ * 

initial t ime 

* / 

wh i 1 e 

( ( cc = b_char ( ) ) == 

EOF) 

/* 

read char from interrupt* 




/* 

handler buffer, EOFs 

*/ 




/* 

start timer = tn - tO 

* / 


timer = time(&tn) - 

tO ; 

/* 

tn is current time 

* / 


if (timer > T2LIMIT ) 


/* 

Exit if T2LIMIT is 

* / 


r 


/* 

exceeded => No Commun. 

. */ 


printf("No communication for 3 minutes; exit ing ! \n " ) ; 
fault [1J = -1; /* write fault flags */ 

f aul t [ 0 ] = 1 ; 
wf (FAULTDAT) ; 
f inish ( ) ; 
i 

J 

\ 

j 

return ( cc ) ; 

} /* end get */ 
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PROCESS. C 


^define DCHGLIMIT 30 
^define CHGLIMIT 50 

*inc 1 ude ” hs t . h" 

int no_druns , no_cruns ; 

/************************************************************************/ 
df_init() /* Initialization of buffers*/ 

r 

l 

int i , j ; 

/* printf(”in df_init\n"); */ 


f aul t [ 0 ] 

= 0; 

/* 

FAULTDAT buffer 

*/ 

f aul t [ 1 ] 

= 0; 

/* 

FAULTDAT buffer 

* / 

no_cruns 

= no_druns = 0; 

/* 

run counters 

*/ 

for ( i =0 ; 

{ 

batt 

i < 6 ; i++ ) 

/* 

for 6 batteries 

* • 

_avgt [ i ] = DEFIT; 

/* 

CURF2D AT buffer 

* / 

avg_ 

temp bufferfi] = 0.0; 

/* 

SH0WF9D AT work buffer 

*/ 

high 

buf f er [ i ] = DEFLT ; 

/* 

SHOWF2D AT work buffer 

*/ 

ahoo [ i ] = 0.0; 

/* 

SHOWF3D AT work buffer 

*/ 

ahi [ 

i] = 0.0; 

/* 

SH0WF3DAT work buffer 

*/ 

rc_orb i t [ i ] = 0 ; 

/* 

SHOWF12DAT buffer 

* . 

trickle[i] = 0; 

/* 

S IIOWF 1 1 D AT work buffer 

* / 

for 

r 

( j = 0 ; j < 12 ; j + + ) 

/* 

for 12 orbits 

* / 

X 

eod_voltage [ i ] [ j 1 = DEFLT; 

/* 

SHOWF1DAT buffer 

*/ 


hc_voltage [ i ] [ j ] = DEFLT; 

/* 

SHOWF2D AT buffer 

* i 

• 

rc_rat io [ i ] [ j j = DEFLT; 

/* 

SHOWF3DAT buffer 

* / 


cv_eod_hv [ i ] [ j ] = DEFLT; 

/* 

SH0WF4D AT buffer 

* / 


cv_eod_hv [ i ] [ j ] = DEFLT; 

/* 

SH0WF4DAT buffer 

*/ 


cv_eod_av [ i ] [ j ] = DEFLT; 

/* 

SH0WF4DAT buffer 

*,•' 


cv he hv[i][j] = DEFLT; 

/* 

SH0WF5DAT buffer 

*/ 


cv_hc_ 1 v [ i J [ j] = DEFLT; 

/* 

SH0WF5DAT buffer 

* / 


cv_hc'_av [ i ] [ j ] = DEFLT; 

/* 

SHOWF5DAT buffer 

*/ 


avg temp[i][j] = DEFLT; 

/* 

SHOWF9D AT buffer 

*/ 


ahof i ] [ j] = DEFLT; 

/* 

SHOWF13DAT buffer 

* 

} 

for 

( 

V 

t ime_tc [ i ] [ j j = DEFLT; 

/* 

SHOWF 1 ID AT buffer 

X ’ 

(j=0;j<23;j++) 

/* 

for 23 cells per batt 

*/ 

cellv eod[i][j] = DEFLT; 

/* 

SHOWF6D AT buffer 

* / 


ce 1 1 v_hc [ i ] [ j] = DEFLT; 

/* 

SH0WF7D AT buffer 

*/ 


cp_eod[i][j] = DEFLT; 

/* 

SHOWF10DAT buffer 

*/ 


cp_eoc[i][jJ = DEFLT; 

/* 

SHOWF10D AT buffer 

*/ 


} 
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for ( j=0; j<48; j++) 

r 


avgt [ i ] [j] = DEFIT; 

/* 

SHOWF8D AT buffer 

*/ 

bc_drc [ i ] [ j ] = 0.0; 

j 

/* 

SHOWF12DAT buffer 

*/ 

wf ( FAULTDAT) ; 

/* 

write no fault yet 

* / 


/*end df_init */ 
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/***********************************************************************/ 
df _ini t 1 ( phase_s ignal ) /* After each telemetry burst is read, or after*/ 

/* EOC or EOD, buffers used for processing need*/ 
int phase_signal ; /* to be re- in i t i al iz ed . */ 

{ 

int i , j ; 


/* printf(”in df_initl \ n ” ) ; */ 

switch (phase_s ignal) 

{ 

case EOC : 

no_cruns - no_druns = 0; 
for ( i = 0 ; i < 6 ; i++ ) 


/* 

start switch 

*/ 

/* 

EOC 

case 

*/ 

/* 

run 

counters 



avg_t emp_buf f er [ i ] = 0.0; 
high_buf f er [ i ] = DEFLT; 
ahi [ i ] = 0.0; 
ahoo [ i ] = 0.0; 


/* 

SHOWF9DAT 

work 

buffer 

*/ 

/* 

SHOWF2D AT 

work 

buffer 

* / 

/* 

SHOWF3D AT 

work 

buffer 

*/ 

/* 

SHOWF3D AT 

work 

buffer 

*/ 


for ( j=0 ; j <23 ; j++) 

{ 

cel lv_hc [ i ] [ j ] = DEFLT; 
cp_eoc[i][j] = DEFLT; 
cp_eod[i][j] = DEFLT; 
cel lv eod [ i ] [ j] = DEFLT; 

} 

for ( j = 0; j <48 ; j++) 

avgt [ i ] [ j] = DEFLT; 

J 

break ; 


/* for 23 cells per batt * / 


/* SHOWF7D AT buffer */ 

/* SHOWF10DAT buffer */ 

/* SHOWF1 0DAT buffer */ 

/* SH0WF6D AT buffer */' 

/* SHOWF8D AT buffer */ 

/* end EOC case */ 

/* No EOD case */ 


case EVMIN: /* Every Minute case */ 

for ( i =0 ; i < 6 ; i++ ) /* for 6 batteries */ 

bat t _avgt [ i ] = DEFLT; /* CURF2D AT buffer */ 

hidtohid(T; /* switch hid columns */ 

break; /* end EVMIN case */ 

} /* end switch */ 

} /* end df_initl */ 
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/*********************************:M************************************/ 
hidtohid() /'* After each telemetry burst is read and processed, the */ 

/* data in hid[l] is put in hid[0] to prepare for next */ 

/* data burst which will be read into hidfll. */ 

{ 

int i , j ; 


/* printf("in hidtodhid, ready to exchange buffers\n"); */ 

hid[0].year = hid[lj.year; /* switch header items 

hid[0].day = hid[l].day; 

hid[0].hour = hid[l].hour; 

hid [0]. min = hid[l].min; 

hid[0].sec = hid[l].sec; 

hid[0]. orbit = hi d [ 1 ] . orb i t ; 

hid[0]. phase = hid [ 1 ]. phase ; 

hid[l]. phase = DEFLT; 

hid [ 0 ] . day_min = hid [ 1 ] . day_min ; 

hid [ 1 ] . day_min = DEFLT; 

hid [ 0 ] . night_min = h i d [ 1 ] . n ight_min ; 


for ( i =0 ; i < 6 ; i++ ) /* switch battery items * / 

{ 

hid [0 ] . bat d [ i ] . battno = hid [ 1 ] . batd [ i ] . battno ; 

for ( j = 0; j<23; j++) /* for 23 cells per batt */ 

{ 

hid [ 0 ] . bat d [ i ] . cel lv [ j ] = hid [ 1 ]. batd [ i ]. cel lv [ j ] ; 

hid[l] . batd [ i ] . cel lv [ j ] = DEFLT; 

hid [ 0 j. batd [ i ]. cel lp [ j ] = hi d [ 1 ] . bat d [ i ] . ce 1 lp [ j J ; 

hid [ 0 ]. batd [ i ]. batv = hid [ 1 ]. batd [ i ]. bat v ; 
hid [ 0 ]. batd [ i ]. bat c = hid [ 1 ] . bat d [ i ] . b at c ; 
hi d [ 0 ] . bat d [ i ] . bprcc = hid [ I ] . batd [ i ] . bprcc ; 

> 

for ( i=0 ; i < 1 3 ; i++ ) /* switch SPA items * */ 

hid [ 0 ] . spac [ i ] = hid [ 1 ] . spac [ i ] ; 
hid [ 1 ] . spac [ i J = DEFLT; 


for (i=0;i<3;i++) 

{ 

hid[0] .bd[i] .busv 
hid[0] .bd[i] . busc 
hidfl] . bd [ i ] .busc 

J 

} /* end hidtohid */ 


/* switch bus items 

hid[l] . b d [ i ] .busv; 
hid[l] .bd[i] . busc ; 

DEFLT; 


*/ 
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/***********************************************************************/ 
move_buffers(phase_signal) / * Prepares the global buffers for move */ 
int phase_s i gnal ; 


/* 

printf(” in move_buffers, 

ready to move buf fers \n ") ; */ 


switch (phase_si gnal) 

/* start switch 

* / 

case 

EOC : 

/* start EOC case 

* 


move ( hc_vo 1 t age ) ; 

/* SH0WF2D AT buffer 

* / 


move ( avg_t emp ) ; 

/* SH0WF9D AT buffer 

* / 


move ( t ime_ t c ) ; 

/* SH0WF1 ID AT buffer 

* •' 


move ( cv_hc_hv ) ; 
movef cv_hc_lv) ; 
move ( cv_hc_av ) ; 

/* SH0WF5DAT buffer 

* •' 


break; 

/* end EOC case 

*/ 

case 

EOD : 

/* start EOD case 

*/ 


move(eod_ voltage) ; 

/* SH0WF1DAT buffer 

* / 


move(rc_ratio) ; 

/* SH0WF3D AT buffer 

* / 


move ( aho ) ; 

/* SH0WF13DAT buffer 

*/ 


move ( cv_eod_hv) ; 
move(cv_eod_lv) ; 
move ( cv_eod_av ) ; 

/* SH0WF4D AT buffer 

* / 

} /* 
} /* 

break ; 

end switch */ 

end move_buffers ( ) */ 

/* end EOD case 
/* no E VMIN case 

*/ 

* . 


/ft**********************************************************************/ 

move(buffer) /* If the Data-Handler contines after 12 completed orbits*/ 
/* the buf fers ( 6 , 12 } will lose their first column, the */ 
/* remaining data will be shifted one column to the left,*, 
/* the next orbits data will be put in the 12th column. */ 

float buffer [6] [ 12 ] ; 

{ 

int i , j , j j 5 


f or ( i = 0 ; 

r 

i <6; i++ ) 

/* 

for 6 battery rows 

*/ 

i 

for 

( j = 0 ; j< 10 ; j++) 

/* 

for 11 orbit columns 

* / 

r 


/* 

column 1 is dropped 

* / 


jj = j+1; 

/* 

shift each column one 

*/ 

■* 

buffer [ i ] [ j ] = buf fer [ i ] [ j j J ; 

/* 

column to the left 

* 

buf f er [ i ] [11] = DEFLT; 

/* 

initialize 12tn column*./' 


} 

} /* end move */ 
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/***********************************************************************/ 
process_data ( phase_s i gnal , orbitno) /* Perform mathematical operations */ 

/ * to summarize telemetry. Keeps * 

int phase_sig'nal , orbitno; /* results in global buffers. */ 

{ 

int current_min , col ; 

int j j , i , j ; 
float x , xl , x2 , x3 ; 
float sum , avg ; 

/* printf(”in process_data\n" ) ; */ 
if (orbitno > 11) 

r 

i 

col = 11; 

move_buffers ( phase_s i gnal ) ; 

}else col = orbitno; 


switch ( phase_s i gnal ) /* start switch */ 

r 

l 

case EOC : /* start EOC case */ 

jj = no_cruns + no_druns; 

if (jj < ( CHGLIMI T~+ DCHGLIMIT) ) /* Enough runs? 

break ; 

for ( i =0 ; i < 6 ; i++ ) /* for 6 batteries */ 

{ 

avg_temp [ i ] [ col ] = avg_temp_buf f er [ i ] / j j ; /* SH0WF9DAT buffer*/ 

if (ahoofi] != 0.0) /* SH0WF3D AT buffer */ 

rc_rat io [ i ] [ col ] = ahi f i ] /ahoo [ i ] ; 

if ( trickle [ i ] != 0) /* SH0WF11DAT buffer */ 

t ime_tc [ i ] [ col ] = trickle[ij; 

f or ( j =0 ; j < 23 ; j++ ) /* for 23 cells per batt #/ 

cp_eoc[i][j] = hid [ 0 ] . batd [ i ] . cellp [ j ] ; 

/* SHOWF10DAT buffer */ 

} 

break; /* end EOC case */ 
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case EOD : 

if (no.druns < DCHGLIMIT ) 
break ; 

f or ( i = 0 : i < 6 ; i++ ) 


/* start EOD case */ 

/* Enough discharge runs?*/ 


/ * for 6 batteries 


eod_voltage [ i ] [ col j = hi d [ 0 ] . bat d [ i 1 . bat v ; 

/* SH0WF1D AT buffer 

aho [ i ] [col] = ahoo [ i ] ; /* SHOWF 13DAT buffer 


xl = DEFLT ; 
x2 = -DEFLT; 
x3 = 0.0; 
for( j=0; j<23; j++) 


/* high, low, avg of ceil*'' 
/* voltage at EOD * 

/* for 23 cells per batt * ,/ 


x = hid [ 0 ] . bat d [ i ] . cel Iv [ j ] ; 
if (x > xl) xl = x; 
if (x < x2 ) x2 = x; 
x3 = x3 + x; 

cel lv_eod [ i ] [ j ] = hid f 0 J . batd [ i ] . cellv [ j ] ; 

/* SHOWF6D AT buffer 

cp_eod [ i ] [ j ] = hid[0] .batd[ i] . cellp[ j] ; 

} “ /* SHOWFIODAT buffer 

cv_eod_hv[i] [col] = xl; /* SHOWF4DAT buffer 

cv_eod_lv [ i ] [ col ] = x2; 

cv_eod_av[i] [col] = x3 / 23.0; 


break ; 


/* end EOD case 
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case EVMIN: 

current_min = hidfl] 

jj = -i; 

if ( ! ( cur rent_min % 
for ( i = 0 ; i < 6 ; i++ ) 


/* start EVMIN case 
night_min + hid [ 1 ] . day_min ; 

2)) jj = current_min/2 ; 

/* for 6 batteries 


*/ 


* / 



} 



sum = 0.0; 

for ( j=0 ; j <6 ; j++ ) /* for 6 batt temp sensors*/ 

sum = sum + hid [ 1 ] . bat d [ i ] . bat emp [ j ] ; 


av g = sum / 6.0; 

/* 

avg temp per battery 

*/ 

batt_avgt[i] = avg; 

/* 

CURF2D AT buffer 

*/ 

avg_temp_buf f er [ i ] = avg 

t emp buf f er [ i ] 

+ avg; 



/* 

SHOWF9DAT buffer 

*/ 

if ( j j ! = -1) 

{ . 

avgt[i] [jj] = avg; 

/* 

every 2 minutes 

*/ 

/* 

SHOWF8D AT buffer 

* / 

if ( hid [ 1 j . batd [ i ] . 

batrecond) /* 

SHOWF12DAT buffer 

% / 


{ 

if ( ( rc_orb it [ i ] != hid [ 1 ]. orbit ) && 

(rc_orbit [ i] ! = 01) 

{ 

rc_orbit[i] = hid [ 1 ] . orb it ; 
for ( j=0; j<48; j++) 

bc_drc [ i ] [ j ] = 0.0; 

bc_dr c [ i ] [ j j ] = hi d [ 1 ] . b at d [ i ] . bat c ; 

J 

} 
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( hid [ 1 ] . phase ) 

/* 

CHARGE PHASE 

* / 

no_cruns++ ; 

/* 

increment charge runs 

* / 

for ( i=0 ; i <6 ; i++ ) 

/* 

for 6 batteries 

* , 

xl = high buffer[ij; 

/* 

high, low, avg of cei 

1 * -■ 

x2 = -DEFLT ; 

/* 

volts at hig'h-charge 

* 

x3 = 0.0; 




f or ( j = 0 ; j < 2 3 ; j++ ) 

/* 

: for 23 cells per bat 

t * / 


1 

x = hid[l].batd[ij.cellv[j]; 
if (x > x 1 ) xl = x; 
if (x < x2) x2 = x; 
x3 = x3 + x ; 

} 

cv_hc_lv[i] [col] = x2 ; /* SH0WF5DAT buffer *./ 

if ( high_buf f er [ i ] < xl) 

{ 

cv_hc_hv [ i ] [ col ] = xl; 

cv _hc_av [ i ] [ col ] = x3 / 23.0; 

hc_voltage [ i ] [col ] = hid [ 1 ] . batd [ i ] . batv ; 

/* SHOWF2D AT buffer */ 


} 

] else { 
if 

{ 


} 

} 

break ; 


for ( j=0; j<23; j++) /* SHOWF7D AT buffer */ 

cel lv_hc [ i ] [ j J = hid [ 1 ] . bat d [ i ] . cel 1 v [ j ] ; 

J 

x = hid [ 1 ]. batd [ i ]. bate ; 

ahifi] = ahi[ij * (x/60.0); /* SHOWF3DAT work buffer *' 

printf("ahi[%dj = *f \n" , i , ahi [ i j ) ; 
if ( x < 3.0) 

trickle[i] = trickle [ i ]++; /* SHOWF11DAT work buffer*/ 


/* end CHARGE case * 

( ! hid [ 1 ] .phase) /* DISCHARGE PHASE */ 

no_druns++; /* incr discharge runs */ 

for( i=0 ; i<6; i++) /* for 6 batteries */ 

ahoo[i] = ahoo[iJ - ( hid [ 1 ]. bat d [ i ]. bat c/60 . 0 ) ; 

/* SHOWF3D AT work buffer */ 

/* end DISCHARGE case * 

/* end EVMIN case */ 


/* end switch */ 

/* end process_data( ) */ 
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/t ********************************************************************** ^ 
check fault() /* check telemetry for a fault */ 

int i , j ; 
float x , sum; 

for ( i=0 ; i < 13 ; i++) 

x = hid[l].spac[i]; 

if ( ( hid [ 1 ] . day_min < 5) && (x < 5.0)) 
return ( FAIL) ; 
j = (i + 1 ) / 2 ; 
if ( j ) 

{ 

if (x > = 16.0) r eturn ( FAIL ) ; 

} else{ 

if (x >= 8.0) return ( FAIL) ; 

J 

if (! hid [ 1 ]. phase && (x > 5.0)) return ( FAIL ) ; 

J 

for ( i = 0 ; i < 6 ; i++ ) 

{ 

sum = 0.0.; 

for ( j = 0; j < 6 ; j++) 

{ 

sum = sum + hid f 1 ] . batd f i ] . bat emp [ j ] : 

} 

if ((sum > 25.0) I I (sum < -10)) return ( FAI L ) ; 

for (j = 0;j<23;j->-+) 

/ 

x = hid [ 1 ]. batd [ i ]. cell v [ j ] ; 

if ((x <= 0.0) ! ; (x > 1.55)) return (FAIL) ; 

} 

} 

sum = 0.0; 

for (i=0;i<3;i++); 

{ 

x = hid [ 1 ] . bd [ i ] . busc ; 
sum = sum + x; 

if (! hid [ 11 . phase && (x < 5.0)) return ( FAIL) ; 

> 

if (sum > 99.0) return (FAIL) ; 
return (SUCCESS) ; 

} /* end check_fault */ 
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WRIT-FIL. C 


^include "hst . h” 


■ t***********t*t********tt******ttt*ttt*t****t*ttt*t*^**t**t*t *********** / 
wr i t e_ f i 1 e ( phas e_s i gna 1 ) /* Determines which data files should *- 

int phase_s i gnal ; /* be written at this time. */ 

{ 

int err; 


/* printf("at wr i t e_f i 1 e\n " ) ; */ 

swi t ch ( phas e_s igna 1 ) /* start switch */ 

{ 


case EOC: 

if ( ! (err 
if ( ! ( err 
if ( ! (err 
if ( ! ( err 
if (! (err 
if ( f (err 
if ( ! ( err 
if ( ! ( err 
if ( ! ( err 
break ; 


wf ( SH0WF2DAT ) ) ) printf(” 
wf ( SH0WF3D AT ) ) ) printf(" 
wf ( SH0WF5D AT ) ) ) printf(" 
wf ( SH0WF7D AT ) ) ) printf (" 
wf ( SH0WF8D AT ) ) ) printf(" 
wf ( SH0WF9DAT ) ) ) printf(" 
wf (SH0WF11DAT) ) ) printf( 
wf ( SH0WF12D AT ) ) ) printf( 
wf (SHOWFIODAT) ) ) printf( 


/* 

Couldn ’ 
Couldn ’ 
Couldn * 
Couldn ’ 
Couldn ’ 
Couldn ’ 
"Couldn 
"Couldn 
"Couldn 
/* 


start EOC case 


*/ 


open %s \n " , " showf 2 . da t ’’ ,) ; 
open %s\n” , "showf3 . dat" ) ; 
open S»s\n" , ”showf5 . dat" ) ; 
open %s\ n" , " showf 7 . dat ” ) ; 
open Ss\n" , " showf 8 . dat " ) : 
open -»s \ n " , " s howf 9 . dat " ; : 
’t open %s\n" ," showf 1 i . da t " 

* t open ?*s \n” , " showf 12 . dat " ) 
’ t open *£s\n" , "showf 10 . dat" ) 
end EOC case */ 


case EOD: 

if ( ! ( err 
if ( ! ( err 
if ( ! ( err 
if ( ! (err 
break ; 


wf (SHOWF1DAT) ) ) printf ( "Couldn ’ t open % s \ n " , ” s howf 1 . dat " : 

wf (SH0WF4DAT) ) ) printf ( "Couldn ’ t open %s\n" , "showf4 . dat " ) 
wf ( SHOWF6D AT ) ) ) printf ( "Couldn ’ t open ?»s \n" , " s howf 6 . dat " ) : 
wf ( SHOWF1 3D AT ) ) ) printf ( "Couldn ’ t open ?*s \n " , " s t at f 1 . da t ” ' : 

/* end EOC case */ 


case EVMIN: 

if ( ! ( err 
if ( ! ( err 
if ( ! ( err 
break; 


/* 

wf ( CURF2D AT ) ) ) printf ( "Couldn ’ t 
wf (CURF3DAT) ) ) pr in t f ( " Cou 1 dn ’ t 
wf (CURF1DAT) ) ) printf ("Couldn’ t 

/* 


start EVMIN case */ 

open \n" , " cur f 2 . dat " ) ; 

open S»s\n" , "curf 3 . dat " ) ; 
open fcs\n" fault . dat ") : 
end EVMIN case */ 


} /* end switch */ 

} /* end write_file */ 
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/ **t ******************************************************************** / 
wf(filename) /* Write output files for Expert System from buffers */ 

/* containing summarized and condensed telemetry. */ 

int filename; /* Files are written at EOC, EOD or Every Minute. * 


FILE *s f p : / * Show File# Pointer*./ 

int i , j , err ; 


switch 

f 

v. 

case 


( f i lename ) 

/* 

Write 

data outpu t file 

* 

SH0WF1DAT : 

/* 

Write 

showfl.dat, contains 

* 

err = wr i te_ f 1 (” showf 1 . da t " 

, eod_voltage, 1) ; 

/* EOD voltages per 

* 

return ( err ) ; 



/* orbit, per bat. t 

* 


case SH0WF2DAT: /* Write showf2.dat, contains * 

err = wr i te_f 1 ( " showf 2 . dat ” , hc_ v o It age , 2 ) ; /* hi gh_ vo 1 t ages per * 
return(err); /* orbit, per battery* 


case SH0WF3DAT: /* Write showf3.dat, contains * 

err = wr i t e_f 1 ( " showf 3 . dat " , r c_r at i o , 3 ) ; /* recharge ratio per * 

return(err); /* orbit, per battery * 


case SH0WF4DAT: /* Write showf4.dat, contains * 

err = write_ i f2( "showf 4.dat",cv_eod_hv,cv_eod_lv,cv_eod_av,4); 
return(err); /* high, low & avg cell volts * 

/'* per orbit and per battery * 

case SH0WF5D AT : . /* Write showf5.dat, contains*, 

err = write_f 2 ( ”showf5 . dat ” , cv_hc_hv , cv_hc_lv , cv_hc_av , 5 } ; 
return ( err ) ; /* high, low & avg of 23 cell* 

/* voltages at high charge * 
/* per orbit, per battery * 

case SH0WF6DAT: /* Write showf6.dat, contains* / 

err = write_f3( "showf6. dat" , cellv_eod, 6) ; /* 23 cell voltages at*/ 

return(err); /* EOC, per orbit, per battery*/ 

case SH0WF7DAT: /* Write showf7.dat, contains*/ 

err = wr i te_f 3 ( " showf 7 . dat ” , cel lv_hc , 7 ) ; /* 23 cell voltages at */ 

/* high charge per */ 

return(err); /* orbit, per battery *./ 
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case 


case 


case 


SH0WF8DAT: /* Write showf8.dat, contains* ; 

if ((sfp = f open ( ” showf 8. dat ” , " w" ) ) >= 0) /* average temp per */ 

{ /* batt every 2 mins */ 

fprintffsfp, "showf (8, [ " ) ; 
for (i=0; i<6; i++ J 

r 

fpr intf ( sfp , 

for ( j = 0 ; j < 48 ; j++) 

f 

\ 

if (j != 47) fprintffsfp, " % f , " , a v g t [ i ] [ j j } ; 
else fprintf f sfp, "?®f", a vgt[ij [j] ) ; 

} 

if (i ! = 5) fprintf (sfp, " ] , \n” ) ; 
else fpr int f ( sfp ,"]]).”) ; 

} 

fclose(sfp) ; 
return (SUCCESS) ; 

}else return ( FAIL ) ; 

SH0WF9DAT: /* Write showf9.dat, contains*, 

err = wr i t e_f 1 ( "showf 9 . dat " , avg_ temp , 9 ) ; /* average temperature * 

return(err); /* per orbit, per batt * / 

SHOWFIODAT: /* Write showfl0.dat, contains* 

if ((sfp = f open ( "showf 10 . dat ", "w" ) ) >= 0) /* 23 cell pressures * 

{ /* at EOC & EOD for last full * 

fprintf(sfp, " showf ( 10 ,[") ; /* orbit, per battery * 

for (i=0; i<6; i++ ) 

{ 

fprintf (sfp, 

for ( j = 0 ; j < 23 ; j++) 

fprintf(sfp, "%f , ", cp_eod [ i ] [ j ] ) ; 

for (j=0; j<23; j++) 

{ 

if (j != 22) fprintf (sfp , "^f, ", cp.eoc [ i ][ j ]) ; 
else fprintf(sfp,”S»f”,cp_eoc(ij [j]); 

} 

if (i != 5) fprintf(sfp, "],\n”); 
else fprintf(sfp, ”]])."); 

} 

f close ( sfp) ; 
return (SUC CESS) ; 

}else return ( FAIL) ; ' 
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case 


c'as e 


case 


case 


SHOWF 1 1DAT: /* Write showfll.dat, contains* 

err = write_f 1 ( "showf 1 1 . dat" , t ime_tc , 11 ) ; /* trickle time per * 

return ( err ) ; /* orbit, per battery 

SH0WF12DAT: /* Write showfl2.dat, contains* 

if ( ( s f p = f open ( "showf 12 . dat ", ”w” ) ) > = 0) /* battery current * 

i; /* during reconditioning at * 

fprint f ( sfp , "showf (12, [ " ) ; /* 2 min intervals per battery 

for (i=0; i<6; i++ ) 

r 

fprintf(sfp, ”[S»d,",rc_orbit[i]); 

for (j=0; j<48; j++ ) 

if (j != 47) fprintf (sfp , "S>f , " , bc_drc [ i] [ j ] ) ; 
else fprintf(sfp, ”Sf”,bc_drc[i] [j] ) ; 

} 

if (i != 5) fprintf(sfp, "],\n”); 

else fprintf(sfp, "]])."); 

} 

fclose(sfp) ; 
return(SUCCESS) ; 

}else return ( FAIL) ; 

SHOWF13DAT : /* Write showl3.dat, contains* 

err = write_f 1 ( "showf 13 . dat” , aho , 13 ) ; /* AHO per orbit, per battery * 

return ( err ) ; 

CURF1DAT: /* Write curfl.dat, contains * 

if ((sfp = f open ( " curf 1 . dat ” , "w" ) ) >= 0) /* -orbit number k * 

{ /* reconditioning flags per ■■ 

fprintf (sfp , ”curf(l, [%d, [”,hid[l] .orbit) ; /* battery r 

for ( j = 0 : j < 6 ; j++) 

{ 

if (j 5) ‘fprintf (sfp , "%d, ", hid [ 1 ] . batd [ j ]. batrecond ) ; 

else fprintf(sfp, "Ssd] ]).", hid [ 1 ]. batd [ j ]. bat recond i ; 

} 

fclose(sfp) ; 
return (SUCCESS) ; 

}else return ( FAIL) ; 
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case 


case 


CURF2D AT : /* Write curf2.dat, contains*/' 

if ( ( sf p = fopen("curf2.dat" , ”w" ) ) >= 0) /* AHI.AHO, phase,*:/ 

{ /* 13 SPA currents, average */ 

/* temperature for each batt*/ 

fprintf (sfp, "curf (2 , [%d, %d, \n [ " , hid [ 1 ] . phase, hi d [ 1 ] . day_min) 
for (j = 0; j < 1 3 ; j ++ ) 

if (j ! = 12) fprintf (sfp, "S»f, ", hid f 1 ]. spac [ j ]) ; 
else fprintf(sfp,”S»f",hid[l] . spac [ j ] ) ; 

f p r i n t f ( s f p , ” ] , \ n [ " ) : 

for ( j=0; j<3; j++) 

{ 

if (j ! = 2) fprintf (sfp , "*f , " , hid [ 1 ] . bd [ j ] . busc ) ; 

else fprintf(sfp, "3»f",hid[l].bd[j].busc); 

fprintf (sfp, "],\n["); 
for (j=0; j<6; j++) 

r 

i 

if (j ! = 5) fprintf (sfp, "Sf, ”,batt_avgt [j] ) ; 
else fprintf(sfp, ’’Sif] ] ) . ”,batt_avgt[j] ) ; 

} 

fclose(sfp) ; 
return(SUCCESS) ; 

} e 1 s e return ( FAIL ) ; 

CURF3D AT : /* Write curf3.dat, contains *• 

if ((sfp = f open ( " curf 3 . dat " , ”w” ) ) >= 0) /* day mm, night min «• 

{ /* 23 cell voltages per b a 1. 1* 

fprintf(sfp, "curf(3,[”); 

for (i = 0; i < 6 ; i++) 

{ 

fprintf(sfp, 

for (j=0; j<23; j ++ ) 

if(j != 22) fprintf (sfp, "Sf, ", hid [ 1 J. batd [ i ]. cel lv [ j 
else fprintf ( sfp , "JSf" , hid [ 1 j . batd [ i ]. cel lv [ j ] } ; 

} 

if (i ! = 5) fprintf(sfp, ”],\n”); 
else fprintf(sfp, ”]])."); 

} 

f close ( sfp ) ; 
return(SUCCESS) ; 

} else return ( FAIL ) ; 
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FAULTDAT: /* fault.dat has fault */ 

if ((sfp = fopen( "fault . dat" , "w" ) ) >= 0) /* flagl = 0 or 1 */ 
l / * f 1 ag2 = 0 or - 1 * / 

f print f ( s fp , " f aul t ( f *d, %d j ) . ", fault [0] , faul t [ 1 j ) ; 
fclosefsfp) ; 
return ( SUCCESS ) ; 

}else return (FAIL); 

} / * end of switch * / 

} / * end of wf */ 




case 
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/*********************************************************************** 
wr ite_ fl ( f i lename , buf f er , number ) /* Write output file for Expert System* 
char *filename; / * buffers containing summarized telemetry. * 

float buf f er [ 6 1 [ 12 ] ; 

int number; 

FILE *file; 

int i , j ; 

if ((file = f open ( f i lename , "w” ) ) >= 0) 

fprintf(file, ” showf ( %d , [ ” , number ) ; 
for (i=Q; i<6; i ++ ) 

{ 

fprintf(file, 

for ( j = 0 ; j < 12 ; j + + ) 

{ 

if (j != 11) fprintf(file, ” % f , " , buffer[i][j]); 
else fprintf(file, ” S» f " , b u f f e r [ i ] [ j J i ; 

} 

if (i != 5) f print f ( f i le , ”],\n"); 
else fprintf(file, "]])."); 

\ 

j 

fclose(file) ; 
return (SUCCESS) ; 

}else r eturn ( FAIL ) ; 

} / * end write_fl * / 
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/***********************************************************************• 
write_f2( filename ,bufferl,buffer2,buffer3, number ) 

/* Write output files for Expert System from buffers *. 


char ^filename; /* containing summarized and condensed telemetry, 
float bufferl [6] [12] ,buffer2[6] [ 12 ] , buf f er 3 [ 6 ] [121 ; 
int number; 


r 


FILE *file; 

int i , j ; 


if ((file = f open ( f i 1 ename , " w" ) ) >= 0) 

r 

fprintf(file, "showf (%d , [ " , number ) ; 
for (i=0; i<6; i++ ) 

{ 

fprintf(file, ” [ " ) ; 
for (j = 0; j < 12 ; j + + ) 

fprintf(file, "3>f,",bufferl[i] [j]); 

f or ( j = 0 ; j < 12 ; j«-+) 

fprintf(file,"%f,",buffer2[i][j]); 



\ 

J 



for ( j = 0 ; J < 1 2 ; j++) 

{ 

if (j ! = 11) fprintf ( file, "Sf , " ,buffer3 [ i] [ j ] ) ; 
else f print f ( f i le , "%f " , buf f er3 [ i j [ j ] ) ; 

if (i ! = 5) fprintf(file, ” ] , \ n ” ) ; 
else fprintf(file, "]]).”); 

\ 

fclose(file); 
return (SUC CESS) ; 

}else return ( FAIL ) ; 

/* end write_f2 */ 
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/*************************************************** ************* ******* 
wr i t e_f 3 ( f i lename , buf fer , number ) /* Write output file for Expert System* 
char ^filename; /* from buffers containing summarized telemetry. * 

float buffer [6] [231; 
int number; 

FILE *f i 1 e ; 

int i , j ; 

if ((file = f open ( f i lename , "w" ) ) >= 0) 

<. 

fprintf(file, ” showf ( ?id, [" , number ) ; 
for (i=0; i<6; i++ ) 

{ 

fprintf(file, 

for ( j = 0 ; j < 23 ; j++) 

{ 

if (j .' = 22) fprintf (f ile , "3f , " , buf f er [ i ] [ j ] ) ; 
else fpr int f ( f i 1 e , "%f " , buf f er [ i ] [ j ] ) ; 

% 

1 

if (i ! = 5) fprintf ( f i le \n” ) ; 
else fprintf(file, "]]).” ) ; 

} 

fclose(file) ; 
return ( SUCCESS ) ; 

}else return ( FAIL ) ; 

} / % end write_f3 * / 
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APPENDIX B 


CODE AND DOCUMENTATION FOR EXPERT SYSTEM 



PROLOG. INI 

This is the initialization file for ARITY PROLOG. The directives 
in this file are run every time the Prolog interpreter is called, 
by typing "api" from the NICBES directory. Prolog.ini consults 
all the PROLOG programs needed to run NICBES (start. prg, 
faultd.prg, utility. prg, showpak.prg, grafpak.prg, status. prg', 
advice. prg). Any additional routines added to this package should 
be included in prolog.ini. All the data files necessary for 
the execution of the Expert System are copied to the NICBES 
directory. The current data files ( curf (N) . dat , N = 1 to 3) and 
fault. datare loaded. begin/0 is called to initiate the Expert System 



[ ’ utility. prg* , * showpak . prg’ , ’grafpak.prg’ , ’ start . prg’ , ’advice. prg’ , 
’status. prg, ’faultd.prg’J. 

shell(’copy c : \usr\showf * . dat 

shell(’copy c : \usr\curf* . dat .’), 
shell(’copy c:\usr\fault.dat .’). 

[’fault.dat’ , ’curfl.dat* , ’curf3.dat* , ’curf2.dat’ ] . 

begin . 


2 




Start. prg is the main driver for the 
if a fault flag is set. Otherwise is 
that the user, by making selections, 
the system to view next. 


NICBES. 

prints 
can det 


It calls faultd.prg 
menus to the screen so 
ermine what portion of 


PREDICATES AVAILABLE IN START. PRG 


curf ( 1 


t o 


begin/0 — calls fault to read the fault flag. It also calls 

get the orbit number and reconditioning flags, which are both 
asserted to the data base for use by other routines. begin/1 then 
calls eval_flag/l. begin/0 is called by prolog.ini. 


eval_f lag( Flag) — if Flag is 1, contro 
fault diagnosis. When finished, 
information is desired. The use 
and then continue/1 is called, 
the show files to load the data 
to insure that enough values are 
analysis. If there are too many 
allowed to view the graphics por 
is called next to prepare the ho 
drawn later. wri te_message/l th 
eval_flag/l called by begin/0. 


1 is passed 

the user is 
« 


to faultd.prg for 
asked whether more 

r's response is validated by checkl./l 
If Flag is 0, 
as facts. c< 
present in 
missing dati 
tion of the 
rizontal axis 

en brings up the Main Menu. 
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the 
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cont inue ( yes/no ) — If yes, eval_flag(0) is called. If no, all data 
are deleted from the NICBES directory and PROLOG is halted. 
The user is returned to DOS. 


files 


write_message(N) — brings the Main Menu to the screen and waits for user 
input which is checked for validity by check2/2. N is a flag which 
determines the allowable selections by the user. write_message/l 
operates in a repeat-fail loop so that the user can always come back to 
the Main Menu. There are 4 selections in the menu: 


NICBES MAIN MENU 

1. PLOTS AND GRAPHS 

2. BATTERY STATUS 

3. ADVICE ON RECONDITIONING, 

4. Quit NICBES 
SELECTION (e.g. 1<CR> ): 


WORKLOAD AND CHARGE 


ENTER YOUR 

wri te_message/l is called by eval_f lag( 0) . 


It calls battery/1 


battery (Choice) — Choice (1 - 4) is the item selection from the Main Menu. 
battery(4) calls continue(no) to halt the Expert System. battery/1 
asks the user to select a battery. The user’s response is validated by 
check/2. execute/2 is then called to execute the Main Menu Choice for 
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the chosen battery. battery/1 is in a repeat-fail loop so that the the 
user can make multiple battery selections. 

execute ( Choice , Bat ) — is called by wr ite_message/l to pass control to the 
Choice selected in the Main Menu. There are 3 cases. 

execute(l,Bat)— Plots and Graphs. A Graphics Menu is written on the 
screen with 12 plots per battery. The user may also make another 
Battery Selection or Quit to the Main Menu. User’s choice is checked 
for validity by calling check/2. After the graph has been drawn on the 
screen the user is asked to enter *<CR>’ when ready to continue. 
execute/2 is in a repeat-fail loop so that the Graphics Menu will always 
return until the user opts to return to the Main MEnu. show_view/3 
( showpak . prg) is called to generate the plots. 

GRAPHICS MENU FOR BATTERY N 

1 — Battery Voltage at EOD for last 12 orbits 

2 — Battery Voltage at high in-charge for last 12 orbits 

3 — Recharge ratio for last 12 orbits 

4 — Cell Voltages at EOD; high, low, average for last 12 orbits 

5 — Cell Voltages at high in-charge; high, low, avg, last 12 orbit 

6 — Cell Voltages at EOD for latest orbit 

7 — Cell Voltages at high in-charge for latest orbit 

8 — Average Battery Temperature for latest orbit, each 2 min 

9 — Average Battery Temperature for last 12 orbits 

10— Cell Pressures at EOC and EOD for latest orbit 

11— Time on Trickle Charge for last 12 orbits 

12— Battery Current during reconditioning, 1 orbit, each 2 min 

13— Quit for Another Battery Selection 

14— Quit to Main Menu 

Enter choice (e.g. 5<CR> ): 

execute ( 2 , Bat ) — activates the Status analysis portion. A header is 
written to the screen, STATUS FOR BATTERY N, LATEST ORBIT NN. Then 
battery_status/l is called to analyze Battery Status after which the 
user can opt to make another Battery selection and return to Status or 
Quit to the Main Menu. 

execute ( 3 , Bat ) — An Advice Menu of five items is displayed to the user. 
User’s response is checked for validity by check/2. more/2 is called 
to activate the Advice portion of the Expert System. execute/2 is in a 
repeat-fail loop so that the Advice Menu will always return until the 
user opts to return to the Main Menu. 

BATTERY ADVICE MENU FOR BATTERY N 

1. RECONDITION BATTERY? 

2. CHANGE CHARGING REGIME? 

3. CHANGE WORKLOAD? 

4. QUIT FOR ANOTHER BATTERY SELECTION 

5. QUIT TO MAIN MENU 

ENTER CHOICE (e.g. 1<CR> ): 
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more ( Choice , Bat ) — writes a header to the screen, ADVICE FOR BATTERY N, 
LATEST ORBIT NN. advice/2 (advice. prg) is then called, more/2 is 
by execute ( 3 , _ ) . 


called 
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ORIGINAL PAGE IS 
Of POOR QUALfTY 


start .prg is the main control unit for the Expert System. It calls 
fau.tt_diagnosis if fault flag is set to 1. Main Menu and support in 
menus are drawn on the screen. Control is passed to other portions 


the Expert System depending on the user’s selections. 
See start .doc for further documentation. 


begin : •- 

c a 1 1 ( f a u 1 1. ( [ F 1 a g , _ ] ) ) , 
c a 1 1 ( c u r f ( 1 , [ 0 r b i t , R e c o n d ] ) ) , 
a s s e r t a ( o r b .i. t ( 0 r bit) ) , 
a s s e r t a ( r e c o n dt R e c o n d ) ) , 
e va . 1 . ... f .1 ag ( F 1 ag ) . 

eva 1 f 1 ag ( .1 ) : - 
(. .i. s i n i . 

w r i t e ( $ * * * * * * * * * * * * * * + * * * * * * * * * * * * * * * * * * * * * * * * * ** * * * **:* * * * * * * * $ ) 
ril , write ($NICBES FAULT DIAGNOSIS FOR ORBIT $) , 
call (orbit (Orb) ) , write(Orb) ,nl , 
fan. It... diag.nl , 

wri te ($********************#**#***#*************#*********#****$) 
n 1 , 

repeat , 

wri te($ Do you want more information (yes or no) ? $) , 

reader (Ans) , 
check! (Ans) , 
con tinue ( Ans) . 

continue (yes ) : - 

eval_.f lag (0) . 


continue (no) : - 

shell ( ’del showf* . dat 1 ) , 
she! 1 ( ’ del cur f '+ . dat ’ ) , 
shell (’del fault.dat’) , 
halt . 

eval_f lag (Flag) 

[’showfl.dat’ , ’showf2.dat’ , ’showf3.dat’ , ’showf4.dat’ , ’showf5.dat ’ 

’ showf6 . dat ’ , ’ showf 7 . dat ’ , ’ showfS . dat ’ , ’ showf 9 . dat ’ , ’ showf 10 . dat ’ 

’ showf 11 . dat ’ , ’ showf 12 . dat ’ , ’ showf 13 . dat ’ ] , - 

complete ( N ) , 

c a .1 1 (or b i t ( 0 r b i t ) ) , 

s equ e nee ( 0 rbi t ) , 

wr i te_message ( N ) . 
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wr i te.. .message ( N ) 
c .1. s , n 1 , n 1 , 
repeat , 
wr i te ($ 
w r x t e ( $ 
wr i te ($• 
wr i te ( S 
wr x te ($ . 
wr i te ( $ 
n 1 , 

wr i te (S 


NICKEL CADIUM BATTERY EXPERT SYSTEMS ) , n 1 , 

M I SEES M A I N MENUS ) , n 1 , n 1 , 

1 PLOTS AND GRAPHS $),ni, 

2 BATTERY STATUS $),nl, 

3 ADVICE ON RECONDITIONING, WORKLOAD OR CHA 


4 QLJ I T N I CB ESS ) , n 1 , n 1 , 
writeCS ENTER YOUR SELECTION (e.g. 1<CR> ) : $) , 

r e a d e r ( C h o i c e ) , 
check 2 ( C h o i c e , N ) , 
b a 1 . 1 e r y ( C h o i c e ) , 
fail. 


battery (4) 

c on t i ri ue ( ri o ) 


b a 1 1 e r y ( C h o i ce) : - 
c 1 s , n 1 , n 1 , n 1 , 
repeat , 

ri 1 , wr i te ( $ NICKEL CADIUM BATTERY EXPERT SYSTEM $) , 

n .1. , ri I , 

write (S ENTER BATTERY SELECTION 1 - 6 (e.g. 2<CR> ) : $) , 

r eade r ( B a t ) ,nl ,nl , 
check (Bat , 6 ) , 
e x e c u t e (Choice, B a t ) , 
c 1 s , ri 1 , n 1 , 
i , fail . 

e x e c u t e ( 2 , B a t. ) 

W r i t e ( $ * * * * :+: * * * * * * * * * * * * * * * * * * * * + * * * * * * * * * * * * * * * * * * * * * * * * * * $ ) , . r 1 .1 
writeCS STATUS FOR BATTERY $) , wr i te (Bat ) , 
call (orbit (Orbit ) ) , 

wr i te ( S , LATEST ORBIT $ ) , wr i te ( Orbit ) , n 1 , n 1 , 
batt.ery_st.atus (Bat ) , 

n 1 , wri te ( f*****************^********^************:*:***********-*:*:*:**® ) 

ri 1 , n 1 , 

writeCS 1 QUIT FOR ANOTHER BATTERY SELECTION S),n.i , 

writeCS 2 QUIT TO MAIN MENU $),ril, 

writeCS ENTER CHOICE (e.g. 1<CR> ) : $) , 

r e a d e r ( C h o i c e ) , 
c h e c: k ( Ch o i c e , 2 ) , 

((Choice -- .1, !, fail ); (Choice == 2, els, ni, nl) ) . 
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e x e c u t e ( 1 , B a t ) : - 

c a 1 1 ( o r b i t ( 0 r b i t ) ) , 
c i s , n 1 , n .1 , 
repeat. , 

MENU FOR BATTERY $>:, write (Bat) ,nl,nl. 
Battery Voltage at EOD for last 12 orbits 3) 
Battery Voltage at high in-charge for last 1 
Recharge ratio for last 12 orbitsS) , ril , 

Cell Voltages at EOD; high, low, average for 
Cell Voltages at high charge; high, low avg. 


write($GRAPHIC! 


wr i te ($ 
wr i te ($ 
write (3 
wr i te ($ 
wr i te ( $ 
nl, 

write (3 
wr i te (3 
write (3 
wr i te ( 3 
write (3 
wr i te ( 3 
write (3 
nl, 

wr i. te(3 
write (3 


last 


bits 

12 

12 


u r U ± I 

orb 1 1-: 


9 

10 

11 

12 


14 


Cell Voltages at EOD for latest orbit 3), nl. 

Cell Voltages at high in-charge for latest orbit 3), ni. 
Average Battery Temperature for latest orbit, each 2 mm3. 
Average Battery Temperature for last. 12 orbitsS) ,nl , 

Cell Pressures EOC and EOD for last orbit3),nl. 

Time on Trickle Charge for last. 12 orbits3) , nl. 

Battery Current during reconditioning, 1 orbit-each 2 min 


wr i te ( 3ENTER CHOICE 
read e r ( C h o i c e ) , 
c h e c k ( Cl h o ice, 14) , 

( ( Ch o i c e = = 13, ! , fail) ; 

( ( ch o i c e = = 1 4 , c 1 s , n 1 , n 1 ) ; 
( s h o w _ v i e w ( C h o i c e , 
tmove ( 0 , 0 ) , 
write (3ENTER <CR> 
reader (A, ns ) , 
gc (full ) , 
c 1 s , f a i 1 ) ) ) . 


Quit for Another Battery Selections) , ri 1 , 
Qu i t to Ma in Menu 3 ) , n 1 , 

(e.g. 1 <CR> ): 


3). 


Bat. , Orbit) , 
WHEN READY 


TO CONTINUE: 3) 


BATTERY ADVICE MENU FOR BATTERY 3) , write (Bat ) 


e x e c u t e ( 3 , B a t ) 
repeat , 
write (3 
ri 1 , n 1 , 
write (3 
write (3 
wr i te (3 
wr i te ( 3 
wr i te (3 
write (3 
r e a d e r ( C h o ice ) , 
c h ec k ( Ch o i c e , 5) , 

( (Choice == 4, ! , fail) ; 

( ( Ch o ice = = 5 , c 1 s , n 1 , n 1 ) ; 
(more (Choice , Bat.) , fail ) ) ) . 


ENTER CHOICE 


RECONDITION BATTERY? 3>,nl, 

CHANGE CHARGING REGIME? 3),nl, 
CHANGE WORKLOAD? 3),ril, 

QUIT FOR ANOTHER BATTERY SELECTION 
QUIT TO MAIN MENU S),nl,n.l. , 

(e.g. .1 < CR > ): 3), 


3 ) , n .1 , 
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in o r e ( C h o i. c e,Ba t ) : — 

wr i te ( $jmc*****************************************************$ ) , n 1 , 
c a 1 1 ( o r b i t ( 0 r b i t ) ) , 

write <$ADVICE FOR BATTERY $} , write (Bat ) , 
w r i t e ( $ , L A T E S T ORBIT $ ) , w r i t e ( 0 r b i t ) , n i , n 1 , 

ad v .1. c e ( B a t , Ch o i c e ) , 

uj r i t e ( $ * * * * rr * 4 * * * .+ * * * * * * *** * * * * t. * * * * * * * * # ** * * # * * * i: 4 4 t * + 4. * $ ) , 

n i , n 1 , 

write($ TO CONTINUE ENTER <CR> $) , 
reader (Ans) , els, ! . 
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Faultd.prg is the fault diagnosis section of NTCBES. It uses the 
current data contained in curfl.dat, curf2.dat, and curf3.dat to 
diagnose alarms to the system. Faultd.prg is automatically 
invoked by the system when the C lan gu age data handler detects an 
anomaly. It loads values into the database, diagnoses the 
problem, and cleans up the database, then quits. It is not 
accessible by the user. If more than one fault exists, it will 
diagnose them all. 


ALARMS FOR NICBES SYSTEM: 

More than 18 faults trigger the suggestion to "Check the control 
computer." This is a means to quantify the statement 

"lots of crazy data coming out means the control computer is faulty.” 

The data handler monitors the telemetry as it comes in, and looks for 
values that indicate a serious malfunction of the system, for which the 
DEC LSI-11 has sent an alarm to the engineer. The alarms are in five 
categories : 

1. Power supplies 

Current from any SPA less than 5 A. during first five minutes of 
the day is a failure of power supply or control circuitry. 

Current greater than 8 A. for 1-SPAs (1,3,5,7,9,11) or 16 A. for 2-SPAs 
(2,4,6,8,10,12,13) is a failure in control circuitry. 

a. current < 5 A during first 5 minutes of charge period 

— check the 13 SPA currents (SPA I), phase is 1, day min. <= 5. 

b. current >= 8 A for 1-SPAs, >= 16 A for 2-SPAs 

— check the 13 SPA currents. 

c. current > 5 A during discharge period 

— check the 13 SPA currents, phase is 0. 

2. Batteries 

Cell voltage for any cell less than 0 means a failure in BPRC, 
and voltage of 0 means a short to a cell. Cell voltage for any 
cell greater than 1.55V. is a failure in control limit circuitry. 

a. Cell voltage <= 0 for any cell in any battery 

— check cell voltages 

b. Cell voltage > 1.55 for any cell in any battery 
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— check cell voltages 

3. Load banks 

Tf the sum of the current on busses A, B and C (VI, V2, V3 in 
code) is greater than 99 A., bus overload has occurred, a failure 
of the control circuitry or load banks. If any of the bus currents 
is less than 5 A. during discharge period, a failure of control 
circuitry has occurred. 

a. Sum of three bus currents > 99 A 

— sum the three bus currents (Bus I). 

b. Load < 5 A on any single bus during discharge period 

— check the three bus currents when phase = 0. 

4. Temperature 

If 1 battery has high temperature ( > 25C) during charge cycle, thermal 
runaway is suggested. If 2 or more have high temperature, it is 
suggested that the chamber be checked for cooling malfunction. 

If 1 or more batteries are cold ( < -10C) , it is suggested that 
the chamber be checked for failure of coolant control (failed ON). 

Average of the six temperature sensors > 25 degrees or < -10 degrees 
— take the average of the six temp sensors. 

5. Communication 

More than two telemetry runs missed mean either a system shutdown 
or a communication failure to the PC and will trigger an alarm 
diagnosis by NICBES. NICBES resets the fault.dat flag file to 0. 
after an alarm diagnosis, so it is ready to start again when the 
system is restarted. 

No data coming into PC. 


CURRENT DATA FILES: 

curf ( 2 , [Phase , Day_min , Cur list, Vlist, Tempi is t ] ) . 
where Curlist contains 13 SPA Currents 
Vlist contains 3 Bus Currents 

Templist contains Average Temperature for 6 Batteries 
curf(3, [Voltlist]). 

where Voltlist contains 6 sublists, one for each battery, 
with 23 Cell Voltages each. 

fault([fault flag, Type]). 


THE PREDICATES AVAILABLE IN FAULTD.PRG 





f au 1 t_ di ag— uses the current data files loaded in the database by 

prolog.ini to find the faults which cause alarms to the testbed system. 
Five classesof faults are diagnosed by calls to comm_fail/0, 
power_ci r_f ai 1 /0 , 1 oadb ank_f ai 1 /0 , chamber_f ai 1 /0 , and cell_fail/0. 

A message is written to the screen as each of the above are checked. 

If any of these five predicates discover a fault, they write a message 
to the screen. After diagnosis is done, diagnose/0 checks the 
counter 20. If more than 18 fault lines have been written in fault, 
diagnosis, it suggests that the control computer be checked for 
sending bad data. It would be unlikely for that many faults 
to occur otherwise. diagnose/0 is called automatically from 
eval_flag(l) (start. prg). 

comm_f ai l--checks f aul t ( [ _ , Type ] ) . The Data-Handler will set this 

to -1 if it detects a failure in communication, 3 missed telemetry 
runs. comm_fail/0 is called by diagnose/0. Fault diagnosis can still 
analyze last set of current data files. This may however be 
misleading so handle with caution. If no problem found, a message 
will be written to the screen so stating. 

chamber_fail — checks battery temperatures from Templist. If 

only one battery is hot ( > 25C), it suggests checking for thermal 
runaway. If more than one battery is hot, it suggests checking 
coolant control. If any battery is cold ( < -IOC), it suggests 
checking coolant control failed on. chamber_ f ai 1/0 will always 
succeed, so that other faults can be diagnosed as well. If no 
problems found, a message will be written to the screen so stating. 
chamber_f ai 1 /0 calls check_temp,/2 . It is called by f ault_diag/0 . 

check_temp(Templist , Phase) — calls check_t to recursively check the list 
of temperatures. Counters 1 and 2 are used to store the number of 
high or low values found, respectively. Counter 3 is used to keep 
track of which batteries temperatures are being checked. checkl_t/3 
is called to check highs and lows. Diagnostic messages are printed. 
check_temp/2 is called by chamber_fail/0 . 

check_t ( [ H 1 T ] ) — checks the head of the list (H) for missing data 

(-9999), or temperature out of range, prints a message if either 
condition holds, then calls itself on the tail of the list (T). 

Counter 3 keeps track of which battery is being checked. Counter 
1 keeps track of how many high values, and counter 2 counts low 
values . 

checkl_t ( Highs , Lows , Phase ) --checks number of Highs and Lows. Diagnostic 
messages are written to the screen. 

loadbank_fail — uses Phase and Vlist from curf(2,_). check_bus/2 and 
checkl_bus/2 are called to check for problems in the load bank. 
find_sum/2 is called to get the sum of the 3 bus voltages. 
Ioadbank_fail/0 always succeeds, so that other problems can also 
be diagnosed. If no problems found, a message will be written to 
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the screen so stating. load_bank/0 is called from f ault_di ag/0 . 

check_bus ( V 1 is t , Phase ) --checks bus values for missing data (-9999). 

If any bus has less than 5 A during discharge phase, failure 
of load bank or control circuitry is diagnosed. check_bus/4 
is called by loadbank_fail/0. 

checkl_bus ( Sum, Phase ) — if the sum of the bus currents is over 99A, 
failure of load bank or control circuitry is diagnosed. 

Called by loadbank_ f a i 1/0 . 

power__ci r_f ai 1 — uses Day_min, Phase and Curlist from curf(2,_), then 

calls check_cur/3 to diagnose failures in SPAs. Counter 3 is used 
to keep track of which bus is being checked. power_ci r_f ai 1 ,/0 
always succeeds, so other faults can also be diagnosed. If no 
problems found, a message will be written to the screen so stating. 
power_cir_f ail is called by fault_diag/0 . 

check_cur ( [H I T] , Phase, M) — checks data from the list of SPA currents (Curlist 
recursively, checking an item, then calling itself on the 
tail of the list. It checks for missing data (-9999), then checks 
for malfunctions in SPAs. Possible malfunctions are current 
less than 5 A in first five minutes of charge period, current 
from any SPA in discharge period, current of greater than 8 A 
from 1-SPAs (1,3,5,7,9,11), or current of greater than 16 A 
from 2-SPAs (2,4,6,8,10,12,13). check_cur/3 writes diagnostic 
messages to the screen. It is called by power_ci r_f ai 1 /0 . 

cell_fail — uses Voltlist from curf(3,_). Calls check_volt/l to check 
each battery cell for problems. If no problems found, a message 
will be written to the screen so stating. 
cell_fail/0 is called by f ault_diag/0 . 

check_volt ( Voltlist ) — sets counter 3 to keep track of which cell is being 
checked. It recurses through the 6 lists of cell voltages, one for 
each battery. The recursive function read_volt/2 is called to check 
the voltages. check_volt/2 is called from cell_fail/0. 

read_volt ( [H ! T] , N) — checks the head of the list (H) for missing 
data (-9999), voltage of 0 (diagnoses hard short), negative 
voltage (diagnoses BPRC failure), and’ voltage over 1.55 
(diagnoses overcharge). read_volt/2 then calls itself on the 
tail of the list (T). Counter 3 is used to keep track of which 
cell is being checked, for use in the failure messages. read_volt/2 
is called from check volt/1. 
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% Faultd — fault diagnosis section of NICBES. Checks five fault. 
% conditions. See fault.doc for further documentation. 


fault_diag 

ctr set (20,0) , 

n i , wr i te ( $ CHECKING COMMUNICATIONS : $) , 
comm .fail , 

n 1 , wr i te ( $ CHECK ING SPA CURRENTS : $) , 
powe r c i r _ f a i 1 , 

nl, write?$ CHECKING BUS CURRENTS: $), 
loadbank_fai 1 , 

nl, write ($ CHECKING BATTERY TEMPERATURES :$) , 
chamber _f a i.1 , 

r i .1. , wr i te ( $ CHECKING CELL VOLTAGES: $), 
cel.L_.fai 1 , 
c t r .... i. s ( 2 0 , F a u 1 1. s ) , 
i f t h e n ( F a u 1 1 s = = 0 , 

( n 1 , wr i te ( $**No cause found for alarm — check alarm 
i f t h eri ( F au 1 ts > IS, 

(nl , write ($ Many failures — check control computer 



omm_fai 1 : - 

c a 1 .1. ( f au .1 t ([_,-!] ) ) , ri 1 , 
wri te($+* Communication failure — receiving no data from sys 
wri. ts($ Fault diagnosis is limited to last set of current 
ctr _inc ( 20 ,_) . 


comm... fail write ($ No Communications Problems Found! $) . 


chamber _f ail 

call (curf (2, [Phase Templist] ) ) , 
c h ec k_ temp ( T emp 1 i s t. , Ph ase ) . 


chamber _ fail 

( retract ( flag) ; writ.e($ No Chamber Problems Found! $)). 


check_temp (Templist, Phase) : - 

ct.r_set (1 , 0) ,ctr_set (2,0) , ct.r_set (3, 1) , 

c h ec k _ t. ( Temp 1 i s t ) , 

ctr_is(l , Highs) , 

ctr_is(2,Lows) , 

c h e c k 1 .._ t ( H i g hi s , L o w s , P h a se) . 
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tern . $ ) ) ) , 

$ ) ) ) . 

tern . $ ) , n 1 , 
data f i Les .. t ' 
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check_t ( [--9999 ! T] ) : -- 

c t r ... i n c ( 3,N) , 

n 1 , write ($ Temp data item missing on battery $ ) , write (M! 
check.t (T ) . 

check _.t ( [HIT] ) 

H > 25, 
c t r ._ i n c (3, N) , 
c t r _ i n c (1 , _) , 

n .1 , w r i t e ( $ * * B a 1 1 e r y $ ) , w r i. t e ( N ) , 

write($ is overheated — $) ,write(H) , write($ degrees. $) , 
o ne_ a s s e r t ( f 1 a g ) , 
ctr_i.no (20 , _ ) , 
check. t.(T) - 

check... t ( [H ! T] ) 

H < -.10, 
c tr ... i.ric (3 , N) , 
ctr_inc(2,_) , 

ril , write ($ : +;*Bat.tery S) ,write(N) , 

write ($ is too cold at S) , write (H) , write ($ degrees. $) , 
drie_assert (f lag) , 
c tr_.inc (20 , ) , 
check. t (1 ) . 

check., t ( [H IT] ) : — 

c tr.inc (3, ... ) , 
cheek_t(T) . . 

c h ec k 1 ._ t ( _ , L o ws , _ ) 

Lows >0, 

nl , write ($**Chamber too cold — check coolant control. $), 
c tr_i.nc (20 , _ ) , 

! , fail . 

checkl.t (1 ,_, 1 ) : — 

ril , write ($**Check for thermal runaway. $) , 
ctr.inc (20 ,_) . 

checkl.t (Highs, _,_) :~ 

Highs > 2, 

n .1. , w r i t e ( $ **0verhe a t e d chamber — che c k c o o 1 a n t c o n t r o 1 . $ ) 
ctr.inc (20, _) . 
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I o a d ban k __ f a i 1 : - 

c a 1 1 ( c u r f ( 2 , [ Ph a s e , _ , , V 1 i s t. ,_] ) ) , 

ctr_set (3, 1 ) , 

chec k_bus ( VI i st , Phase ) , 

find sum(Vlist 7 Sum) , 

chec kl_bus ( Sum , Phase ) . 

1 oa db a n k _ fa i .1 : - 

(retract (flag) ; write ($ Mo Loadbank Problems Found!®)) 
check_bus ( [j , _) . 

check.bus ( [-9999 I T j , Phase ) : - 
ctr_.is(3,N) , 

ril , write ($Miss.irig data on bus $),write(N) 

ctr_inc(3,_) , 

chec k_bus ( T , Phase ) . 


chec k_ bus ( [hi J T] , Phase) : - 
Phase == o, 

H < 5, 

ctr_is(3, N) , 

ri 1 , wr i te ( $ : +' : +'Load on Bus $ ) , 

write(N) ,write($ 1 < 5 amps — load bank or control 
e t r .... i n c ( , c t r .._ i n c ( 2 0 , _ ) , 
o n e ... a s s e r t ( f 1 a g ) , 
chec k_bus ( T , Phase ) . 


failures) 


chec k_ bus ( [H ! T] , Phase) : - 
ctr_inc (3,_) , 
check_bus ( T , Phase ) . 

checkl._bus(Sum , Phase) : - 
Sum > 99 , 

nl , write ($**Bus currents > 99 amps — failure in load bank* 
circuitry .$) , 
ctr_inc(20, ) . 


con 


power_cir_fail : - 

c a 1 1. < c u r f (2, [ Ph ase , Da y min ,Curlist , ,_])), 
ctr._set.(3, 1 ) , 

check_cur (Cur list. , Phase , Day min) . 
p o w e r _ c i r _ f a i 1 : - write($ No SPA Current Problems Found!®) 
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check, ..cur ( [.] . 

c h e c k _. c u r ( [ - 9 9 9 9 ! T ] , P h a s e , D a y m i n ) : - 
c t. r i n c ( 3 , i... } , 

ril., write ($Missing data from SPA $ ) , wr i te ( L ) , n 1 , 
chec k ..cur ( T , Phase , Day m in ) . 

check_cur ( [HIT] , F'hase , Daymin ) : - 

etr_is(3,L) , 

Phase == 1, 

Day min < 7, 

H < 5, 

H > -50, 
c tr_inc (3 , _ ) , 

nl , wr i te ($*+Cur rent .low from SPA-$) , write(L), 
write($ in first 5 minutes of day — $),nl, 
write($ failure of power supply or control.®), 
c tr_inc (20 , _ ) , 
check_cur (T , Phase , Daymiri ) . 

check_cur ( [H IT], Phase , Daymin ) : - 
Phase = = 0, 

H > 5, 

c tr_inc (3, L ) , 

nl , wr i t e ( $ * * C u r r e n t. f r om S P A - $ ) , wri t e ( L ) , 
write($ during night — control circuitry failure®.), 
c t r inc ( 20 , _ ) , 
check _cur ( T , Phase , Daymin ) . 

c h e c k c u r ( [HI T] , P h a s e , D a y m i n ) : - 
ctr._.i.s(3, L) , 

L =\= 13, • 

I< is L mod 2 , 

K == 1, 

H > - 8, 

ctr_inc (3,_.) , _ . 

nl , write (S't^Cur rent from SPA-$) , write(L), wri te ($ > - H amps . $ ) , ri ... , 
write($ Check for failure in control limit circuitry .$) , 
dt.r_.inc (20 , _) , 
check_cur (T , Phase , Daymin ) . 
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heck ..car ( [H 5 Tj , Phase , Daymin) = “ 
c 1 r .... i s v 3 , L ) , 

K is L mod i. 5 _ ^ 

( « r.: - O ; L. = = i } ’ 

H / -• 1 1 < 7 

M r wri?p(ii*curr 8 nt from SPA-*) , u r i t«(L),«r ite<* > = 
nr i te ( * Check current limit circuitry..), 

c t r ___ .i. n c v 20 . 

chec k_cur ( T , Phase , Day mm ) - 

check.cur ( [H IT] , Phase , Daymin ) : ~ 

rtr inc(3,_)i . , 

chec k_cur ( T , Phase , Daymin ) - 

cel 1 fail 

"call ( cur f( 3, Volt list) ) , 

c t r s e t C* 1 7 1 ) * 
c h e c k _ v o 1 1 ( V o 1 1 1 1 s t ) - 

C9ll ~(*etract(flafl) ; nrite(* No Cell Failures Found!*)) 

c h © c k _volt([j) - 

C h & c k ... V o 1 1 ( C H ! T ] ) s- 
ctr .inc ( 1 * N ) , 
r; ■{"I'" Sei- \ O <? I /* i 
Y' 0 cl d .... V O JL "t C H , N) 7 

r: h ec k_v o .1 1 (T) - 

v" e<::l d V 0 1 1- s i_ -J •? — > - 

r ead_vo 1 1 < [-9999 ! T ] , N ) : ~ 

nl , write ($Missing data on cell $) ,w»ite(L , 
w r i te ( % o f ba t te r y $) , write (N ) 7 
r e a d _ vo 1 1 ( T , N ) - 

read„vol t ( [0 ! T] , N ) : - 

ctr_inc(3,L) , t i n rP H $) ,write(L) 

n 1 , wr i te ( $**Suspec t hard shor t ~r. ofeu 

w r i t e ( $ o f b a 1 1 e r y S ) , write (N, , 
ct.r_.inc (20, ) ■? 
one asser t (. i’ iaQ ,■ ■> 
r ead_vol t- C T •> • ^ ) - 


* 
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r e a d _ v o .1 1 ( [ H ! T j , N ) 

~H < 0, 

r: t r iric (3 , L ) ■> . H ,• ■. ■» 

nl,write($**BPRC failed in cell $;,wixteu/, 

w r i t e ( $ o f b a 1 1 e r y $ ) , w r i t e (. N ) , 
tf ... i nr (20 , _J , 
ope ;'riSS 0 r t ( f IclQ • f 
read_vol t. ( T ? N ) - 

r e a d _ v o 1 1 ( [Hi T J ,N) : - 
H > 1 - 55 , 

nllwrite($**Volta9e too high in cell $),writ evU, 

write($ of battery t) , write(N) ,n . ’ charge rate.S) 

write($ check overcharge or too hxg <- 

ctr _iric ( 20 , _) , 
one assert (flag) ■> 

read_volt (T 5 N) - 


r ead _vo.l. t ( [H ! T J , N ) : ~ 
ctr _inc (3 , L) , 
r e a d . _vol t (T , N) . 
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Status. prg has the STATUS portion of NICBES. Selecting Battery Status 
from the Main Menu for battery n (n is 1 through 6) will trigger the 
checking' of battery n for reconditioning, temperature, workload, 
charging scheme, and divergence. Messages will be printed on the screen. 

The STATUS section uses averages of data from some of the showf (n ) . dats 
in its analysis. If nothing is wrong with the battery, a message 
is written to this effect. Problems with the battery often 

produce messages to use "Plots and Graphs" or "Advice" for more information. 


The following conditions are examined: 

Condition Status — 

Temperature average over last orbit < 0, Cold 

and over last 12 orbits < 0 


Temperature avg. over last orbit > 11, Hot or possible overcharging 

and over last 12 orbits > 10 

Temperature avg. over last orbit > 11, Possible overcharging 

and over last 12 orbits < 10 


AHO average over last 12 orbits < 3 Underwork 

(3.0 Ampere-hours), and avg. EOD 
voltage < 27 

AHO average over last 12 orbits > 14 Overwork 

(14.0 Ampere-hours) 


Recharge ratio average < 1.020 Insufficient charge 

and avg. high-charge voltage < 32.5V 

High-charge voltage > 33.8 Possible overcharge 


Divergence average > ,8V between High divergence 

high and low cells 


High cell minus average greater than Too many cells to 

(avg. minus low cells) + 1 low values 


PREDICATES AVAILABLE IN STATUS. PRG 

battery_status (Bat ) — First, battery_status/I checks for reconditioning. 
If the battery is being reconditioned, status analysis stops 
at that point, since data will be misleading. Otherwise, four 
procedures are called, temp_status/l , work_status/ 1 , 
charge_status/l , and div_status/l . Each of these procedures 







writes messages to the screen and increments counter 15 if 
problems are found. If counter 15 is still zero, a message is 
written to the screen that the battery seems healthy. 
battery_status/l is called from execute(2 , Bat ) ( s t ar t . pr g ) . 

temp_status ( Bat ) — showf(8,_) (temperature over last orbit) and showf(9,_) 
(average temperature over last 12 orbits) are used in analysis. 
get_data/5 is called to get the appropriate battery data and then to 
find the its average. Finally, check_ t erap/3 is called to analyze the 
data. temp_status/l will always succeed, even if no problem is found, 
so that status analysis can continue. If no problem was found a 
message will be printed to the screen. 
temp_status/l is called from b a t t ery_ s t at us/ 1 . 

work_status ( Bat ) — showf(13,_) (AHO for last 12 orbits) and showf(l,_) 

(EOD battery voltage over last 12 orbits) is used to analyze workload 
status. get_data/5 is called to get the appropriate battery data and 
then to find the its average. Analysis is done by calling check_work/3 
find_avg/2 is called to find averages. work_status/l will always 
succeed, so that status analysis can continue. If no problem was found 
a message will be printed to the screen. 

It is called from battery_status/l . 

charge_status (Bat ) — showf(3,_) (recharge ratio) and showf(2,_) (high 
cell voltage during charge) is used to analyze charge status. 
get_data/5 is called to get the appropriate battery data and then to 
find the its average. Analysis is done by calling check_charge/3 . 
charge_status/l will always succeed. If no problem was found a 
message will be printed to the screen. It is called from 
batter y_stat us /l. 

get_data ( Bat , N1 , N2 , Avgl , Avg2 ) — N1 and N2 reference the show files, Bat 
is the selected battery. get_list/3 is called to get the 
data list for the appropriate show file and battery. find_avg/2 
is then called to calculate the average for each data list. 

Cut is used to prevent backtracking. get_data/5 is called by 
temp_status/ 1 , work_s t at us/ 1 , charge_ s t atus/ 1 . 

div_status (Bat ) — showf(4,_) (EOD divergence) is used to analyze divergence. 
get_list/3 is used to get the appropriate battery data, find_div/5 is 
called to find highs, lows and averages, then eval_div is used to 
analyze the divergence. div_status/l will always succeed. If no 
problem was found a message will be printed to the screen. 

It is called from battery_status/'l . 

check_recond(L.ist)--checks the single value passed it by batter y_ status/1. 

If this value is 1, the battery is in reconditioning, and a message 
is printed on the screen. No further status analysis is done. 

check_temp(Oavg, Avgl2 , Bat ) — If average temperature in last orbit and 
average temperature over last 12 orbits are both below zero, 
a message that the battery is cold is given. If average temperature 
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over the last orbit is over 11, and average over last 12 
orbits is over 10, a message is given that the battery is hot. 
If average over last orbit is greater than 11, and average over 
last 12 orbits is less than 10, a message is given to check 
chamber or check for overcharging. check_temp is called by 
temp_status/l . 


check_work ( Bat , Ahoavg , Eodv ) — analyzes the battery for overwork or underwork. 
If the average AHO is less than 3 ampere-hours, and EOD voltage 
is less than 27 V, a message that the battery may have memory 
effect is given. If the average AHO is greater than 14 ampere- 
hours, a message is given that the battery is overworked. check_work/3 
is called from work_s t a t us/ 1 . 

check_charge ( Rr , Cavg , Bat ) — analyzes the battery for overcharging or 

undercharging. If recharge ratio is greater than 1.020, and high 
battery voltage during charge is less than 32.5, a message is given 
that the charging scheme may be insufficient. If in-charge high 
voltage is greater than 34 volts, a message is given that 
possible dangerous overcharging may be occurring. In-charge high 
voltage average greater than 33.6 volts triggers a message to 
check for overcharging. check_charge/3 is called from charge_s tatus/ 3 . 


f ind_div ( List , Avg, Div , Highs , Lows ) — Given a List taken from 

showf(4,_) or showf(5,_), find_div/5 calls break_list/4 to break 

the list into a list of high values, a list of low values, and a 

list of average values for cell voltages. find_avg/2 is then 

called on each of these lists to find the average high reading (Highs), 

low reading (Lows) and average reading (Avg) over the 12 orbits. 

Div is instantiated to Highs minus Lows. find_div/5 is called by 
div_s tatus/ 1 . 


eval_div (Div , Avg, Highs , Lows , Bat ) — evaluates two measures of 

divergence, using variables passed out from find_div. eval_div/5 
writes a message that divergence is high if Div is greater than 
0.8. If Avg is closer to Lows than to Highs by a noticeable margin 
(1), a message is given that too many cells are migrating to low 
values. eval_div is called by di v_status/l . 
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% This section has status portion of NICBES 

% status (Bat) examines the battery status from five viewpoints 
% See Status.doc for documentation 


b a 1 1 e r y _ s t a t u s ( B a t ) : - 

call ( recond ( Recorid) ) , 

fin d _ n t h ( Recond, 1 , B a t , L i s t ) , 

w r i te ( $CHEC K I MG RECOND I T I ONI NG :$) , 

c h e c k _ r e c o n d ( L. i s t ) . 


b a 1 1 e r y _ s t a t u s(Bat) : - 

wr.ite($ Battery is not being reconditioned! 

ctr_set.( 15,0) , 

t e m p s t. a t u s ( Bat) , 

wor k s t a t. u s ( B a t ) , 

char ge_status (Bat) , 

d i v s t. a t u s ( B a t ) , 

c t r i. s ( 1 5 , A ) , 

ifthen((A < 1), 

(nl , write ( $ * * Battery seems heal t h y ; no obv i o u s 


'inp 


status (Bat) : - 

n I , wr i te ( $ CHECK. TNG TEMPERATURE: $) , 


9 e t _ d a t a ( B a t , 8 , 9 , A v g , A v g 1 ) , 


c he c k _ t e m p ( A v g , A v g 1 , B a t ) . 


$) , 


problems .$))). 


ternp_status (Bat ) write(S No Battery Temperature Problems ! $ ) , n j 


w o r k ... s t a t u s ( B a t ) : - 

” ri .1 , wr i te ( SCHECK ING WORKLOAD : $ ) , 
get. .data (Bat , 13, 1 , Avg , Avgl ) , 
check., wor k (Bat , Avg , Avgl ) . 

wor k_status (Bat ) write($ No Battery Workload Problems! $),nl. 


c h a r g e ._ s t a t. u s ( Bat) : - 

ri 1 , w r i te ( SCHEC l< I NG CHARGE : $ ) , 
get...data (Bat ,3,2, Avg , Avgl ) , 
check_.charge (Avg , Avgl , Bat ) . 


c h a r g e _ s t a t u s ( B a t ) 


wr i te (S 


No Battery Charging Problems 


: ) , n 1 . 


get_data ( Bat , Nl , N2 , Avgl , Avg2) : - 
get_list (Bat. , Nl , List.) , 
f ,ind_avg (List , Avgl ) , 
get_list(Bat,N2,Listl ) , 
f ind_avg (List! , Avg2) , ! . 
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d j. v _ s t a t. u s ( B a t ) : " 

w r i t e ( $ C H E C KING DIVERGENCE :S) , 

9 e t _ I i s t ( B a 1 7 4 , L i s t ) , 

find div (List , Avg , Div , Highs , Lows ) , 

eval d i v ( Div , Avg , Highs , Lows , Bat ) - 

div status(Bat) writ.e($ No Battery Divergence Problems ! $),n 


c h e c k ... r e c o n d ( [ 1 ] ) :* 
ri 1 , 

wr i te ($**Ba t ter y is in reconditioning. $) ,nl, 
writ.e($**No further Status, as analysis would be mi 


ading . **$ ) - 


c heck... t e m p ( 0 a v g , A v g 1 2 , B a t ) : - 

Oavg < 0, 

Avg 12 < 0, 
c t r _inc(15,_) , n 1 , 

write($**Battery is cold — check chambers) , nl , 

write($ Select Plots and graphs from Main Menu for Battery $) , 
write (Bat), nl, 

write ($ with choice 8 from Graphics Menu $),n.l. 


chec k_ temp ( Oavg , Avgl 2 , Ba t ) 


Oavg > .1.1, Avgl2 > 10, 
c t r _ i ri c ( 1 5 , _ ) , n 1 , 

write (S : + :; * ; Battery is hot — check chamber or overcharging$) ,n.i, 
write(S Select Plots and Graphs from Main Menu for Battery $) , 


write (Bat) 
wr i te ( $ 
wr i te ( S 
n .1 . 


, ri .1 , 

with choice 8 from Graphics Menu to see 
and with choice 2 from Graphics Menu for 


temp , $ ) , nl , 
o v e rcharg i rig 


c h ec k_ temp ( 0a vg , A vg 1 2 , B a t ) : - 
Oavg > 11, Avgl2 < 10, 
c t r _ i n c ( 1 5 , _ ) , ri 1 , 

w r i t e ( $ * * B attery is ho t t h is cycle — c h e c k f o r o v e r c h a r g e $ ) , n l , 
write($ Select. Plots and Graphs- from Main Menu for Battery S) , 

write(Bat) , ril , ... 

wr i te ( $ with choice 2 from Graphics Menu for overchargi.ng$) , nl . 

c h e c k w o r k ( B a t , Ah o a v g , E o d v ) : - 
Ah oavg < 3, 

Eodv < 27, 
c t r ... i r i c ( 1 5 , _ ) , ri 1 , 

write($**Batt.ery may have memory effect .$) ,ril , 

write (S Select Plots and Menus from Main Menu for Battery S) , 
write (Bat) ,nl, . . 

write ($ with choice 1 from Graphics Menu for plot of EOD voi iciv^e . $ 
write($ Select Advice from Main Menu for Battery S) , 

write (Bat) ,nl , ... 

write($ with choice from Advice Menu for reconditioning. $) ,ril . 
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c h e c k __ w o r k ( B a t , A h o a vg,E o d v ) : - 
Ahoavg > 14, 
c tr _inc ( 15 ,_) , n 1 , 

w r i t e ( $ * * B a 1. 1 e r y m a y b e o v e r w o r k e d . $ ) , n 1 , 
w r i t e ( $ S e 1 e c t A d v i c e f r o m M a i n H e ri u f o r B a 1 1 e ry $ ) , 
w r . ! t ■ © ( B a !~ ) , f i .1. , 

wri te($ with choice from Advice Menu for workload .® ) ,ni 

chec k charge ( Rr , Cavg , Bat ) : -* 

Rr < 1.020, 

Cavg < 32.5, 
c t r _ i n c ( 1 5 , _ ) , n 1 , 

write($**Battery charging scheme may be .insufficient. $),n 
write($ Select Advice from Main Menu for Battery $) , 
write (Bat ), n.l , 

write($ with choice from Advice Menu for workload.®) ,nl 

c h ec k_c h a r ge ( _ , Cavg , B a t ) : - 
Cavg > 34 , 
c t r i n c ( 1 5 ,_) , ri 1 , 

wri te($**Possible dangerous overcharging . **$) ,nl , 
writeCS Select Advice from Main Menu for Battery $) , 
wri te(Bat) , ril , 

write($ with choice from Advice Menu for charging.®) ,nl 


chec k_char ge ( _ , Cavg , Bat ) : - 

<3 V ) »."> O - 6 •} 

c It r i n c ( 1 5 ? _ ) ,n .1 

w r i t e ( $ * •+: c h e c k f o r o v e r c h a r g i n g . •+ * ® ) , n 1 , 
w r i t e ( $ S e 1 e c t Ad v i c e from M a i n M e n u f o r B a 1 1 ery $) , 
wri te ( Ba t ) , n 1 , 

write($ with choice from Advice Menu for charging . ® ), ril 

f ind_di v (List , Avg , Div , Highs , Lows ) : - 
b r e a k _ 1 i s t ( L. 1 s t , L i s t .1 , L i st.2,12) , 
break_list.(List2, List3, List.4, 12) , 
f irid_avg (Listl , Highs) , 
firid_avg(List3,Lows) , 
f irid_avg (List4 , Avg) , 

Div is Highs - Lows, ! . 

eval_di v (Div , Avg , Highs , Lows , Bat ) : - 

A is Highs - Avg, 

B is Avg - Lows, 

C is B + 1, 

C < A , 

n 1 , w r i t e ( $ * Too m a n y c e 1 1 s m i. g r ating t o 1 o w v a lues. *$) , 
c tr_inc ( 15 , _ ) , fail . 
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d .1 v ( D i v , A v g , H i g h s , Lows , B a t ) : - 

Div > 0.8, 

ril , write ($**Divergerice is high, may need reconditioning. 
nl,write($ Select Advice from Main Menu for Battery $) , 
w r x t e ( B a t ) , n i , 

wr i te($ with choice from Advice Menu for reconditioning 

c t r _ i n c ( 15, _) - 


3-7 -? 


Advice. prg contains the ADVICE section of NICBES which goes into 
further detail on three subjects: whether a battery needs reconditioning, 

changes in charging scheme, or changes in workload. Select Advice from 
the Main Menu for battery n brings up a menu of three choices for the user. 
The ADVICE section finds trends in many of its data files, and uses them, 
along with averages, for analysis in greater depth. Messages to 
the user give trends and average values, as well as advice. 

ADVICE section, using trends in voltage, recharge ratio, divergence 
and temperature — 

AHO average (Ahoavg) is used for advice on workload. It is considered to 
be overly high above 14, high above 9, average around 7 and low below 3. 

EOD pressure (Eodpress) is considered to be high above 100, and fairly high 
above 90. EOC pressure (Eocpress) is considered to be high above 110, and 
fairly high above 100. 

Divergence (Div) is considered to be high above 0.8V (difference 
between high and low cells). 


PREDICATES AVAILABLE IN ADVICE. PRG 

advice ( Bat , Choice ) --with choice 1, advice/2 calls get_one/6 to get data, 
and recond/5 to analyze it. With choice 2, advice/2 calls 
get_two/7 to get data and chgchg/6 to analyze it. With choice 3, 
advice/2 calls get_three/4 to get data and chgwrk./3 to analyze it. 
advice/2 is called by more/2 (start. prg). 

get_two(Tren d2, Trend3, Div, Trend9, Eodpress, Eocpress, Bat) — The 

hand-crafted deviation factors are called ( Dev2 , Dev3 , Dev5 , Dev9 ) . 

Then get_list/3 is called to read data from showf(2,_) 

(in-charge cell high voltage), showf(3,_) (recharge ratio), 
showf(5,_) (in-charge divergence), showf(9,_) (temperature), 
and showf(10,_) (cell pressure). Trend analysis is done on voltage, 
divergence and temperature by calling trend_analys is . find_div/'5 is 
called to find average divergence . find_avg is called to find average 
cell pressure EOC (Eocpress) and EOD (Eodpress). get_two is called by 
advice/2 . 

get_three(Trendl, Trend4, Ahoavg, Bat) — Hand-crafted deviation 

factors are called (Devl,Dev4). Then get_list/3 is called to read 
data from showf(l,_) (EOD battery voltage), showf(4,_) (EOD 
divergence) and Statfl,_) (AHO). 

trend_analysis/3 is called to derive Trendl of voltage. disp_trend/3 
is called to find divergence trend. find_avg/2 is called to find 
average AHO. get_three is called by advice/2. 





get_one(Trendl,Trend2,Trend3,Trend4,Trend5,Bat) — Hand-crafted 

deviation factors are called ( Dev 1 , Dev2 , Dev3 , Dev4 , Dev5 ) . Then 
get_list/3 is called to read data from showf(l,_) (EOD 
battery voltage), showf(2,dat (in-charge cell high voltage), 
showf(3,_) (recharge ratio), showf(4,_) (EOD divergence) and 
showf(5,_) (in-charge divergence). 

trend_ analysis/ 3 is called for voltages and recharge ratio. 
disp_trend/3 is called to find trends for divergence data. 
get_one is called from advice/2. 

disp_trend(List, Trend, Dev) — break_list/4 is called to break the 

divergence List into highs (Listl), lows (List3), and averages 
(L.ist4). get_disp/3 is called to figure the divergence list (Divlist). 
Then trend_analys is/3 is called on Divlist. disp_trend/3 is 
called by get_three/4 and get_one/6. 

get_disp ( List 1 , List2 , Lis t3 ) — is a recursive function to find the 

differences between corresponding values of two lists, and put the 
differences in a third list. It is used to find the divergence 
between high-voltage cells and low-voltage cells in the data from 
12 orbits. Missing data (-9999) in either of the comparands means 
that -9999 is written to the list of divergence. Comparison is 
done of the heads of the two lists, and the result is put in the 
head of the third list. Then get_disp/3 is called recursively on 
the tails of the lists. Halting condition for the recursion is the 
empty list. get_disp is called by disp_trend/3 . 

trend_analys is ( List , Trend , Dev ) — analyzes trends in data using two 

simple mathematical functions. find_wl/2 is called, which weights 
recent values more heavily, and sums the weighted values. Then 
find_w2/2 is called, which weights all the values equally, with 
the weight ((n+l)/2) chosen so that for a list of constant values, 
find_wl and find_w2 give the same sum. Then eval_trend/3 is called 
using the difference of the two weight functions (Diff), and the 
hand-crafted deviation factor (see dev/6 below) appropriate to the 
set of data being analyzed, and returns the Trend. trend_analysis/3 
is called by get_two/7, get_three/6, and get_one/4). 

f ind_wl ( Weightl , List ) — given a list of n values, find_wl/2 

calls multwt/2 to multiply the first by 1, the second value by 2, 
the third by 3, and so on till the nth value is multiplied by n. 
find_wl/2 uses counter 12 to keep track of the number of 
items in the list. Weighttl is the sum of the individual products. 
find_wl is called by t rend_ana 1 ys is./3 . 

multwt ( Wt , [H ! T] ) — checks for missing data (-99), then calls 

itself recursively on the tail of the list ( T) until the empty list is 
reached. Each recursive call that does not find a missing data item 
increments the counter. Then as multwt/2 climbs out of the recursion, 
the value of the counter is multiplied by the head of the 
list (H) at that time, and is added to the weight, which has 
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been is initialized to 0 for the empty list. multwt/2 is 
called by find_wl/2. 

f.ind_w2 ( Wei ght2 , List ) — initializes a counter, then calls find_sum 

on List to find the sum of the values. Then a weight factor is 
figured which will make Weight2 (from find_w2) equal to Weightl 
(from find_wl) if the List contains all one constant value. This 
weight factor is figured by (n+l)/2 for a list with n elements, 
and is multiplied by the sum. find_w2 is called by t r end_ analys i s / 3 . 

eval_t rend ( Dev , D i f f , Trend ) — takes the deviation factor, and the 

difference (Diff) between Weightl and Weight2, and determines a 
trend. If the absolute value of Diff is less than Dev/3, the 
trend is none. For positive Diff, Diff greater than Dev/3 and 
less than Dev is considered to indicate the trend is s 1 i ght 1 y_up . 

Diff greater than Dev but less than 2 times Dev indicates the 
trend is up, and Diff greater than 2 times Dev indicates the 
trend is strongly_up. Similarly, negative Diff is s 1 i ght ly_down 
for Diff less than -Dev/3, down for Diff less than -Dev, and 
st rongly_down for Diff less than 2 times Dev. Prolog syntax is 
used to advantage in this code, with atoms giving the trends. 
eval_trend/3 is called by trend_analys is/3 . 

recond ( T1 , T2 , T3 , T4 , T5 ) — writes the trends (EOD voltage, in-charge 
high voltage, recharge ratio, EOD divergence and in-charge 
divergence) on the screen, and calls tell_l/5 to analyze the 
trends to determine whether reconditioning is advisable, 
recon d/5 is called by advice/2. 

chg'chg(T2,T3,Div,T9,Eodpress,Eocpress) — writes the trends (in-charge 
voltage, recharge ratio, and temperature) on the screen, along 
with the averages (in-charge divergence, EOD pressure and EOC 
pressure). Then tell_2/6 is called to advise as to whether the 
charging scheme should be changed, chgchg/6 is called by advice/2. 

chgwrk ( T1 , T4 , Ahoavg) — writes EOD voltage and EOD divergence 

trends to the screen, along with the AHO average. Then tell_3/3 
is called to advise if the workload should be changed. chgwrk/3 
is called by advice/2. 

tell_3 ( T1 , T4 , Ahoavg) — Six conditions are used by tell_3 to advise 
as to workload change. 

1. EOD voltage down, EOD divergence down, AHO avg. less than 5 

— loads probably too -light. 

2. EOD voltage strongly down, EOD divergence strongly down, AHO 
average less than 5 — losing capacity due to memory effect, 

try heavier load. 

3. AHO average greater than 14 — possibly destructive overwork. 

4. AHO average greater than 9 — heavy workload. 

5. EOD voltage down or strongly down, EOD divergence down or 
strongly down, AHO average above 6 — losing capacity. Consider" 

reconditioning. 
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6. none of the above — see no need to change workload. 

te.ll_2(T2,T3,Div,T9,Eodpress,Eocpress) — Nine conditions are used by 
tell_2/'6 to determine if the battery should have its charging' 
scheme changed. 

1. Temperature strongly up, EOD pressure greater than 100, 

EOC pressure greater than 110 — overcharging causing high 

pressure . 

2. In-charge voltage strongly up --check for overcharging. 

3. EOD pressure greater than 80, EOC pressure greater than 
100 — charge rate may be too high. 

4. EOD voltage strongly down, in-charge voltage strongly down 

— undercharging may be seriously affecting capacity. 

5. EOD voltage down, in-charge voltage down — appears to be 

undercharged. 

6. EOD voltage strongly down or down, in-charge voltage strongly 
down or down — undercharging may be causing loss of voltage. 

7. EOD voltage down, in-charge divergence greater than 0.8 

--undercharging may be causing loss of voltage. 

8. EOD voltage strongly down, in-charge divergence greater 

than 8 — undercharging may be causing serious loss of voltage. 

9. none of the above — see no need to change charge scheme. 

tel 1_ 1 ( T1 , T2 , T3 , T4 , T5 ) — uses 11 conditions to determine if a 
battery needs to be reconditioned. 

1. EOD voltage strongly down, in-charge divergence strongly up 

— reconditioning advised. 

2. EOD voltage strongly down, in-charge voltage strongly down 

— recondition soon. 

3. EOD divergence strongly up, in-charge divergence strongly 
up — strongly rising divergence indicates reconditioning. 

4. Voltage trends strongly down or down, recharge ratio and 
divergence trends strongly up or up — reconditioning is 

indicated by all five trends. 

5. Recharge ratio strongly up — reconditioning indicated. 

6. In-charge voltage is down or strongly down, in-charge 
divergence is up or strongly up — reconditioning recommended. 

7. EOD voltage down or strongly down, in-charge divergence 
up or strongly up — reconditioning indicated. 

8. Both divergences up — consider reconditioning. 

9. Recharge ratio up — consider reconditioning. 

10. Both voltages down — consider reconditioning. 

11. none of the above — see no need to recondition at present. 

dev ( Dev 1 , Dev2 , Dev3 , Dev4 , Dev5 , Dev9 ) --Data structure, listed as a fact, 
which contains the hand-crafted deviation factors. Each one 
is associated with a data file. 

Devl=2 - showf(l,_) - EOD battery voltage 

Dev2=2 - showf(2,_) - in-charge high battery voltage 

Dev3=0.02 - showf(3,_) - recharge ratio 

Dev4=0.15 - showf(4,_) - EOD cell voltage divergence 

Dev5=0.09 - showf(5,_) - in-charge cell divergence 


b'4 -fV 



Dev9=13 - showf(9,_) - battery temperature 

These values are used in the weighting scheme to analyze trends. 

To INCREASE sensitivity to trends', lower the values. These need 
not be integers, but must be greater than 0. 

To change these values, put advice. prg in your favorite editor and 
revise the data structure dev/5. No other changes are necessary. 
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'V. advice . prg derives trends to 9 
% battery should be reconditions 
% w o r k 1 o a d c: h a n g e d . S e e a d v 1 c e - 


ive further analysis on whether a 
d, charging sheme changed or 
doc for documentation . 


adv .1 c e ( B a t , 1 ) 

write($** ADVICE ON RECONDITIONING BATTERY: **$) ,nl, 
g © t. _ o ne(Trendl, Trend 2 , Trend3, Trend4, Trends, B a t ) , 
recond ( Trend! , Trend?. , Trends , T rend4 , T rends ) . 

a d v .1 c e ( B a t , 2 ) : - ' , 

write ($*=♦= ADVICE ON CHANGING CHARGING REGIME OF BATTERY: ** $',n 
get two ( T rend? , T r erid3 , Di v , T r end9 , Eodpress , Eocpr ess , Bat ) , 
chgchg (Trend? , Trends , Div , Trend'? , Eodpress , Eocpr ess ) . 

a d v i c e ( B a t ,3) . 

wri te ($** ADVICE ON CHANGING WORKLOAD OF BATTERY: **$> ,nl, 

get_ three (T reridl , T rerid4 , Ahoavg , Bat ) , 
c h gw r k ( T r eri d 1 , T r eri d4 , Ah oa vg ) . 

get_two( Trend? , Trends , Div , Trend? , Eodpress , Eocpr ess , Bat ) : - 

| c a 1 1 ( dev ( , Dev? , DevS Dev? ) ) , 

' ge t _ 1 i s t ( B a t , 2 , L i s 1 2 ) , 

trend_analysis (List? , T rend? , Dev? ) , 

g e t .... 1 i s t ( B a t , 3 , Lists) , 

trend _ana lysis ( Lists , T rends , DevS ) , 
g e t _ 1 i s t ( B a t , 5 , L i s 1 5 ) , 

r in d d 1 v ( L .1 s i S , 7 Div -i 7 ) 7 

get_Iist(Bat.,9,List?) , 
t r ert d_ari a 1 y s i s ( L i s 1 9 , Trend? , Dev9) , 
ge t 1 i s t ( B a t ,10, L i s t .1 0 ) , 
b r e a k _ 1 i. st(L i. s 1 1 0 , D .1 i s t. , C 1 i s t. , 2 3 ) , 
find_avg(Dlist, Eodpress) , 
find„avg (Clist , Eocpress ) . 

get._ three ( Trend 1 , T rerid4 , Ahoavg , Bat ) :- 
call (dev (Devi Dev4, , 
g e t _ 1 i s t ( B a t , 1 , L i s t i ) , 
trend._analys.is (Listl , T rendl , Devi) , 
get_list. (Bat , 4 , List4) , 
di s p t r e n d ( L 1 s 1 4 , T r e n d 4 , D e v4) , 
get_Iist (Bat , 13, Listx ) , 
f ,i r 1 d . a v g ( L .1. s t. x , A h oa v g ) . 
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one ( Tr endl , T r end2 , Trend3 , T rerid4 , T reridb , Bat . ) : ~ 

c a 1 1 ( dev < Dev .1 , Dev2 , Dev3 , Dev4 , Dev 5 , ) ) , 
cj e t . 1 i s t ( b a t , 1 , L i s t .1 ) , 


y p fj 

an a 1 y s i s ( L. l s t .1. 

, T reridl , Iievl ) , 

;”j r. L 1 

s t ( B a t , 2 , L i s 1:2 

J , 

tren d 

n ct 1 y s i. s ( L i. s t.2. 

,Trend2, Dev2) , 

ae t j. i 

st (Bat , 3 , L.ist.3 

i 

f •) 

trend . 

a n a 1 y s i s ( L i s 1 3 

, T ren do , Dev 3 ) , 


g e t _ I. i s t ( B a t ,4, Li s 1 4 ) , 

di sp_ trend ( Li. st.4 , T r end4 , Dev 4 ) , 

g e t _ 1 i s t ( B a t , 5 , L i s 1 5) , 

d i s p „ t r e ri d ( L i s 1 5 , T r e n d 5 , D e v 5 ) - 


d i sp _ t ren d ( L. i s t , Trend, Dev ) : - 

break_list (List. , Listl , L. i. s 1.. o ,12) , 
break_l ist. (List.2 , List3 , List4 , 12) , 
get d isp ( L i st.1 , List3, Di v list.) , 
t r eri d ... an a .1 y s i s ( D i v 1 i. s t , T r eri d , Dev ) . 


get._d.isp(. [],[],[])- 

get_disp( [] , List , [] ) . 

g e t d i s p ( L i s t , [ ] , [ ] ) . 

get disp ( |.H1 ! T 1 j , [H2 ! T2] , [Hr ! ! 4] .) 
(HI == —9999 ; H2. -9999) , 

Hr = -9999 , 
ge t ... d i sp ( T1 ,T2,T4) . 


get._di.sp ( [HI ! Tl] , [H2 ! T2] , [Hr ! Tr] ) :~ 

Hr is HI - H2 , 
get_disp (Tl , T2 , T r ) . 

trend.., .analysis (List, Trend, Dev) : — 
f ,i ri d ._ w 1 ( W e i g h 1 1 , L. i s t ) , 
f i. rid~w2 (Wei gh t2 , Li s t ) , 

Diff is Weigh tl - Weigh t2, 
eval trend (Dev , Diff , Trend) . 


f 1 n d w 1 ( W e i g h 1 1 , L ist.) 
"ctr_set. (12,1) , 
mu 1 1 w t ( We i gh 1 1 , L i s t ) . 
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fjiultwt. (Wt , [] ) 

Wt is 0. 

multwttWt, CHIT] J tr inc(12,A) ) , 

\ f then vH ■■■■■■> 

multwt (Temp , i ) , 

B i f A ? is Temp ’ Wt 1 

if then else v.H — '•••••’ 

f i n d _w2 ( We i gh t2 , L i '-=> t 
c; c se t ( 10 •> 0 ) ■) 
f ind_sum (List , sum ) , 
c t r _ i s ( 1 0 , A > , 

R is A + li 
C is E/2.0, 

Weights is C * bum. 

trend(Dev,Diff , Trend) 
a is abs(Diff), 

A < Dev , 

Trend = none. 

eV al _trend(Dev, Diff , Trend) 

"Dif f > 0, 

Diff < (2 * Dev), 

Trend = slightly_up- 

eval_ trend (Dev , Diff , Trend) : ~ 


Ditt ) 
Diff : 
Trend 


( ■*, "4 s - 0 ft- 'v ) ■? 

up - 


p v a } t r e ri d ( D e v , D i ff,Tt e n d ) 

’Diff > 0 , 

Diff > (3 * Dev), 

Trend = strongly _up - 

eval_ trend (Dev , Diff , Tf end) 
Diff < 0, 

A is -2 * Dev, 

D i. f f > A , 

Trend = si ight.ly_dowri . 

v a 1 _ t r e n d ( Dev, D i f f , Trend) :■ 
Diff < 

A is -3 •+' Dev, 

Diff > A ' 


Trend = down. 


Temp) - 
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e v a .1 ... t r e ri d ( D e v , I) i f f , 7 r e n d ) : - 

Diff < 0, 

A is -3 * Dev, 

Diff =< A, 

T r e n d = s tron g 1 y _ d o w n . 

r e eo rid ( 71 , 72 , 73 , 74 , 75 ) : - 

w ri t e ( $ E 0 D v o 1 1 a g e t r e nd i s $ ) ,wr i t e ( T 1 ) , n 1 , 

wr i te ( $In-char ge voltage trend is $) , wr ite (72) ,ril , 

wr.ite($Recharge ratio trend is $), write (T3) , rtl , 

wr i t e ( $ E 0 D diverge n c e trend i s $ ) , wr i t e < T 4 ) , n 1 , 

wr i te ($In-charge divergence trend is $) , write (T5) ,nl , 

te 1 .1. _ 1 ( 7 1 , 72 , T3 , 7 4 , 75 ) . 

chgchg (72, 73 , Div , 79 , Eodpress , Eocpress ) : - 

wr i te ( $In-charge voltage trend is $ ) , wr i te (72 ) , n 1 , 
wr i te (SRecharge ratio trend is $) ,write(73) ,n.l. , 
wr i te ($ In-charge divergence avg. over last 12 orbits is $) , 
write (Div) ,nl , 

wri te($Avg. of cell pressures at EOD over last 12 orbits is $) , 
wr i te ( Eodpress ) , ri 1 , 

write($Avg. of cell pressures at EOC over last 12 orbits is $) , 
w r i t e ( E oc p r e s s ) , ri .1 , 

wr i te ($7emperature trend is $) ,write(79) , ril , 
tel 1...2 (72, 73, Div, 79, Eodpress , Eocpress ) . 

oil g w r k ( 7 .1. ,74, A h o a v g ) : - 

w r i t e ( $ E 0 D v o 1 1 a g e t r e n d .1 s $ ) , w r i i e ( T 1 ) , n 1 , 
write($E0D divergence trend is $) , write (T4) ,ril , 
write($Average amp-hour s-out. over last 12 orbits is $) , 
wr i te ( Ahoavg ) , n 1 , 
t e 1 .1. _ 3 ( 7 .1 , 7 4 , A h o a v g ) . 

t e .1 1 _ 3 (down, d o w n , A h o a v g ) : - 
Ahoavg < 5, 

wr i te ( $**Battery loads are probably too light — memory effects), ri 
write(S is indicated. $) . 


tell ...3 (strongly_down , strongly_dowri , Ahoavg) : - 
Ahoavg < 5, 

wr i te ( $**Bat tery is losing capacity due to memory effect. Try a$ 
n .1 , w r i te ( $ deepe r DoD f i r s t $ ) , n 1 . 

t e 1 1 ... 3(_,_, Ahoavg) 

Ahoavg > 14, 

wr i te ( $*+Ba ttery is overworked — may result in destructive damage 
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te 1 1 _3 ( Ah oavg ) : ~ 

Ahoavg >9, 

w r i t. e ( $**B at ter y workload seems heavy at $) •, 
w r j. t e ( A i'i o a v g ) , w r i 1 e (. $ A H u . $ ) , 
n 1 . 


3 C r 1 , T4 , Ahoavg ) : ~ 

( t 1 = down; T 1 = strongly.. down) , 

( T,j. ~ down; T 4 = st r ongly __down ) , 
Ahoavg > 6 , 

write (c&^Battery is losing capacity . 


Consider reconditioning - $ ) , nx . 


te 


write($See no 


need to change 


woe k 1 . oad . $ ) , n .1 . 


tel. 1_2 strong ly_up , Eodpress , Eoc press) : 

Eodpress > 100 , 

E o c p r e s s > 110 , 

w r i t e ( $ * * 0 v e r c h a r g i n g c a u s ed high pre s s u r e . 
wr i te (.$ circuitry, or change charge regime 


C h e c k c h a r g e 1 1 m i t. $ ) 
i turned i a te 1 y - $ ) , n i . 


, n i , 


t e 1 1 _2 v T r snd2 , ^ 

(Trend2 is stror.gly_up) , 

wr i te ($**Gheck for overcharging; voltage is nigri.'J/ 


2 (_ Eodpress , Eocpress ) : 

Eodpress > SO, 

Eocpress > 100, 

write ($+'+Charge rate may be too high 
n 1 , wr .i te ( S gas recombination lowet b 


Pressure high at EOC, thou 
p r essu r e d u r i n g d i sc h a r ge - $ > 


girt 


, n J. 


tell 2 (strongly _down , strongly.down 

"write ($**Battery undercharging may ^ be seriously ariei. 

r 1 1 , w r i t e ( $ c a p a c i. t. y f or wo r k . $ ) , n 1 . 


ing batteryt) , 


te .1 1 _2 ( down , down ) : . , 

write{$**Bat.tery appears to be under chat ged.S; ,n.L 


tel 1 _2 < T2 , T3 . 

(T2 is strongly _down ; T2 is down), 

( T3 is st.rongly_dowri; T3 is down), 

wr i te($**tJndercharging may be causing loss of voltage. $/ 


t e 1 1 _2 ( down , _ , D i v ) : 

~Di v < 0.3, 

write($**Undercharging may be causing lo 


of vol tag© . $ ) , n i - 


tell_2 (strongly_down , _ , Div 

Div <0.8, . 

wr ite ($*#Under charging may be causing set ious 


loss 


o f vo 1 1 age . $ ) , n 1 
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write <$See no need to change charging scheme .$), nl . 

tel 1 1 (strongly ...down strong ly_up ) re- 

write ($** Reconditioning advised to correct failing (capacity . $ J , ru - 


t e 1 1 _ 1 ( s t r o n g 1 y _ d o w n , s t r o n g .1. y _ d o wn,_,_ , „ ) - - 

wr i t. e ( $*+ Rec ondi t. i on soon before battery fails.®) ,ni . 

tell _1 ,strongly_up , st. rongly_up) 

wr i te ($ : + :: + : Strongly r .1 si. rig divergence suggests recondi tioni ng soon . ® } , 
n .1 . 

tel 1_1 < T1,T2, T3 , T4 , T5 ) : - 

( T1 is strongly _dowri ; T1 is down), 

(T2 is strongly _down ; T2 is down), 

(T3 is strorigly_up; T3 is up), 

( T 4 is strongly_up; T4 is up), 

( T5 is strongly_up; T5 is up), 

wri te($+*Reconditioning is indicated by all five trends . S ), ni. . 
te 1 1. _ 1 st rongly_up,_,_) 

wr i te ( $ : tt Recharge ratio indicates serious loss of efficiency ®) , 
nl,write($ suggest reconditioning soon . $ ) ,nl. 

tell ,T2,_,_,T5) : ~ 

( t 2 is down; i 2 is strongly_down ) , 

(15 is up; T5 is strongly .up) , 

wr i t e ( $** Rec onditioriing recommended. Too many cells are migrating®.) • 
nl, write ($ to low values . $), nl . 

tel .1...1 ( T .1. T 5 ) 

( T 1 i s d o w n ; T 1 .1 s s t rorig 1 y_down) , 

(T5 is up; T5 is strongly_up) , 

wr i te (SttCorisider reconditioning to correct poor performance . $) ,n l . 
tel. l_.l up, up) 

w r 1 1 e ( $ * c o n s i. d e r r e c ondi ti o ning t. o c o r r e c t d i v ergen c e . $ ) , n 1 . 
t.ell_.l up, :~ 

write ($ : + : *Consider reconditioning to correct rising recharge ratio.®) 
nl. 

tel 1_1 ( down , down . 

wr i te ( St -rBattery is losing capaci. ty ; consider recondi t;i.on ing . ® ! , nl 

t e 11 1 \ , , , , ) 

wri. te ($See no need to recondition at present®) , nl . 
dev (2.0,2.0,0.02,0.05,0.09,4.3). 


■■ B- li " k 



Showpak.prg contains the DECISION SUPPORT system, which gives 
the 12 plots which are available to the user upon selecting 
Plots and Graphs from the Main Menu. The graphics] primitives 
from Grafpak.prg are used in the plots. Showpak contains a few 
supporting data handling primitives and the titles, horizontal 
and vertical captions used in the plots. 

12 Plots in the DECISION SUPPORT package: 


1) 

EOD battery 

low plotted 

26.4V. 



voltage 

high plotted 

28.6V 




scale 

. 2 V 


2) 

In-charge high 

low plotted 

32.0V 



cell voltage 

high plotted 

34.2V 




scale 

. 2 V 


3 ) 

recharge ratio 

low plotted 

1.010 




high plotted 

1.055 




scale 

.005 


4) 

EOD divergence 

low plotted 

1.05V 




high plotted 

1.55 V 




scale 

. 05V 

(cell) 

5) 

In-charge 

low plotted 

1.35V. 



divergence 

high plotted 

1.55V. 




scale 

. 02V. 

(cell) 

6) 

EOD cell V. 

low plotted 

1.05V. 



divergence 

high plotted 

1.55V. 




scale 

. 05V. 

(cell) 

7 ) 

In-charge cell 

low plotted 

1 . 36V. 



divergence 

high plotted 

1 . 56V. 




scale 

. 02V. 

(cell) 

8) 

Temperature 

low plotted 

-4C 



this orbit 

high plotted 

11C 




scale 

1C 


9) 

Temp. avg. 

low plotted 

A ** 

-HO 




high plotted 

lie 




scale 

1C 


10) 

Cell Pressure 

low plotted 

30 



B' l 



EOD and EOC 


high plotted 140 

scale 10 units 


11) Time on 

trickle charge 


low plotted 
high plotted 
scale 


8 

28 minutes 
2 minutes 


12) Reconditioning 
current 


low plotted 0A. 
high plotted 14 A 
scale 1 A. 


The above numbers can easily be manipulated to change range or 
scale of the plots, but must remain integers. These changes can be 
made in Showpak.prg by editing the correct ’show’ data structure, 
’show’ has the following format: 


1. 

N 

- 

the number of the plot (1 to 

12) 

2. 

Vertcap 

- 

vertical caption 


3. 

Horizcap 

- 

horizontal caption 


4. 

Title 

- 

graph title 


5. 

H 

- 

heigth of graph 


6. 

W 

- 

width of scale 


r* 

/ . 

Base 

- 

lower bound for data points 

plotted 

8. 

Top 

- 

upper bound for data points 

plotted 

9. 

Hscale 

- 

horizontal scale of values 


10. 

Vscale 

- 

vertical scale of values 


11 . 

S tart 

- 

start point of graph 



Changing the number of orbits plotted from 12 is more difficult, and 
will require chasing down many ll’s and 12’s, as well as rewriting the 
file handler portion of the code. 

Colors are 

1 = blue 

2 = green 

3 = light blue 

4 = red 

5 = pink 

6 = gold 

7 = white 

8 = gray 

Symbols for the graphs are referrenced by their ASCII code. 


Colors are set in predicates 

"wa ( n , col or ) " , where n is the number 

of characters of that color, and 

color is one of the numbers to 

the lef t . 


THE PREDICATES AVAILABLE IN SHOWPAK.PRG 

show_view ( N , Bat , Orb i t ) — N specifies which plot is done, and Bat 

specifies which battery is displayed. Orbit gives the current 
orbit number. The data structure ’show’, discussed above, is 
called to get all the captions and values needed for plot N. 
show_view/3 calls get_list/3 to read the appropriate data file. 


V2-D 



Then graphplus/6 and plot/8 (in Grafpak.prg) are called to plot 
the display. Besides the information in ’show’, only the color 
and symbol need to be passed. The present orbit is then written 
to the plot, if applicable. show_view/3 is called by execute/2, 
the Graphics Menu. 

show_view/3 — There are four cases; plot 10, plot 12, plot 4 and 5, 
and the rest. Plot 10 displays 2 lines, plot 12 has the orbit no. 
preceeding its data, plot 4 and 5 display 3 lines each. For plots 
4,5 and 10, breaklist/4 is called to break the data into 
lists, each one to be plotted on the same graph. plot/8 is then 
called for each list with a different color and symbol. 

write_orbit/0 — writes ’Orbit* on the horizontal axis. Called by 
show_v iew/3 . 

wr i t e_orb i t /I — writes ’Orbit’ and the last orbit on the horizontal 
axis. Called by show_view/3. 

show/ll — a data structure containing the captions and values needed 
to display the graphs. Itemized above. 
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% Showpak* — contains the decision support system, i . e . , show view/ 
% build the 12 plots which support analysis in Status arid Advice. 
Z See showpak . doc fro further documentation. 

sh o w _ v i e w ( 1 0 , B a t , 0 r b i t ) : - 
ge t ... 1 i s t ( B a t , 1 0,Li s t. ) , 

cal 1 ( show ( 10 , Vcap , Heap , T.i. tie , H , W , Base , T op , Hscale , V scale , Star t 

g r a p h p 1 u s ( H , W , V c ap , He a p , T i 1 1 e , B a t) , 

b r e a k _ 1 i s t ( L. ist. , L i stl ,Li s 1 2 ,23) , 

plot. (Listl , Hscale , Vscale , Base , Top , 3 , Start , 254 ) , 

plot (List2 , Hscale , Vscale , Base, Top, 4 , Start. , 4 ) , 

w r i t e _ o r b i t ( 0 r b it),!. 

s h o w _ v .1 e w (12, B a t , 0 r b .1 1 ) 

ge t._ 1 i s t ( B a t, 12, [R, List]), 

c a .1. 1 ( s h a w ( 1 2 , V c a p , H c a p,Ti 1 1 e , H , W , B a s e , T o p,H s c a 1 e , V s c a 1 e , S t a r t 
graphplus(H, W, Vcap , Heap , Title, Bat) , 
plot (List , Hscale , Vscale , Base , Top , 4 , Star t , 254) , 
write_orbit(R) . 


show_view( V, Bat. , Orbi t. ) 

( V - = 4 ; V == 5) , 
g e t _ 1 i s t ( B a t , V , L i s t ) , 

call (show ( V , Vcap , Heap , T i tie ,H,W, Base , Top , Hscale , Vscale , Star t ) 

g r aph p 1 us (H, W, Vc ap , He ap , Tit 1 e , B a t ) , 

b r e a k _ 1 i. s t i L i s t , L i s 1 1 , L .1 s 1 2 ,12), 

b r e a k _ 1 i. s t ( L i s 1 2 , L i s 1 3 , L i s 1 4 , 1 2 ) , 

plot ( List 1 , Hscale , Vscale , Base , Top , 3 , Star t , 254 ) , 

P 1 o t (Lists, Hscal e , V s c a 1 e , B a s e , T o p,4, S t a r t ,4) , 
plot (L.ist.4 , Hscale , Vscale , Base , T op , 7 , Star t. , 42 ) , 
write_orbit, ! . 


::;how_view ( V , Bat , Orbi t. ) : - 

gei_ 1 ist (Bat , V , List ) , 

call (show ( V , Vcap , Heap , Title , H, W, Base , Top , Hscale , Vscale , Star t ) 
graphplus (H , W , Vcap , Heap , Ti tie, Bat) , 

P 1 o t (List, Hsc ale, Vscale, Base, Top, 4, S ta r t , 254 ) , 

< ( ( V = = ,9 ; V == 6; V = = 7 ) , wr .1 te_orbi t ( Orbi t ) ) ; 
wri te_orbit) , ! . 


write_orbit 

tmove (24,3) , 
wa (5,2) , 
write($0r bit.$) - 
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Of poor Qutm 

w r i t e _ o r b i t ( 0 r b i t ) : - 

f.niove ( 24 , 0 ) , 
wa ( 10 , 7) , 
w r i. t e ( $ 0 r b i t $ ) , 
wr .1 te ( Orbi t ) . 


ST f OW ( 1 , [ b ' - 6 1 b O . 4 , b '.'J - 2 -1 2b - 0 , 2 / .. ' - O 7 2 / - 4 , .--V / - U i 

He a pi , ’ EOD VOLTAGE FOR LAST 12 ORB ITS ’,13,65,26.4,28 

SO O W ( 2 7 [ 4' _ b n v > . 0 7 JO _ b 7 bb - O , Ob - 4 , bb .il 7 bb • 0 , . O 7 04 - 6 7 

Heap!,’ HIGH VOLTAGE DURING CHARGE FOR LAST 12 ORBITS 
.1 3 7 6-5 7 32 . 0 , 34 . 2 , 5 , 0 . 2 7 .1 1 ) . 



show (3 , [1.055,1.050,1 .-045,1.040,1 . 035 , 1 . 030 , 1 . 025 , 1 . 020 ,1.015,1 . 010] , 
Heap!,’ RECHARGE RATIO FOR LAST .12 ORB ITS ’,11, 65, 1.010,1. 055 ,5,0. 005 ,11) 

show ( 4 , [ 1 . 55 , 1 . 50 , 1 . 45 ,1.40,1.35,1.30,1.25,1.20,1.15,1.10,1. 05] , Heap! , 

7 EOD HIGH (blue) - LOW (red) - AVG(wh) FOR LAST 12 ORBITS’, 

1 2 , 65 , 1 . 05 , 1 . 55 , 5 , 0 . 05 , 11 ) . 


show ( 5 ,[1.55,1.53,1.51,1.49,1.47,1.45, 
’ HIGH (blue) - LOW (red) - AVG(wh) - 
12 , 65 , 1 . 35 , 1 . 55 , 5 , 0 . 02 , 11 ). 


1.43,1.41,1.39,1.37 
DURING CHARGE LAST 


.1. . 35] , He a pi , 
12 ORBITS’, 


show ( 6 , [ 1 . 55 ,1.50,1.45,1.40,1.35, i. 30, 1.25, 1.20, 1.15,1. 1 0 , 1 . 05 1 , 

T , 1 ’ 1 ’ j_ i ’ ’ 16 1 ’ ■:>! 7 1 p f, 1 ’ J 

’ 23 EOD CELL ’ VOLTAGES FOR LAST ’ ORBIT ’ ,13, 62 , 1 . 05 , 1 . 55 , 2 , 0 . 05 , 1 1 ) 


1 C 


show ( 7 , [ 1 . 56 , 1 . 54 , 1 . 52 , 1 . 50 , 1 . 48 , 1 . 46 ,1.44,1.42,1.40,1. 38 , 1 - 36] , 
f 1 7 7 6 7 7 1 i 7 7 1 6 7 7 2 j 7 7 26 7 7 ] , 

7 ’ 23’ HIGH CHARGE CELL FOR LAST ORBIT ’ , 13,62, 1 . 36 , 1 .56 , 2, 0 .02 , 11 ) . 

sh o w ( 8 , [ 7 1 1 C 7 , 7 1 OC 7 , 7 9C 7 , 7 8 G 7 , 7 7C 7 , 7 6 C 7 , 7 5C 7 , 7 4C . 7 , 7 3C 7 , 7 2C 7 , 7 1 C 7 , 7 OC 7 , 7 

7 -20 7 , 7 -3C 7 , 7 — 4C 7 ] , [0, 10,20,30,40,50,60,70,80,90, 100] , 

7 TEMPERATURE IN LAST ORBIT, EACH 2 MIN 7 , 17 , 60, -4 , 11 , 1 , 1 , 12) . 

show (9, [ 7 11C 7 , 7 IOC 7 , 7 9C 7 , 7 8 C 7 , 7 7C 7 ,’6C 7 , 7 5C 7 , 7 4C 7 , 7 3C 7 , 7 2C 7 , 7 1C 7 , 7 DC 7 , ’-1C 
7 -2C 7 , 7 — 3C 7 , 7 -4C 7 ] , Hcapl , 

7 AVG. TEMPERATURE IN EACH ORBIT, LAST 12 ORBITS’ , 17,65, -4, 11 , 5, 1 , 11 ) . 

Sh ow ( 1 0 , T 1 40 ,130,1 20 ,110,1 00 , 90 , 80 , 70 , 60 , 50 ,40,30], 

[l , 7 ’, 6 , 7 ’, 11 , 7 ’,16, 7 ’, 21, 7 7 , 26 , 7 7 ], 

7 CELL PRESSURE EOD (red), EOC (blue), DURING LAST ORBIT 7 , 

13,62,30, 140,2, 10,11) . 


show (11, [28 , 26 , 24 , 22 , 20 , 18 ,16,14,12, 10 , 8] , Hcapl , 

7 TIME ON TRICKLE CHARGE, LAST 12 ORBITS ’ , 12, 65, 8, 28, 5 , 2 , 11 ) . 
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h ow ( 
0 , 1 0 


1 2 , [ 1 4 , 1 3 , 1 2 , 1 1 ,10,9,8,7,6,5,4,3,2,1,0] , 

, 20 , 30 , 40 , 50 , 60 , 70 , 80 , 90 , 1 00] , 

UR RENT DURING CAP. TEST/RECOND, EACH 2-M.IN, 0 => NO RECOND 
,0, 14,1 , .1 , 12) . 
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Grafpak.prg contains graphic primitives and supporting predicates 
for Showpak.prg to use in giving the user plots and charts to aid 
in decision support. The procedures are written to be 
general -purpose so that plots may be altered easily. New plots may 
be developed by the programmer using these primitive operations. 

Arity/ Prolog graphical primitives 


Parameters in plot are 


(List, 
Hscale , 
Vscale , 
Base , 
Top, 
Color, 
Coord , 


Obj) 


list of values to be plotted 

distance horizontally between plotted points 

vertical scale of values 

lowest value plotted 

highest value plotted 

color of the plotted points 

the horizontal positioning factor, generally starts 
with 11 (columns from left-hand side of screen), 
increments by "num" with each point (Start) 

ASCII code for symbol to be plotted 


Parameters in graphplus are 


(H, 


W, 

Vcap , 
Heap , 
Title) 


Rows in height (1 more than number 
of values wanted) 

Columns in width (5 more than the spaces 

needed to plot values — see "coord" 
Vertical caption 
Horizontal caption 
Title written above plot 


Built-in Graphics Primitives: 

tmove(x,y) — moves the cursor to row x (rows are numbered 0 to 

24 beginning at the top) and column y (columns are numbered 0 to 
79 beginning at left). 

wc(n,char) — prints n copies of char (given as an ASCII code 
or within single quotes) starting at the cursor. 

wa(n, color) — prints the next n characters in color (given as a number). 

1 = blue 2 = green 3 = It. blue 4 = red 

5 = pink 6 = gold 7 = white 8 = gray 

Characters can also be made flashing or inverse, in any of the colors, 
using this command. To see the attributes that are available, with 
the Arity interpreter running, type "[’ table. ari ’ J . " and when Arity 
returns "yes" type "table. " 




THE PREDICATES IN GRAFPAK.PRG 


graphplus(H,W, Vcap, Heap, Tit 1 e , Bat ) — creates the basic graph with 

height H and width W, vertical caption Vcap (usually the scale), 
horizontal caption Heap, and title Title including the battery 
number, Bat. H should be chosen to be 1 greater than the number of 
rows (which is the number of possible, values that will be plotted). 

W should be chosen to be 5 wider than the spaces needed to plot the 
values. See plot/8 for more information. 

graphplus/6 calls graph ( H, W , Heap ) to build the skeleton along 
with the horizontal caption, then calls write_vert ( Vcap ) , and 
wr i t e_ t i t le ( Ti t le , R , Bat ) to fill in the captions. write_header is 
also called to write information about missing data and data out of 
range at the the top of the graph. 

graphplus/6 is called by show_view/3 in Showpak.prg, with the proper 
parameters of size and captions. 

graph ( H , W , Heap) — builds the skeleton graph of height H and width W, using 
vertical/2 for the vertical graphics characters making the left-hand 
side of the graph, and horizontal/2 for the lower edge. graph/2 is 
called by graphplus/6. 

vert i cal ( A , W2 ) --uses graphics characters 180 (ticks) for the vertical axis 
and 250 for the dots across the inside of the graph. 

hor izon t a.l ( B , L i s t ) --uses graphics characters 196 and 194 for the bottom edg 
of the graph. The tick marks are placed every five spaces. The 
elements of the list, containing the horizontal captions, are alined 
with the ticks but one line under. 

wr ite_vert ( L , P ) --writes the list L along the left-hand side of 

the graph. Calls write_scale(L) which writes the list recursively. 

write_scale(L) — writes the vertical caption recursively from the 
list L. The empty list ([]) halts the recursion. 

wri te_t i t le ( S , P , Bat ) — writes the Title along with the Bat, battery 
number, above the graph. S is the Title and P is the row 
position . 

pi ot ( Li st , Hscale , Vsca le , Base , Top , Col or , Coord , Ob j ) — Plots the values in 

list List on the basic graph produced by graphplus/6. Hscale is the 
horizontal distance (in columns of the screen) between the points 
plotted. Vscale is the vertical scale of the graph. Base is 
the lowest value plotted, and Top is the highest value 
plotted. An attempt to plot any value outside this range will 
result in out of range values being plotted at the top of the graph 
and the actual values being written, in order at the top of the 



screen. Top minus Base divided by Vscale will equal the number of 
rows upon which the graph is plotted. Color is the color of the 
plotted point. Coord is the horizontal positioning factor, as the 
plot moves across the screen from left to right. Coord ordinarily 
starts at column 11, and increments by either 1 or 5. Obj is the 
ASCII code for the symbol which will be plotted for each data point. 

plot/8 calls plot_point/8 to recursively plot the points from 
the list L. When plot_point/8 is done, plot/8 changes the color 
back to white and moves the cursor to the top left-hand corner. 

plot/8 is called by show_view/3 in Showpak.prg with the data 
in the List and specific parameters for the plotting. 

plot_point([Head;Tail],Hscale, Vscale, Base, Top, Coord, Color, Obj) — plot_point/8 
recursively plots the list of points, divided into Head and Tail, 
then calls itself on the Tail, having updated Coord for the next 
point. The empty list halts the recursion. Hscale, Vscale, Base, 

Top, Coord, Color and Obj are the same as in plot/8. 

plot_point/8 calls place/7 to find the correct place for the 
plotted point, and plot it there in the correct color. 

p lace ( Head , Vscale , Base , Top , Coord , Color , Ob j ) — figures the correct place 

for the plotted point, moves the cursor there, and plots the point. 

The value for missing data causes an asterisk in gold 
to be plotted on the low boundary line of the graph. place/7 
truncates values as it plots them. When values are out of range, 
place/7 plots a symbol at the top boundary of the graph and writes 
the value at the top of the screen. 

place/7 is called by pi ot_point/8 . 
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ts 

OF POOR QUALITY 


% grafpak.prg — contains graphic primitives and supporting predicates. 
% graphp.lus builds a graph skeleton given height, width and captions.. 
% plot, draws the data points on the graph. 

% S e e G r a f p a k . d o c f o r f u r t h e r d o c u m e n t a t i o n . 


graphplus (H, W, Vcap , Heap , Ti tie , Bat. ) : - 
gr aph ( H , W , He ap ) , 

U is 24 - H, 
w r i t e _ v e r t. ( V c a p , Q ) , 

R is Q - 2, 

writ.e_title(Tit.le, R,Bat.) , 
wri te_beader , 
tmove (1,1) , 

R1 is R + 1, 
ctr_set(2, R1 ) , 
ctr_set( Id , 0) . 


w r i t e _ h e a d er : - 
tmove (2,0) , 


wa ( SO , 6 ) 
tmove ( 2 , 
wc ( 1 , 254 
tmove ( 2 , 
wr i te ( $'= 


, wc ( . 1 . , 236 ) , 

2) , write <$= MISSING DATA, $) , 

) , tge t ( .... , X ) , XI i s X + 2 , 

XI ) , 

VALUE OUT OF RANGE, LISTED ON NEXT 


LINES) . 


a r a oh ( H , W , He a p ) : - 

els, 

HI is 24 - H, 

W1 is W - 1, 
vertical (HI ,W1) , 
horizontal (W, Heap) , 
tmove ( 23 , 5 ) , 
wa (. 1 , 2 ) , 
wc(l , 192) . 

ve r t i e a 1 ( A , W2 ) :~ 
ctr_set(0,A) , 
repeat. , 

etr_inc (0 , Al) , 
tmove (A1 , 5) , 
wa (1,2) , 
wc ( 1 , 180) , 

•Hmrti/rt ( A 1 \ 

L>iii w vc: \ nx -t *-.* / 7 

wa ( W2 ,2) , 
wc ( W2 , 250 ) , 

Al == 22. 
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h o r izo n t a 1 ( B , Li s t ) : - 

B 1 is (B - 2)// 5, 
c tr_set (0 , 1 ) , 
t move (2 3, 6) , 
wa ( 1 , 2 ) , 
wc ( 1 , 196) , 
tmove ( 23 , 7 ) , 
repeat , 

i j. s t i t e r a t e ( H , L .1 s t ) , 
c t r .... i ri c .... 1 ( 0 , C ) , 
wa ( 4 , 2 ) , 
wc ( 4 , 196 ) , 

D is 5 * (C + 1) + 1, 

tmove (23, D) , 

wa (1,2) , 

wc ( 1,194) , 

tmove (24 , D) , 

wa(l 0,2) , 

write (H) , 

E is D + 1, 
tmove (23, E) , 

0 -= Bl, 

F is (B - 2) mod 5, 
wa(l 0,2) , 

1 f t h e ri ( F 8 > 0 , w c ( F , 1 9 6 ) ) . 


: tr_inc 1 ( Cnt , N) 


ctr_inc (Cnt, N) , ! . 


1 ist_i ter ate (H , [HIT] ) - 
1 i s t _ i t e r a t e (H, [ _ I T ] ) : • 
1 i s t _. i t e r a t e ( H , T ) . 
1 i s t ... i t e r a t e (’’,□) - 

wr i te_ver t ( L , F ! ) : - 
tmove ( P , 0 ) , 
write_scale(L) . 

w r i t.e_sc a 1 e ( [ j ) - 
wr i te_scale ( [HIT]) : - 
wa ( 5 , 2 ) , 
write(H), ril, 
w r i te_.sc a 1 e ( T ) . 


w r i t e _ t i 1. 1 e ( S , P , B a t ) : - 
tmove ( P , 6 ) , 
wa ( 80 ,7) , 

wr i. te ( $B ATTERY $) , write (Bat) , 
write(S) . 
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lot (List , H scale , Vscale, Base, Top , Color , Coord, Ob j ) : - 

plot_poin t (Li s t , H s c a 1 e , V s c a le,Base , T o p , C o o rd, Color 
tmove ( 1 , 1 ), 


At a (. 1 


Ob i ) , 


p 1 o t ... p o i n t ( r ] , ... ... ) - 

p 1, o t po i n t ( [. Head ! Tai i J , Hsca le , Vsc a Is , Base , i op , Uoor d , Col or , od j ) 
place ( Head , Vscale , Base , T op , Coord , Color , Ob j ) , 

Coordl is (Coord + Hscaie), 

p j. o t .... p o i n t (Tail , H s c ale, V s c a 1 e , B a s e , T op , C o ordl, C o 1 o r , 0 b j ) . 

P .1. a c e ( -- 9 9 9 9 , ... , .... ,_,C o ord,_,_) 
tmove ( 23 , Coord ) , 
wa (. 1 , 6 ) , 
wc ( 1 , 236 ) . 


place (Head , Vscale , Base , Top , Coord , Color , Ob j ) 
(Head < Base; Head > Top) , 
ctr_is(2, R.1. ) , 

.tmove (R.1, Coord), 
wa ( 1 , 6 ) , 
we(l,Obj), 

Ctr_is ( 10 , R2) , 
tmove ( 3 , R2 ) , 
wa ( 5 , 6 ) , 

w r i t e ( H e a d ) , w r 1 1 e ( $ $ ) , 

t g e t ( ... , R 3 ) , c t r ... s e t ( 1 0 , R 3 > . 


pi ace ( Head , Vscale , Base , T op , Coord , Col or , Ob j ) : - 
Q .1 is (Head -- Base) /Vscale, 

Q is round (Q 1,0), 

Level Is 22 - integer (Q) , 
t m o v e ( L e v e 1 , C o o r d ) , 
wa ( 1 , Color ) , 
wc ( 1 , Ob j ) . 


Splace (Head, Vscale , Base , Top , Coord , Color 
% Q1 is (Head - Base ) /Vscale , 

% Q is round (01 , 0 ) , 

% !.. e v e 1 1 s 22 - i ri t e g e r ( Q ) , 

% tmove (Level , Coord) , 

."o wad. Color), 


wc ( 1 , Ob j ) . 



Obj ) 
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Utility. prg is a collection of general Prolog predicates which are used 
in one or more of the routines comprising NICBES. 


Data Files in the Expert System: 

All data files are loaded into PROLOG as facts in the following 
f ormat : 

showf ( N , Lis t ) . N = 1 to 13 for showf(N).dat 
cur f ( N , List ) . N = 1 to 3 for curf(N).dat 

f ault ([ Flag , Type ]) . for fault.dat 


List for show files contains 6 sublists, one for each battery, 
seperated by commas. For example, in PROLOG, showfl.dat is loaded as: 


showf ( 1 , [ [27.0,26.4,27.0,26.4,-9999,27.0,26.4,27.0,26.4,26.9,26.8,27.2 
[26.8,27.0,26.8,27.0,26.8,26.6,27.0,26.8,27.0,27.0,26.4,27.0], 
[27.2,27.2,27.2,27.0,26.8,26.8,26.8,26.8,26.8,27.1,26.5,26.0], 
[27.5,26.9,26.9,26.9,26.9,26.8,26.8,27.0,26.8,27.0,27.0,26.7], 
[27.4,27.0,27.1,27.0,27.0,27.3,27.2,27.0,26.9,26.8,26.8,26.8], 
[27.0,26.4,27.0,26.4,26.4,27.0,26.4,27.0,27.0,281,27.2,27.5] ] ) . 


Actual data points are reals with 5 decimal places. -9999.0 
represents missing data points. There are 12 columns, one for each 
orbit in chronological order from oldest to latest. 


An example of the current data files is curf2.dat: 

curf ( 2 , [ Phase , Day_min , [ 13 SPA Currents], [3 Bus Currents], 
[6 Average Battery Temperatures]]). 


Built-in ARITY PROLOG Functions: 


nl — gives a new line (like carriage return). 


els — clears the screen. 


Many of the functions use the ARITY PROLOG counters. There are 32 counters 
(0 to 31) which are accessible from any procedure in any file. The 
operations on a counter are: 


ctr_set ( n , m) 

ctr_is(n,A) 

ctr_inc(n,A) 

ctr_dec ( n , A) 


set counter n to value m 

A is instantiated to the value of counter n 
counter n is incremented, and A is instantiated 
to the previous value 

counter n is decremented, and A is instantiated 
to the previous value 


Underscores (_) appearing in place of arguments signifies an anonymous 
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var iab le--we do not care what the value is and will not be using 
it at this time. This'saves the system work. 

ifthen and ifthenelse use the following format: 
ifthen(X,Y) => if X then Y. 
i f t hen e 1 s e ( X , Y , Z ) => if X then Y else Z. 

var(X) — succeeds if X is a variable. 

integer(X) — converts X to an integer. 

round(X,N) — rounds X to N decimal places. 

read_ 1 i ne ( H , S t r ing ) — reads string from file H. When H == 0, reads string 
from keyboard. 


PREDICATES AVAILABLE IN UTILITY. PRG 

get_ 1 i s t ( Bat , N , Lis t ) — N tells which show file. get_list calls the particular 
’showf’ fact which was loaded in eval_flag(0) (start. prg) and which is 
described above. find_nth is called next to find the data list 
corresponding to Bat, the battery of interest. This list is returned to 
the variable List. get_list/3 is called from complete/1, status. prg, 
advice. prg and showpak.prg. 

find_nth/4 — given a list, find_nth/4 finds the Nth element of the list. 

It is called by get_list/3. 

append ( Li s t 1 , Li s t 2 , Li s t3 ) — appends List2 to Listl, with the 

result in List3. append/3 calls itself recursively, with the 
empty list the halting condition. 

break_ 1 is t ( Lis t , Lis t 1 , Lis t2 , N) — breaks List into two lists. The first N 
items are put in Listl, and the remainder into List2. After setting 
counter 1 with N, break_list/4 calls steal_item to recursively build 
the lists. break_list/4 is called by show_view/3 ( graf pak . pr g ) , 
disp_trend/3 (advice. prg) and find_div/5 ( status . prg) . 

s teal_item ([ H ! T] , List 1 , Lis t2 ) — recursively traverses the tail of its first 

argument, decrementing the counter, until the counter is 0, or the first 
argument is the empty list. Then List2 is set to the remaining portion o 
the original list (or the empty list), and Listl is built as steal_item/ 
climbs out of the recursion. steal_item/3 is called by break_ 1 i s t /4 . 

f ind_avg ( L is t , Avg ) --f inds the average of a list of numbers. Counter 10 is 
set to count the number of items in the list. find_sum is called to 
sum them, then the average is found. find_avg is called by get_data/5 
and find_div/5 (status . prg) . 

find_sum( [H ! T] , Sum) — recursively finds the sum of a list of numbers. The 
head of the list (H) is checked for missing data (-9999), counter 10 
is incremented, then find_sum/2 is called on the tail (T) of the list. 
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The empty list halts the recursion, then the sum is found as find_sum 
climbs out of the recursion. A missing data item is left out of the sum 
and a message is printed. find_sum/2 is called from find_avg/2, find_w2 
(advice. prg) and 1 oadbank_f ai 1 /0 ( f aul t d . prg ) . 

s equen ce ( Orb i t ) — creates the horizontal axis for those graphs using the 12 

latest orbits. It calls seq/3 to make this list and hor iz_place,' 1 to pu 
it in the data structure show/ll. It is called by eva l_f 1 ag'( 1 ) 

( s tart . prg) . 

seq ( H, T , [H ! Z ] ) — Builds the list of last 12 orbits recursively from the latest 
orbit, T and the first orbit in this series, H. Called by sequence/1. 

horiz_place( List ) — searches the data structures show/ll to find ones with 

variable Hcapl, the horizontal axis. It retracts these and then asserts 
them to the data base with the List in -the place of the variable. 

Called by sequence/1. 

cut_retract/0 — retract fact and cut. Called by hor iz_pl ace/ 1 . 

check(X,N) — used for checking the user input to menus. If the user 
selection X is >0 and <= N (upper bound), the screen is 
cleared and program continues. If not, a message is written to 
the user saying ’YOUR CHOICE IS OUT OF RANGE!’. The user can 
then input a proper selection. Called by Status, Advice and Graphics 
menus in start. prg. 

check2(X,N) — used for checking the user input to Main Menu in start. prg. 

If the user selection X is >0 and either X <= N (upper bound) or equal 
to 4, the screen is cleared and program continues. If not, a message 
is written to the user saying ’YOUR CHOICE IS OUT OF RANGE!’. The user 
can then input a proper selection. 

checkl(X) — used for checking the user yes/no responses to menus. If the user 
input is either yes or no, the screen is cleared and program continues. 
If not, a message is written saying ’ANSWER IS OUT OF RANGE! ’ . The 
user can then re-enter a correct reply. Called by eval_flag(l) (start. pr 

complete(M) — calls get_list/3 to look at show file 1. missing/2 returns the 
number of missing data points in show file 1. If there are more than 
4 missing values, M is set to 1 so only Graphics may be called from the 
Main Menu. Else M is set to 4 and all options from the Main Menu will 
be valid. Called from eval_flag(0) (start. prg). 

missing( List , N-) — returns the number of missing data points, N from the List. 
Called by complete/1. 

reader (Atos) — reads input from the keyboard and returns it as an atom. 

Called by all menus in start. prg. 

one_assert ( X) — if X is already in the data base, no action. Else 
the fact X is asserted to the data base. 
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Z U B I i 1 ity contains general functions which are used in one or more 
ot the routines comprising NIOBES. 

ge t .... 1 i s t ( B a t , N , L i s t ) : - 
c al 1 ( show f ( N , I... ) ) , 
f i r i d ... n t. h ( L , i , B a t , L .1 s t ) . 

f i ri d ._ n t h ( [ H IT] , N , N , H ) . 

f i n d n t h ( [ H I T ] , N 1 , N , X ) : - 
M is N1 + .1 , 
f ind_n th ( T , M , N , X ) _ 

br eak_i 1st (List , L.ist.1 , List.2 , N) 
ctr_set( 1 ,N) , 

s t e a I _ i t e m ( L i s t ,Lis 1 1 , L i s 1 2 ) . 

steai._.i tern ( [j , List! , List.2) -- 
List l = □, 

List-2 = [] . 

steal,. i tern (List , List! , List.2) - 
ctr is(l,M), 

M == 0, 

List! = [], 

List2 = List. 


s t ® a 1 - i t e m ( [ H I T ,] , L i s 1 1 , L. i s 1 2) : - 
c tr_dec ( .1 , ._) , 
s t e a 1 _ i t e m (T Temp, Li s 1 2 ) , 

List.l = [HI Temp]. 

f i n d ._ a v g ( L i s t , A v g ) : ~ 
c tr_set (10,0) , 
f i n d_sum (List, Sum ) , 

ctr_is (10, Num) , 

Avgi is Sum/ Num, 

A vg is r oun d (Avgl , 4 ) . 

f ind... sum ( [] ,Sum) 

Sum is 0.0. 

f'ind_.sum < [h I T] , Sum ) : - 
c tr_inc ( 10, _ ) , 
tirid_sum(T, Temp) , 

i f thenelse ( H == -9999 , ( Sum i s T emp , c t r _dec ( 1 0 , _ ) ) , Sum 


H + Temp) . 
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e q u e n c e ( 0 r b i. 1 ) - ~ 

H is Orbit - 11, 
Orb is Orbit + I, 
seq ( H , 0 r b , L .1 s t. ) , 
h o r i z ... p 1 a c e ( L i s t ) . 


? Q ( I 


) . 


;;eq ( H , T 
HI 


[H S Z] ) :~ 
is H + 1, 


seq (HI , T, 2 ) - 
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' "call (show(N, Vcap, Heap, Title, H,W, Base, Top, Hscale , Vscale , otat t.) ) , 

cut. * retract < show (N, Vcap , Heap , Title , H , W , Base , Top , Hscale , Vscale,bt 
asser ta (show ( N , Vcap , List , Title, H, W, Base , Top , Hscale , Vscalc , tat 

fail - 


h o r i z _ p 1 a c e (. ) 
c u t ... r e t r a c t (X) 
check. (X,N) :~ 


X > 0, 

X =< N, 


! 

retract(X) , ! - 


check (X,N) 

c 1 s , ri 1 , n 1 , 
wr.ite($ 

! , fail. 


YOUR CHOICE IS OUT OF RANKE ! S) , nl,rix. 


check2(X,M) :~ 

X > 0, 

(X — < N; X == 4), 


check2 (X , N) :~ 

c Is, ri 1 , ri 1 , 
write ($ 

! , fail . 


YOUR CHOICE IS OUT OF RANGE! $) , nl ,ni. 


check! (X) :~ 

(X == ’yes’ ; X ---- ’no’),!. 


c 1 s , r 1 1 , n 1 , 

write ($ ANSWER IS 

nl ,ril , ! , fail . 


OUT OF RANGE ! $ ) , 
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c o m p 1 e t e ( M ) : — 

get_list.(l , .1 , L i s t ) , 
c tr _ set ( 1 , 0 ) , 
missing (List , N ) , 
i f t h e n e 1 s e < N < 5 , M i s 4 , 

(write($ Not enough orbital 
n 1 , wr i te ( $ Check data files 
M i s .1 ) ) . 


d a t. a f o r N I C B t. S a n a. 1 y s i s 
Y o i .t may ru n G raph i c s . $ 


\ 


m i ss i n g ( [] ,N) : - c t r _ i s ( 1 , N ) . 

m i ss i n g ( [ H IT] ,N) : - 

i f t h en ( H = = - 9999 , c t r i n c ( 1 , _) ) , 
missing (T , N) . 

reader (Atom) 

r e a d _ 1 i. n e ( 0 , S t r i n g ) , 

( i n t _ te x t. ( A t ont , S t r i n g ) 
a t o »n _ s t r i n g ( A t o m , S t r i n g ) ) , ! . 

a n e _ a s ser t ( X ) : ~ 

X - 



assert (X) 
asserta(X) . 
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TEST PROCEDURES 



,y 


\ 


1.0 Power on the IBM-PC AT and Printer. 

2.0 Run the Data-Handler . 

A. Data-Handler is in directory USR. 

B. Enter ’data_hdl’ to execute Data-Handler. 

C. Check telemetry flow across RS232. 

D. Messages printed to screen: 

"Starting first full orbit", 

ORBIT = N" - every 96 minutes. 

E. If a fault is detected, message printed to screen: 

"Fault detected; exiting!". 

F. If no communications being received, message printed to 

screen; "No communication for 3 minutes; exiting!". 

G. If 5 consecutive incomplete telemetry runs, message printed 

to screen; "Received 5 consecutive incomplete telemetry 
runs; exiting!". 

G. When ready to quit Data-Handler enter Message printed 

"Interrupt caught; exiting!". 

3. Verify and Archive Data-Handler Output. 

A. View showf#.dat (1 to 1 3 ) , curf # . dat (1 to 3) and fauit.dat 

for correct format. 

B. Archive data files under directory DATAFILES, 

subdirectory ’DATE’ where DATE is current date. 

4. Run Expert System. 

A. Expert System is in directory NICBES. 

B. Enter ’api’ to execute Expert System. 

C. View Fault Diagnosis. 

a. Print Fault Diagnosis Report. 

D. View Main Menu and make all selections. 

E. Select battery. 

F. View Graphics Menu and make all selections for one battery, 
a. Print one Plot. 

G. View Status for all batteries, 
a. Print Status Report for one battery. 

H. View Advice Menu and make all selections for one battery, 
a. Print Advice Report for one battery. 

I. Quit Expert System. 

Verify Expert System Output. 

A. Review printed Reports and Plots. 
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